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Final  Report 

The  Major  goal  of  this  program  is  to  study  methods  for  the  nondestructive  characterization  of 
mechanical  properties  of  two-phase  metal-matrix  composites  such  as  discontinuous  SiC-reinforced 
aluminum  alloys.  The  methods  to  be  studied  are  based  on  the  physical  nonlinear  elastic  behavior 
of  solids  and  utilize  measurements  of  the  elastic  as  well  as  the  acoustoelastic  constants.  The 
objectives  are  to  establish  both  experimentally  and  theoretically  relationships  between  those 
parameters  which  can  be  measured  nondestructively  and  the  percentages  of  second  phase  in 
two-r!'  ase  systems.  Mechanical  properties  of  these  systems  are  derived  from  the  presence  of 
reinforcement  second  phase  in  the  matrix. 

The  research  accomplished  during  contract  No.  DAAL03-88-K-0096  of  The  Army  Research 
Office  program  is  documented  in  details  in  4  theses  (1  -  4)  and  8  papers  (5  -  12).  Abstracts  of 
these  theses  and  reprints  or  preprints  of  these  papers  are  attached  in  appendix  A.  The  major 
accomplishments  obtained  under  this  program  are,  however,  summarized  in  the  following. 

A.  Texture 

Texture  is  the  reason  for  the  direction  dependence  of  materials  properties.  It  is  developed 
during  the  plastic  deformation  and  heat  treatment  during  manufacturing.  It  is  characterized  by  the 
fourth-order  expansion  coefficients  of  the  orientation  distribution  function  (ODF),  C1 14  and  C124 

and  C^4  In  two-phase  materials,  however,  the  situation  is  more  complicated  where  the 

properties,  the  shape,  the  volume  fraction  and  the  orientation  of  each  phase  influence  the 
macroscopic  behavior. 

In  this  study  (1,5,  6),  an  ultrasonic  method  is  developed  for  the  nondestructive 
characterization  of  texture  in  metal-matrix  composites.  The  method  utilizes  the  measurements  of 
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the  six  independent  ultrasonic  velocities  Vij  and  the  formulation  given  by  Bunge.  The  examined 
composites  are  the  siicon  carbide  (Si-C)-particle-reinforced  aluminum  8091,  7094,  and  6061 
metal-matrix  composites.  The  fourth-order  expansion  coefficients  of  the  orientation  distribution 
function  are  found  to  change  linearly  as  a  function  of  the  SiC  content.  Also  in  this  study,  the 
relationship  between  elastic  anistropy  and  texture  in  metal-matrix  composites  has  been  developed 
under  certain  conditions.  Linear  correlations  between  anisotropy  described  by  Young’s  moduli 
and  texture  determined  by  orientation  distribution  function  expansion  coefficients  are  obtained,  and 
found  to  confirm  developed  relationships  for  two-phase  metal-matrix  composites.  These  results 
show  that  ultrasonic  measurements  provide  a  viable  technique  for  the  characterization  of  texture 
and  elastic  anisotropy  in  these  materials. 

B.  Elastic  Properties  and  Elastic  Anisotropy 

In  this  investigation  (2,  7,  8,  11),  longitudinal  and  shear  ultrasonic  velocities  have  been 
measured  in  aluminum  based  metal  matrix  composites  containing  either  SiC-particles  or  alumina 
fibers  and  particles  as  the  reinforcements.  Three  types  (extruded,  pressed,  squeeze  cast)  of 
composites  have  been  investigated.  Also,  based  on  microstructural  studies,  models  have  been 
developed  in  order  to  predict  and  explain  the  elastic  properties  of  these  composites.  From  the 
experimental  measurements  as  well  as  the  comparison  between  measured  and  calculated  elastic 
moduli  the  following  conclusions  are  obtained. 

1.  Extruded  Composites 

a)  For  relatively  low  volume  percentage  of  reinforcement,  up  to  20%,  the  elastic  moduli 
increase  linearly  with  SiC-particle  content. 

b)  The  overall  increase  in  elastic  moduli  due  to  the  reinforcement  can  be  well  predicted  by  a 
model  assuming  dilute  concentration  of  particles. 

c)  The  extrusion  direction  is  elastically  stiffer  than  the  other  two  directions. 

d)  The  higher  value  of  Young's  modulus  along  the  extrusion  direction  is  partly  due  to  the 
aligned  panicle  free  regions  created  during  manufacturing,  and  partly  due  to  texture  in  the 
aluminum  matrix. 
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2.  Pressed  Composites 

a)  For  the  high  volume  fraction  of  SiC -particles,  up  to  40%,  the  increase  in  elastic  moduli 
with  particle  content  is  not  linear  and  deviates  significantly  towards  higher  values. 

b)  The  model  assuming  dilute  concentration  of  particles  fails  completely  to  predict  the 
increase  in  elastic  moduli. 

c)  The  high  values  in  elastic  moduli  can  be  attributed  to  the  increased  interaction  between 
particles  as  the  particle  content  becomes  higher. 

d)  The  elastic  moduli  in  the  compression  direction  is  lower  than  in  the  directions  in  the  plane  of 
the  plate. 

3.  Squeeze  Cast  Composites 

a)  The  elastic  moduli  increase  with  the  alumina  particle  content,  whereas  the  elastic 
anisotropies  between  directions  in  the  plane  of  fibers  and  normal  to  this  plane  decrease  with 
the  particle  content. 

b)  The  overall  increase  in  elastic  moduli  as  well  as  the  elastic  anisotropies  can  be  predicted  by 
a  model  assuming  no  interaction  between  particles  and  fibers. 

c)  The  planar-random  alignment  of  fibers  makes  the  composites  stiffer  in  that  plane. 

C.  Manufacturing 

In  this  work,  the  powder  metallurgy  process  has  been  successfully  implemented  in  the  processing 
of  aluminum  metal  matrix  composites  containing  up  to  30  volume  percent  of  silicon  carbide 
particulate,  and  up  to  20  volume  percent  silicon  carbide  whiskers  (3).  The  composites  processed 
have  been  densified  to  above  97  percent  of  their  theoretical  density,  which  is  comparable  to  those 
obtained  from  commercial  sources.  The  stiffness  properties  of  these  composites  have  been 
evaluated  using  flexure  and  compression  mechanical  tests  as  well  as  ultrasonic  measurements,  and 
are  found  to  match  those  of  commerical  composites. 

Theoretical  calculations  of  the  Young's  modulus  based  on  the  Hashin-Strikman  and  the 
Eshelby  models  are  found  to  be  in  good  agreement  with  the  experimental  results  for  particulate 
content  of  15  percent  or  less.  At  higher  percentages  of  reinforcement,  however,  both  models 
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underpredict  the  stiffnesses.  For  the  whisker  reinforced  composites,  the  Esheiby  model  predicts 
the  stiffnesses  better  than  the  shear  large  model.  The  stiffnesses  predicted  by  the  composite 
cylinders  model  are  also  found  to  agree  pretty  well  with  experimental  results. 

D.  Nonlinear  elastic  Effects 

In  this  study,  the  temperature  dependances  of  nonlinear  effects  effects  are  investigated  by 
measuring  the  acoustoelastic  constants  at  different  temperatures  using  ultrasonic  time  of  flight 
measurements  (4,  9,  12).  These  constants  are  determined  in  the  aluminum  alloys  AlMg  3,  7064, 
and  8091  as  well  as  in  SiC-particle  reinforced  metal  matrix  composites  based  on  the  aluminum 
alloys  7064  and  8091.  In  the  metal  matrix  composite  specimens,  the  acoustic  nonlinearity 
parameter  is  also  measured  using  the  harmonic  generation  technique.  The  results  show  the 
following. 

a)  In  all  materials  investigated,  the  acoustoelastic  constants  show  a  linear  temperature 
dependance.  In  single  phase  alloys,  the  third  order  elastic  constants,  calculated  from 
acoustoelastic  constants  show  an  increase  in  magnitude  for  the  constants  1  and  m  with  temperature. 
However,  the  increase  in  the  constant  1  is  stronger.  The  third  order  elastic  constant  n  is  not 
significantly  influenced  by  temperature.  In  the  case  of  the  two  phase  metal  matrix  composites, 
temperature  changes  affect  the  three  third  order  elastic  constants  decreases  with  increasing 
temperature. 

b)  The  temperature  dependence  of  the  acoustoelastic  effect  confirms  the  linear  relationship 
between  stresses  in  the  material  and  the  temperature  dependance  of  the  longitudinal  velocity  found 
by  Salama  and  Long.  The  parameter  K,  which  characterizes  the  magnitude  of  the 
stress-temperature  dependance  of  ultrasonic  velocities,  is  found  to  be  significantly  influenced  by 
the  base  material  as  well  as  the  alloying  elements  contained. 

c)  The  acoustic  nonlinearity  parameter,  which  is  a  measure  for  the  material's  deviation  from  the 
ideal  Hookean  behavior,  is  calculated  using  the  second  and  third  order  elastic  constants  for  all 
specimens  investigated.  Similar  to  the  behavior  of  the  parameter  K,  the  value  of  the  nonlinearity 
parameter  is  found  to  depend  strongly  on  the  base  material  as  well  as  the  reinforcements. 
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d)  In  two  phase  materials,  it  is  not  sufficient  to  realate  changes  in  the  nonlinearity  parameter  to 
lattice  distortions  only.  Since  the  determination  of  (3  involves  the  ratio  between  third  and  second 
order  elastic  constants,  changes  in  the  elastic  moduli,  due  to  the  presence  of  second  phase,  are  to 
be  taken  into  account.  In  the  case  of  coherent  precipitates,  where  increases  in  the  elastic  moduli  are 
negligible,  the  behavior  of  die  nonlinearity  parameter  is  dominated  by  the  lattice  distortions.  Earlier 
studies  by  Razvi  et  al.  revealed  an  increase  in  the  nonlinearity  parameter  with  the  volume  fraction 
of  second  phase  particles  in  such  materials.  However,  in  the  metal  matrix  composites  investigated 
in  the  present  study,  the  behavior  of  the  nonlinearity  parameter  is  dominated  by  the  increase  of  the 
elastic  moduli. 

e)  Values  of  calculated  nonlinearity  parameters  as  well  as  those  directly  measured  are  found  to 
decrease  linearly  with  increasing  content  of  second  phase  particles.  Also  the  presence  of 
reinforcement  particles  changes  the  temperature  dependance  of  the  nonlinearity  parameter  in  the 
metal-matrix  composites.  In  contrast  to  the  behavior  in  single  phase  alloys,  its  value  decreases 
with  increasing  temperature.  Changes  in  the  distortion  of  the  matrix  are  believed  to  be  responsible 
for  this  opposite  behavior.  A  temperature  increase  leads  to  relaxation  of  the  thermal  stresses  which 
are  always  present  in  metal  matrix  composites,  due  to  the  mismatch  of  the  coefficients  of  thermal 
expansion  of  the  matrix  and  the  reinforcement  materials.  The  relaxation  of  these  stresses  reduces 
the  distortion  of  the  matrix  and,  thus,  the  value  of  the  nonlinearity  parameter  decreases. 

E.  Scanning  Acoustic  Microscope  (SAM)  Studies 

In  addition  to  the  above  macroscopic  studies,  preliminary  acoustic  microscopy  studies  (13) 
on  samples  of  Al-8091  and  Al-7064  show  that  the  SAM  has  no  difficulty  at  1.0  GHz  in  resolving 
the  particles  in  both  composites.  The  resolution  in  the  GHz  range  is  on  the  order  of  a  micron,  so 
that  detection  of  the  particles  is  not  a  difficulty.  In  the  characterization  of  the  interface  the 
assumption  is  made  that  changes  in  the  stress  field  will  produce  measurable  changes  in  the  region 
of  the  interface  and  there  will  be  a  corresponding  change  in  the  stress  between  the  two  phases  due 
to  the  difference  in  the  thermally  induced  strain  between  them.  This  will  cause  changes  in  the 
sound  propagation  properties  in  the  stressed  regions.  Although  these  changes  are  small  they 
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depend  on  the  properties  of  the  interface  and  can  be  measured. 

The  affected  region  can  range  from  nanometers  to  micrometers  in  thickness  even  in  the  GHz 
range  and  the  properties  will  vary  throughout  only  one  acoustic  wavelength.  Point  to  point 
changes  along  the  boundary,  however,  cause  significant  changes  in  the  shape  of  the  fringe  pattern 
due  to  the  changing  of  the  boundary  conditions  at  the  reflecting  interface.  In  addition,  many  of  the 
panicles  exhibit  surface  wave  scattering  fringes  over  several  wavelengths,  indicating  higher  phase 
sensitivity. 
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Abstract 


Many  important  material  properties  are  of  nonlinear  nature.  Among  them  are  the  stress- 
temperature  dependance  of  ultrasonic  velocities,  the  distortion  of  elastic  waves  by  the 
generation  of  higher  harmonics,  and  the  thermal  expansion.  These  effects  are  characterized 
quantitatively  by  measurements  of  the  acoustoelastic  constants,  the  acoustic  nonlinearity 
parameter,  and  the  coefficient  of  thermal  expansion. 

In  this  study,  the  temperature  dependances  of  nonlinear  ultrasonic  effects  are 
investigated  in  the  cemented  carbide  WC-Co,  the  ferritic  steel  24  CrMoV  5  5,  the  austenitic 
steel  X6  CrNi  18  11,  the  aluminum  alloys  AlMg  3,  A1  7064  and  A1  8091  as  well  as  in 
metal-matrix  composites  consisting  of  A1  7064  and  AI  8091  matrices  and  SiC-particles. 
The  results  show  that  the  magnitude  of  acoustic  nonlinearity  depends  on  the  base  material, 
alloying  elements,  and  microstructure.  It  increases  with  the  temperature  in  single  phase 
materials  whereas  it  decreases  in  metal-matrix  composites.  The  various  contributions  to  the 
acoustic  nonlinearity  are  analyzed  qualitatively  and  possible  explanations  for  its  behavior  in 
metal-matrix  composites  are  discussed. 
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Abstract 


Textures  or  preferred  orientations  of  single  crystals  in  polycrystalline  materials  are 
developed  by  heat  treatment  and  deformation  processes  like  rolling,  drawing  and 
extrusion.  Thus  magnetic,  elastic  and  plastic  properties  become  directionally  dependent 
This  can  yield  desirable  behavior,  but  it  can  have  also  negative  effects.  In  both  cases  a 
nondestructive  determination  of  the  texture  and  its  consequences  for  materials  behavior  is 
of  great  interest 

In  this  study  texture  is  analyzed  by  determining  the  fourth-order  expansion 
coefficients  of  the  orientation  distribution  function  from  the  measured  ultrasonic 
velocities  of  different  wave  modes.  The  results  obtained  for  rolled  ferritic  steel  sheets, 
rolled  ferritic  steel  plates  and  for  extruded  metal-matrix  composites  of  the  aluminium 
alloys  Al-8091,  Al-7064  and  Al-6061  with  silicon  carbide  particle  reinforcements  show 
that  ultrasonics  provide  an  efficient  nondestructive  method  of  texture  analysis  in  the  bulk 
as  well  as  on  the  surface  of  these  materials.  The  expansion  coefficients  determined  for 
the  ferritic  steels  agree  qualitatively  with  those  determined  by  x-ray  diffraction  which  is  a 
common  nondestructive  technique  for  texture  analysis.  For  the  metal-matrix  composites 
the  expansion  coefficients  have  been  determined  under  the  assumption  that  only  the 
Al-matrix  is  textured.  The  results  also  show  that  the  fourth-order  expansion  coefficients 
allow  the  evaluation  of  the  elastic  and  the  plastic  behavior  of  the  examined  specimens. 
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ABSTRACT 


The  effects  of  varying  the  reinforcement  type  and  content  on  the  elastic  behavior  of 
aluminum  /  silicon  carbide  composites  have  been  investigated  in  order  to  develop 
relationships  to  be  used  in  controlling  the  integrity  of  the  composites.  A  canless  powder 
metallurgy  (PM)  process  was  successfully  developed  for  the  in-house  manufacture  of  the 
MMC  specimens  .  Composite  specimens  with  up  to  30  v/o  particulates  and  up  to  20  v/o 
whiskers  were  fabricated  using  A1  6061  as  the  base  metal.  The  elastic  moduli  of  these 
MMCs  were  characterized  using  ultrasonic  velocity  measurements  as  well  as  mechanical 
testing  involving  flexural  and  compressive  tests.  The  measured  moduli  were  also 
compared  with  model  predictions. 

The  results  indicate  that  the  addition  of  the  particulate  or  the  whisker  reinforcements 
increase  the  elastic  moduli  as  well  as  the  elastic  anisotropy  of  the  MMCs.  The  composites 
were  elastically  suffer  in  the  plane  perpendicular  to  the  hot  pressing  direction  than  in  the 
pressing  direction  and  also  exhibited  transverse  isotropy  about  the  pressing  axis.  While  the 
anisotropy  of  the  particulate  reinforced  MMCs  is  attributed  to  the  uneven  distribution  of  the 
reinforcement  along  the  different  directions,  the  anisotropy  of  the  whisker  reinforced  MMC 
is  traced  to  the  planar  orientation  of  the  whiskers  in  the  plane  perpendicular  to  the  pressing 
direction.  Also,  it  is  found  that  the  substitution  of  particulate  reinforcement  by  whiskers 
does  not  improve  the  elastic  stiffness  of  the  MMCs  for  this  method  of  composite 
manufacture. 
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ABSTRACT 

Conventionally,  metal  matrix  composites  (MMC)  are  reinforced  with 
either  particles  or  fibers,  but  lately,  a  new  class  of  composites  has  emerged, 
where  a  mixture  of  particles  and  fibers  is  used  as  a  reinforcement.  The  particles 
are  present  to  improve  the  overall  mechanical  and  thermal  properties,  whereas 
the  fibers  introduce  the  directionality  often  desired  in  some  applications. 

The  elastic  behavior  of  three  composites  with  different  matrices  and 
volume  fraction  of  particles  (9,13  and  17%),  but  the  same  fiber  content  (6%), 
has  been  characterized  using  ultrasonic  velocity  measurements.  The  results 
show  that  the  elastic  moduli  increase  with  particle  content  and  that  the 
composites  are  elastically  stiffer  in  the  directions  of  the  plane  of  the  fibers  than 
in  the  direction  normal  to  that  plane. 

A  model  is  developed  to  explain  the  observed  elastic  moduli  of  these 
composites.  This  model  uses  the  results  of  theories  presented  by  Ledbetter  and 
Datta,  for  spherical  inclusions,  and  by  Hashin  and  Rosen,  for  aligned  fibers.  It 
also  includes  an  average  procedure  suggested  by  Christensen  and  Waals.  The 
agreement  between  measured  and  calculated  elastic  moduli  :s  found  to  be  very 
good  and  the  elastic  anisotropies  observed  in  these  composites  could  also  be 
predicted. 

Elastic  moduli  are  also  determined  using  ultrasonic  velocity 
measurements  for  two  series  of  extruded  MMC's  and  one  series  of  pressed 
MMC's.  These  composites  are  reinforced  with  SiC-particles  only.  The  elastic 
moduli  as  well  as  the  elastic  anisotropies  in  these  composites  could  be 
explained  using  combinations  and/or  special  cases  of  the  theories  mentioned 
above. 
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Texture  of  Metal-Matrix  Composites 
by  Ultrasonic  Velocity  Measurements 
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Abstract.  An  ultrasonic  method  is  developed  for  the  nondestructive  characterization  of 
texture  in  metal-matrix  composites.  In  this  approach,  it  is  assumed  that  the  presence  of 
reinforcement  particles  changes  the  elastic  properties  of  the  composite  but  only  the 
texture  of  the  matrix.  The  method  utilizes  the  measurements  of  the  six  independent 
ultrasonic  velocities  l'y  and  the  formulation  given  by  Bunge.  The  examined  composites 
are  the  silicon  carbide  (SiC)-particle-reinforced  aluminum  8091.  7064.  and  6061  metal- 
matrix  composites.  The  fourth-order  expansion  coefficients  of  the  orientation  distribu¬ 
tion  function  are  determined  as  a  function  of  the  SiC  content  in  these  composites.  The 
results  show  that  the  expansion  coefficient.-,  change  with  the  presence  of  SiC  where  the 
coefficients  C.J1  and  CV'  increase  as  the  volume  fraction  of  SiC  is  increased  and  the 
coefficient  C.|:  is  zero  in  all  composites  examined.  The  analysis  of  these  results  indicates 
that  ultrasonics  can  provide  a  promising  technique  for  the  texture  characterization  of 
metal-matrix  composites.  v 


Introduction 

Texture  is  the  orientation  distribution  of  the  single  crystals  in  the  polycrystal¬ 
line  aggregate.  A  textured  polycrystal  is  elastically  anisotropic  because  the 
elastic  properties  of  a  single  crystal  are  directionally  dependent.  Because  of  the 
texture,  the  single  crystal  anisotropies  do  not  vanish  when  averaged,  thus  the 
polycrystal  looses  its  quasi-isotropy.  Most  structural  materials  are  polycrystal- 
linc  aggregates  and  their  exposure  to  plastic  deformation  and  heat  treatment 
during  manufacturing  leads  to  the  alignment  of  single  crystals  relative  to  the 
forming  geometry.  This  development  of  preferred  orientation  is  the  main  rea¬ 
son  for  the  anisotropic  behavior  in  these  materials.  The  texture  is  characterized 
by  the  fourth-order  expansion  coefficients  of  the  orientation  distribution  func¬ 
tion  (ODF).  These  three  coefficients  allow  a  satisfactory  evaluation  of  the 
elastic  as  well  as  the  plastic  behavior  of  polvcrystalline  aggregates,  as  has  been 
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shown  by  several  studies  (e.g.,  1 1—5]).  This  situation,  however,  is  more  compli¬ 
cated  in  the  case  of  two-phase  or  multiphase  materials,  where  the  properties, 
the  shape,  the  volume  fraction,  and  the  orientation  of  each  phase  influence  the 
macroscopic  behavior. 

In  metal-matrix  composites  a  ductile  metal  and  a  high  strength  reinforce¬ 
ment  are  combined  to  provide  a  composite  of  high  strength  and  toughness. 
Since  many  of  these  properties  arc  characteristic  of  the  bulk,  ultrasonics  have 
been  shown  to  provide  promising  nondestructive  methods  for  the  characteriza¬ 
tion  of  metal-matrix  composites  16). 

In  this  study,  an  ultrasonic  method  is  developed  for  the  nondestructive 
characterization  of  texture  in  metal-matrix  composites.  Measurements  of  ultra¬ 
sonic  velocities  are  used  to  determine  the  texture  in  the  silicon  carbide  ( SiC )- 
reinforced  8091.  7064.  and  *n61  aluminum  composites. 


Quantitative  Texture  Analysis  Using  Ultrasonics 

Orientation  Distribution  Function 

Texture  is  mathematically  described  by  the  ODF.  This  function  determines  the 
probability  of  finding  a  single  crystal  in  the  polycrystalline  sample  with  a  certain 
orientation  with  respect  to  the  sample  orientations,  given  by  the  axes  of  the 
sample  fixed  coordinate'  system.  According  to  Bunge  [7 1 .  the  ODF  can  be 
written  as  a  series  expansion  into  symmetrical  generalized  spherical  harmonics 
as 

•  ,•./()>  Mil 

v  v  v  cf'trtg)  in 

I  O  n  I  i'-  I 


where  i>  is  the  orientation,  the  T\'  are  symmetrical  generalized  spherical  har¬ 
monics  and  the  Cf  are  the  expansion  coefficients.  The  upper  limits  Mi  1)  and 
Ml)  depend  on  1  as  well  as  the  crystal  and  sample  symmetry,  respectively.  The 
orthonormal  function  system  7Yr  is  invariant  towards  all  rotations  of  the  sample 
symmetry  (indicated  by  the  right  dot)  and  the  crystal  symmetry  (indicated  by 
the  left  two  dots).  In  first  approximation  only  the  three  fourth-order  expansion 
coefficients  C4\  C4:  and  C{3  need  to  be  considered  for  texture  evaluation  in 
cubic  materials  with  orthorhombic  sample  symmetry  [1,8].  These  three  coeffi¬ 
cients  can  be  used  to  characterize  both  the  elastic  and  the  plastic  behavior  of 
these  materials.  The  fourth-order  coefficients  in  the  Roc-notation  [9]  arc  called 
Wjoo,  W 420.  and  W440  and  are  related  to  Bunge's  coefficients  by  the  following 
expression: 
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Ultrasonic  Velocity  Relationships 

Inserting  the  elastic  constants  of  the  textured  polycrystal,  given  by  Bunge  1 1| 
for  cubic  crystal  structure,  into  Christoffel’s  equation  for  the  orthorhombic 
sample,  the  following  relationships  between  the  velocities  of  ultrasonic  waves 
V ij  and  the  expansion  coefficients  Cl1'  result  as 


pV It  -  cn  -  c 


I  -~i  -  |v-5  c;:  -  i  x  55  cj’l 
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where  p  is  the  density  and  cn.  t  i;.  and  e44  are  the  elastic  constants  of  the  cubic 
single  crystal  and  c  =  t'n  -  ct;  -  2c\u.  In  these  relationships  the  velocity  is 
characterized  by  two  subscripts.  The  first  indicates  the  propagation  direction, 
while  the  second  indicates  the  polarization  directiot]  of  the  wave.  I,  2.  and  3 
designate  the  axes  in  a  right-handed  coordinate  system  as  shown  in  Fig.  I. 

The  Lame  constants  A  and  p  for  the  unlexturcd  polycrystal  (Cl1  =  C]:  = 
C.13  =  0)  can  be  expressed  as 


p  =  <44  +  7t  <••  A  +  2p  =  <•,,  -  ^  <•.  A  =  C i; 


c . 


For  the  isotropic  case  (c  =  0),  the  shear  wave  velocity  Vv  as  well  as  the 
longitudinal  wave  velocity  V\  will  be  the  same  in  all  directions,  and 


v 


FIr.  I.  Velocity  designations  for 
free  ultrasonic  waves. 


Spies  and  Salama 


102 


Table  1.  Chemical  composition  of  aluminium  alloys  and  volume  percentage  of 
SiC  reinforcement  of  the  MMC-specimens 


Alloying  Llemcnt 

S 

Alloy 

Si 

Fe 

Cu 

Mg  Xr  Li 

Xn 

Ci 

Co  Al 

8091 

0.02 

0.01 

1.90 

0.80  0.11  2.70 

_ 

—  rein 

7064 

0.0.' 

0.10 

2.00 

2.30  0.20  — 

7.10 

0.12 

0.22  rern 

6061 

0.71 

0.29 

0.21 

0.86  —  _ 

0.10 

0.06 

—  rem 

Specimens 

8091  + 

or;  sic 

8091  -  tor;  sic 

8091 

-  i5';  sic 

7064  ■*• 

or;  sic 

7064  -*•  15<;  SiC 

7064 

-  2o< ;  sic 

6061  -s  i5r;  sic 

6061 

*  2o<;  sic 

/A'i  -  <11  -  A  -  2 /I. 
f>\';  -  <44  =  m- 

A  further  useful  expression  which  can  he  obtained  from  Eqs.  (3)  is 

/>(  V  ii  —  \  ;;  “  \  >0  =  <  |  |  —  2<44.  (4 1 


Experiment 

The  metal-matrix  composites  examined  in  this  investigation  arc  the  aluminium 
alloys  Al-8091.  Al-7064.  and  Al-6061  containing  up  to  209r  volume  fraction  of 
SiC  particles.  The  chemical  compositions  of  these  alloys  and  the  volume  per¬ 
centages  of  the  SiC  reinforcement  are  shown  in  Table  I.  The  Al-8091  and  Al- 
7064  composites  were  received  as  extruded  rods  of  25  mm  in  diameter,  while 
the  Al-6061  composites  were  received  as  extruded  bars  of  16  mm  in  thickness 
and  50  mm  in  width.  Prismatic  specimens  15  to  18  mm  wide  and  30  mm  long 
were  machined  such  that  the  length  direction  was  parallel  to  the  extrusion 
direction.  Opposite  faces  were  machined  flat  and  parallel  to  within  ±0.025  mm. 
All  these  specimens  were  examined  in  the  as-received  condition. 

In  order  to  determine  the  ODF-expansion  coefficients,  ultrasonic  velocity 
measurements  are  performed  using  the  pulse-ccho-ovcriap  method,  which  is 
described  in  detail  by  Papadakis  |I0).  The  system  used  in  this  investigation 
consists  of  an  ultrasound  apparatus  which  genc.ates  pulses  of  approximately  I 
ps  duration  and  of  variable  repetition  rate.  These  pulses  arc  impressed  on  a 
commercial  transducer  of  a  fundamental  frequency  of  5  MHz  which  is  acousti¬ 
cally  bonded  to  the  specimen.  The  reflected  echoes,  rf.  arc  received  by  the 
same  transducer,  amplified,  and  displayed  on  a  screen  of  an  oscilloscope.  Two 
of  the  displayed  echoes  are  then  chosen  and  exactly  overlapped  by  critically 
adjusting  the  frequency  of  a  C.W.  oscillator  and  the  division  factor  on  a  decade 
divider.  The  frequency.  /'.  accurately  determined  by  tin  electronic  counter  is 
employed  to  compute  the  ultrasonic  velocity  using  the  relationship  V  -  sf, 
where  >  is  the  path  length  of  the  ultrasonic  wave.  The  system  is  capable  of 
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measuring  changes  in  the  ultrasonic  velocity  to  an  accuracy  of  better  than  I 
part  in  I05,  while  the  inaccuracy  of  measuring  both  the  longitudinal  and  the 
shear  velocities  is  estimated  to  be  0.3%. 


Results 

Reduction  of  Data 

The  conversion  of  the  measured  velocity  data  to  the  texture  expansion  coeffi¬ 
cients  is  very  critical.  Equations  (3).  which  relate  the  velocities  Vtfof  shear  and 
longitudinal  waves  with  different  polarization  and  propagation  directions  to  the 
fourth-order  expansion  coefficients,  provide  a  variety  of  relationships  to  deter¬ 
mine  these  coefficients.  These  relationships  are  shown  in  Table  2.  In  order  to 
select  the  relationships  from  which  the  coefficients  are  to  be  determined,  the 
propagation  of  errors  involved  in  the  procedure  of  reducing  the  measured  data 
to  the  physically  significant  quantities  becomes  very  important. 

From  Table  2  one  can  see  that  there  are  three  relationships  which  deter¬ 
mine  Ci1  and  four  equations  to  determine  each  of  C\2  and  C|\  The  usual  proce¬ 
dure  to  determine  these  coefficients  from  the  relationships  shown  in  Table  2  is 
to  use  the  least-squares  analysis.  However,  in  order  to  keep  the  number  of 
velocities  to  be  measured  low  and  to  reduce  the  effect  of  propagation  of  errors, 
the  following  procedure  is  employed: 

1.  We  select  relationships  which  do  not  require  the  absolute  values  of  veloci¬ 
ties  and  only  require  the  differences  in  these  velocities.  Thus  errors  in  the 


Tabic  2.  Different  relationships  for  the  determination  of  the 'expansion  coefficients  Cj‘ 
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velocity  calculations,  which  are  due  to  measurements  of  the  path  length,  are 
cancelled  out, 

2.  We  use  relationships  which  involve  the  lowest  number  of  velocities. 

3.  We  select  relationships  which  require  only  ratios  of  elastic  constants.  If  this 
is  not  possible  (as  in  the  cases  of  C]1  and  C];)  the  equations  which  require  the 
lowest  number  of  elastic  constants  are  used. 

4.  We  use  the  set  of  relationships  which  requires  the  lowest  number  of  mea¬ 
surements  necessary  to  determine  the  coefficient. 

The  relationships  w  hich  best  satisfy  the  above  procedure  are: 


- ,  210  ^ 

C  =  —  ^  [2m  -  p(V-„  -  V  ;,,>].  (51 

C\:  =  ApiVh  ~  VU.  (6) 

4Y.V  '  7 

Tin  ^ 

C\'  -  -^==-  v;7  [6m  -  A  -  ptl'j,  ~  8 Vi;)].  (7) 

Replacing  the  density  p  hv  the  expression  in  Eq.  (4)  yields 
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It  must  be  noted  that  Eq.  (A)  of  Tabic  2  also  satisfies  options  1-4  listed  above. 
Replacing  the  density  p  by  the  expression  in  Eq.  (4),  the  resulting  expression 
for  C’V  leads  to  the  same  values  obtained  from  Eq.  (8). 

Equations  (8-10)  are  valid  only  for  materials  with  cubic  crystal  structure, 
and  in  order  to  apply  these  equations  to  the  composites  examined  in  this  work, 
the  following  assumptions  are  made:  I)  The  presence  of  the  SiC-parlicles 
changes  the  elastic  properties  of  the  composite  but  only  the  texture  of  the 
matrix:  2)  The  SiC-particles  are  randomly  distributed  in  the  specimen  without 
any  preferred  orientation:  .3)  The  ultrasonic  w'avclcngth  (frequency  used  is  5 
MHz)  is  much  larger  than  the  average  particle  dimension  (2-4  pm)  so  that  no 
dispersion  effects  occur. 


Determination  of  the  Expansion  Coefficients  C\' 

The  geometry  of  the  specimens  used  in  this  investigation  allows  the  measure¬ 
ments  of  all  nine  ultrasonic  wave  velocities  V„  (3  different  waves  in  3  orthogo¬ 
nal  directions).  In  these  measurements  the  3-dircction  is  taken  to  be  parallel  to 
the  extrusion  direction  of  the  samples  (see  Fig.  I ).  The  velocities  measured  are 
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t'a 

.’477 

7702 

7X05 

7066 

7440 

7567 

74 IX 

7701 

t  ... 

.’>467 

7728 

7X74 

70X2 

7479 

764 1 

747X 

7771 

f;. 

'4X0 

7770 

7X42 

7057 

7469 

7676 

7442 

7777 

t  r. 

7472 

3733 

7X4X 

7051 

7470 

7617 

7446 

7740 

f.i 

74X2 

7776 

7X44 

7062 

7469 

7640 

7475 

7774 

1  Velocities  are  shown  in  units  of  m  s. 


listed  in  Table  3  and  arc  found  to  be  reproducible  to  within  0.3 r/c.  Because  of 
the  rotational  axis,  the  corresponding  velocities  V'n  and  V';s,  VN:  and  V 2|.  IN] 
and  V|t.  and  V:y  and  VN2  are  averaged  according  to  the  rotational  symmetry 
around  the  extrusion  axis. 

The  ultrasonic  velocities  in  the  textured  Al-matrix  (which  has  a  face-cen¬ 
tered  cubic  crystal  structure)  are  then  computed  using  the  equal  stress  condi¬ 
tion  and  the  velocities  measured.  This  condition  assumes  that  the  Al-crystal- 
lites  and  the  SiC-particIcs  are  subjected  to  the  same  stress.  With  .v  denoting  the 
volume  percentage  of  SiC.  the  equal  stress  condition  can  be  written  as 


3*  comp 


.V/.u  A/sic 


,v,V/A|  +  ( I  -  .vl/V/sjf’ 


(II) 


where  A/comp,  M.\\,  and  Ms\c  are  the  elastic  moduli  of  the  composite,  the  matrix 
and  the  reinforcement,  respectively.  Using  Eq.  (II)  and  the  well-known  rela¬ 
tionships  between  elastic  moduli  and  ultrasonic  velocities,  the  velocities  Y„  in 
the  Al-matrix  can  be  expressed  as 

v.  A I  =  ±  (I  ~  ■V)Pcomp(V/"mi,,):pslcV^,(-  '  (  |2) 

tPAI  (> Sick'sif  -  .'Pcoinpt  V 

In  the  calculations  of  VAI  the  density  of  the  composite  is  determined  according 
to  the  law  of  mixture 


Pcomp  =  (l  -  V)p,\|  +  .vpSic.  (13) 

Only  small  deviations  of  the  densities  calculated  using  Eq.  (13)  and  of  the 
measured  densities  are  found  for  the  specimens  examined  in  this  work  III).  The 
densities  used  in  the  calculations  are  2524  kg/m1  for  Al-8091,  2864  kg/nv1  for  Al- 
7064.  2710  kg/nv1  for  Al-6061  (12).  and  3200  kg/nv1  for  SiC  1 13).  The  shear  and 
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Table  4.  Expansion  coefficients  Cj‘  of  the  AI-8091.  A I -7064  and  AI-6061  MMC' 
specimens' 


8091 

+ir; 

SiC 

S09I 
+  it)';; 
SiC 

8091 

+ 15'; 
SiC 

7064 

+  0'; 
SiC 

7064 
+ 15';; 
SiC 

7064 
+  20'.  ; 
SiC 

6061 

+ 15'; 
SiC 

6061 

+  20' ; 
Sic 

2.205 

1.573 

3.532 

2.358 

3.155 

5.352 

3.295 

7.404 

cf 

0 

0 

0 

0 

0 

0 

0 

0 

cr 

2.194 

-0.431 

1 .034 

1 .962 

1.267 

2.175 

1.913 

4.f»3h 

1  3-direction  |!  extrusion  direction. 


the  longitudinal  wave  velocities  V’s  and  V)  tor  SiC  are  obtained  from  the  rela¬ 
tionships  -  /i//).  Vj  --  (X  -t  2/ti)/p  and  v  =  X/|2(X  +  /x)|.  where  v  is  the 
Poisson's  ratio  =  0.19  and  /jl  =  165.5  GPa  (13). 

The  single  crystal  values  used  in  Eqs.  (8)-(I0)  are: 

Ci  -  111. I  GPa.  <  ,:  -  52.5  GPa.  <•«  =  32.5  GPa.  r  =  -6.4  GPa  for  AI-8091  [14]. 

d  =  107.3  GPa.  c;  =  60.3  GPa.  c44  =  28.0  GPa.  <  =  -9.0  GPa  for  Al-7064  [15]. 

Ci  -  108.0  GPa.  <•,:  -  OlO  GPa.  <  44  -  28.3  GPa.  <•  -  -  10.6  GPa  for  AI-6061  [161. 

The  expansion  coefficients  CS‘  are  then  determined  using  Eqs.  (8-10).  and  their 

results  are  shown  in  Table  4.  Because  of  the  rotational  symmetry  of  the  sam¬ 
ples.  the  coefficient  C\:  is  found  to  be  zero  for  all  specimens. 


Discussion 

The  expansion  coefficients  Cl1  and  Ci‘  are  plotted  versus  the  volume  percent¬ 
age  of  SiC  in  Fig.  2.  From  the  figure  one  can  sec  that  there  is  first  a  slight 
decrease  in  these  coefficients  followed  by  a  linear  increase  up  to  20  volume 
percent  of  SiC.  These  plots  also  show  that  the  presence  of  SiC  leads  to  consid¬ 
erable  changes  in  these  expansion  coefficients  and  thus  in  the  texture  of  the  Al- 
matrix. 

The  determination  of  Ci1  is  known  to  be  critical  because  it  requires  absolute 
velocity  measurements.  Therefore  it  must  be  established  whether  the  deter¬ 
mined  values  for  Ci1  are  reasonable.  This  can  be  done  by  considering  the  two 
extreme  cases  of  a  <1 1  l)-fiber  texture  and  a  (100)-fibcr  texture,  which  are  usu¬ 
ally  observed  in  extruded  aluminium  [17],  Because  of  their  rotational  symme¬ 
try.  fiber  textures  are  sufficiently  described  by  one  expansion  coefficient, 
namely.  Cl1,  whereas  Cl:  and  CV  are  zero  1 18).  In  the  case  of  an  ideal  <11!  >-fiber 
texture,  where  the  crystallographic  directions  with  the  highest  ultrasonic  veloc¬ 
ity  lie  along  the  extrusion  direction,  the  value  of  Cl1  is  found  to  be  equal  to 
-4.575.  On  the  other  hand,  for  ani  deal  <10())-fibcr  texture,  where  the  crystallo¬ 
graphic  direction  with  the  lowest  ultrasonic  velocity  lies  along  the  extrusion 
direction,  the  value  of  Ci1  is  +6.887.  These  values,  which  arc  obtained  from  the 
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Fig.  2.  Expansion  coelli- 
cients  (  anil  (T  (hulk 
values!  determined  for  the 
MMC  specimens  plotted 
versus  the  volume  per¬ 
centage  of  SiC. 

s  io  15  :n  25 

Volume  percentage  of  SiC 


results  of  pure  aluminium,  agree  very  well  with  Ci1  =  -4.570  and  C]1  ~  +6.770 
given  by  Bunge  [I8J.  These  two  “extreme"  values  constitute  the  upper  and  the 
lower  bounds  for  Ci1  of  any  texture  in  extruded  aluminium. 

The  values  obtained  in  this  work  (Table  4)  lie  within  these  bounds,  except 
for  the  AI-6061  +  20%  SiC  composite.  This  agrees  with  the  results  reported  in 
reference  [II]  which  show  (hat  the  elastic  behavior  of  all  Al-SiC  composites 
examined  in  this  work  can  be  approximately  described  by  the  equal  stress 
condition  represented  by  Eq.  (II).  with  the  exception  of  the  AI-6061  +  20%  SiC 
composite  which  shows  a  large  deviation  from  that  condition.  The  velocities  \ 
in  the  Ai-malrix  were  also  determined  using  the  equal  strain  condition 

44 comp  —  1 1  A  )/W  \|  ***  A.V/sjO 

which  assumes  that  the  Al-crystallitcs  and  the  SiC-particles  undergo  the  same 
strain.  The  coefficients  resulting  from  these  calculations  are  found  to  lie  far 
outside  the  bounds  given  above,  which  again  agrees  with  the  results  in  refer¬ 
ence  1 1 II. 

The  departure  of  the  elastic  behavior  of  the  AI-6061  +  20%  SiC  composite 
from  the  equal  stress  condition  may  be  due  to  some  degree  of  hydrostatic  stress 
generated  within  the  MMC  specimen.  In  this  case,  the  matrix  material  in  the 
composite  tends  to  deform  at  a  lower  stress  than  the  reinforcement  because  the 
matrix  material  is  softer  than  the  reinforcing  material.  If  the  matrix  material  is 
rigidly  coupled  with  the  reinforcement,  it  is  restricted  from  deforming  in  the 
way  it  would  if  it  were  alone  and  causes  hydrostatic  stresses  to  be  generated 
within  the  matrix  material.  The  exact  magnitude  of  this  effect  is  unknown  and 
rather  complex,  but  it  is  a  function  of  the  mean  free  path  in  the  matrix  and  the 
ratio  of  the  clastic  constants  of  the  constituents  |I9],  In  general  the  mean  free 
path  in  the  matrix  varies  from  region  to  region,  and  therefore  no  simple  ratio 
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can  be  predicted.  Thus  the  extent  of  deviation  of  the  elastic  moduli  form  the 
equal  stress  condition  can  vary  from  specimen  to  specimen. 

Since  the  mean  free  path  in  the  matrix  decreases  with  increasing  volume 
percentage  of  the  reinforcing  phase,  the  limitation  of  the  approach  presented  in 
this  work  to  composites  with  "small"  fractions  of  reinforcements  is  apparent. 
Nevertheless,  in  these  cases  reasonable  results  arc  obtained  up  to  a  volume 
fraction  of  209?.  From  Table  4  and  Fig.  2  one  can  see  that  the  coefficients 
determined  lie  within  the  bounds  given  above  and  that  the  coefficients  of  these 
specimens,  which  contain  the  same  volume  percentage  of  SiC  lie  close  to¬ 
gether.  indicating  that  these  specimens  contain  similar  textures.  Figure  2  also 
indicates  that  the  textures,  which  are  only  due  to  the  extrusion  process  (09? 
SiC).  are  modified  by  the  addition  of  the  SiC-rcinforcemcnts  in  a  similar  way.  If 
it  is  assumed  that  the  examined  specimens  contain  a  mixture  of  a  (111)-  and  a 
(lOO)-fibcr  texture,  which  is  usually  observed  in  extruded  aluminium  117).  the 
increase  in  the  expansion  coefficient  CV  with  increasing  volume  percentage  of 
SiC  indicates  that  the  <  I00)-directions  of  the  Al-cryslalliles  align  more  and  more 
towards  the  extrusion  direction  as  the  SiC'-conccntration  is  increased. 

From  above,  it  can  be  seen  that  ultrasonic  velocity  measurements  can  be 
used  for  the  nondestructive  characterization  of  texture  in  metal-matrix  com¬ 
posites.  The  results  also  indicate  that  the  assumptions  made  concerning  the 
effects  of  reinforcement  particles  on  the  texture  in  these  materials  arc  reason¬ 
able. 
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Figure  9  Experimental  (•)  and  simulated  (  )  operational 

impedance  characteristics  for  a  three-layered  lateral  structure. 
Frequency  range  400  kHz  4  MHz 

Concluding  remarks 

A  model  has  been  presented  for  predicting  the  response 
from  multilayered  piezoelectric  structures  which  demon¬ 


Five  port  lattice  model:  G.  Hayward  and  D.  Gillies 

strate  two.  loosely  coupled,  compressionai  modes  of 
vibration.  Good  agreement  between  theoretical  and 
experimental  results  was  achieved  over  a  range  of  structural 
configuration  ratios.  Provided  that  the  underlying  physical 
assumptions  are  appreciated,  the  model  is  useful  for 
providing  a  first  order  approximation  to  electrical 
mechanical  and  piezoelectric  interaction  in  ultrasonic 
array  structures.  For  example,  the  model  will  provide  a 
good  approximation  for  thickness  mode  behaviour  within 
individual  elements  in  a  diced  array.  The  influence  of 
bond  lines,  mechanical  matching  layers  and  electrical 
loading  on  transducer  response  may  be  investigated  in 
both  time  and  frequency  domains.  In  addition,  lateral 
mode  propagation  may  be  studied  in  conjunction  with 
factors  such  as  element  periodicity,  inter-element  loading 
and  electrical  cross  coupling.  This  approach  is  also 
suitable  for  modelling  of  pure  thickness  mode  structures 
which  are  subject  to  independent  electrical  loading  by 
setting  the  appropriate  mechanical  and  piezoelectric  cross 
coupling  constants  to  zero.  An  extension  of  the  method 
for  the  simulation  of  stacked  piezoelectric  structures 
which  possess  common  or  parallel  electrical  loading  has 
also  been  developed  and  will  be  reported  at  a  later  date. 
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Relationship  between  elastic  anisotropy  and 
texture  in  metal-matrix  composites 
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The  relationship  between  elastic  anisotropy  and  texture  in  two  phase  metal  matrix  composites 
has  been  developed  under  certain  conditions.  Using  measurements  of  the  six  independent 
ultrasonic  velocities  V,t  in  samples  of  the  aluminium  alloys  8091  and  7064  containing  up  to 
20  %  SiC  particles  and  the  formulation  given  by  Bunge,  the  fourth-order  expansion 
coefficients  of  the  orientation  distribution  function  are  determined.  The  Young's  moduli  in 
different  directions  are  also  obtained  from  ultrasonic  velocity  measurements.  Linear  correlations 
between  anisotropy  described  by  Young's  moduli  and  texture  determined  by  orientation 
distribution  function  expansion  coefficients  are  obtained,  and  confirm  developed  relationships 
for  two  phase  metal  matrix  composites.  This  result  shows  that  ultrasonic  measurements 
provide  a  technique  for  the  characterization  of  texture  and  elastic  anisotropy  in  these 
materials. 

Keywords:  metal-matrix  composites;  elastic  anisotropy;  texture 


Introduction 

Most  structural  materials  arc  polycry  stalline  aggregates 
and  their  exposure  to  plastic  deformation  and  heat 
treatment  during  manufeluring  leads  to  the  alignment  of 
the  single  crystals  relative  to  the  forming  geometry.  This 
development  of  preferred  orientation  or  texture  is  the 
mail,  reason  for  the  anisotropic  behaviour  of  these 
materials.  I  he  texture  is  mathematically  described  by  the 
orientation  distribution  function  (ODE  I.  which  determines 
the  probability  of  finding  a  single  crystal  in  the  aggregate 
with  a  certain  orientation  with  respect  to  the  sample 
geometry  Ihc  elastic  behaviour  of  polycrystalline 
materials  is  determined  by  the  three  lourth-order 
expansion  coellicients  of  the  ( )l)l  I  he  situation,  however, 
is  more  complicated  in  the  ease  of  two-phase  materials, 
where  the  macroscopic  properties  ire  strongly  influenced 
by  the  properties,  the  shape,  the  volume  fraction  anil  the 
orientation  of  each  phase 

In  a  recent  study  an  ultrasonic  method  has  been 
developed  for  the  nondestructive  determination  of  texture 
m  metal  matrix  composites'  Measurements  of  ultrasonic 
veloiities  have  been  used  to  determine  the  three  fourth- 
order  <  >|  )|  -expansion  coellicients  m  SO1)!.  ’064  and  6061 
aluminum  matrix  composites  I  lie  examined  composites 
.tie  reinforced  with  up  to  20".,  silicon  carbide  t  SiC  ) 
particles  ot  about  '  4  :;m  in  diameter  and  ol  approximately 
spherical  shape  I  ndcr  the  assumptions  o|  homogeneous 
distribution  of  the  Si(  particles  and  that  the  A I  crystallite' 
and  the  Si<  particles  are  subiccted  to  the  same  stress,  the 
expansion  coellicients  offo-arth -order  have  been  determined 
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as  a  function  of  the  volume  percentage  ol  SiC  -reinforcement, 
and  found  to  lie  within  given  bounds. 

In  this  study  the  fourth-order  expansion  coellicients 
obtained  by  this  technique  are  used  to  examine  relationships 
between  texture  and  Young's  modulus  in  the  SiC  reinforced 
KIWI  and  7064  aluminium  composites.  Measurements  of 
ultrasonic  velocities  are  used  to  determine  the  Young's 
modulus  m  different  directions  of  the  composites. 

Ultrasonic  texture  analysis 

Orientation  distribution  function 

Texture  is  described  by  the  (ODE),  which  determines 
the  probability  of  finding  a  single  crystal  in  the  poly¬ 
crystalline  sample  with  a  certain  orientation  with  respect 
to  the  sample  orientations,  as  show  n  in  I'ii’iirc  I.  According 
to  Bunge’  the  ODE  can  be  written  as  a  series  expansion 
into  symmetrical  generalized  spherical  harmonies  as 

.  Mill  Vi!> 

ft./ 1-  XII  <T  7  III 

i  il  „  !  ■  I 

where  </  is  the  orientation,  the  7 are  symmetrical 
generalized  spherical  harmonies  and  the  are  the 
expansion  coellicients  The  upper  limits  A/l/l  and  \i/i 
depend  on  /  as  well  as  the  crystal  and  sample  symmetry. 

respectively.  The  orthonormal  function  system  TV  is 
invariant  towards  all  rotations  of  the  sample  symmetry 
(indicated  bv  the  right  dot  I  and  the  crystal  symmetry 
t  indicated  bv  the  lei t  two  doisi.  In  liist  approximation 
only  the  three  lourth-order  expansion  cocllincnts  (4:. 

(  ,  and  (  j '  need  to  be  considered  for  texture  evaluation 
in  cubic  materials  with  orthorhombic  sample  symmetry  '  4 
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I  he  modified  relationships  which  allow  the  determina¬ 
tion  ol  the  lourth-order  expansion  coeHieients  C4‘  for  the 
SiC  partiele-reintoreed  aluminium  matrix  composites 
are  described  in  detail  in  reference  1 .  With  c  =  c, ,  c, , 

2i  44.  w  here  <  .  and  <  44  are  the  elastic  constants  of  the 
aluminum  single  crystals.  the  expansion  coeHieients  (  4’ 
are  related  to  the  ultrasonic  velocities  T-""'1’  as  fouous 


(2) 


<  V 

210  3 

4\  -V  \  7 

t  4  ->  i  1,1  ^ 

tl  i  i  -«44»  j-  ;  S  :  -  . 
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( 3 1 

(■V 

210  3 

\  Me  y  7 
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•  o  -  2. 
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44i  -. ,  , ,  i , 

If,  4  I  :  +  1  )  , 
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w  her 

c 

t  r„i 

-I  1 

'  1 1  y;m,,i 

Mhsmp  \  ■  (  ->  s 

M,  1 

(5) 

M,_-i 

(  ' 

'  1  l 

»’  '/,  > 

16) 

and 

\E 

l‘  s.i  -  1  "s,< 

f  /s„  l'\,  /  -  1.  2.  3 

1  7 ) 

1  p  Is 

the  densitv  and  M  is  the  clastic  modulus) 

In  these  evpiations  \  designates  the  volume  fraction  of 
Si(  .  ii  v.  /  M.  «v  and  are  the  l  ame  constants  for  \l 
and  SiC  .  respectivelv.  Also  in  the  calculation  of  the 
coeHieients  C4'.  the  density  of  the  composite  can  be 
determined  according  to  the  law  of  mixture  g  „np  - 

I  1  V  l/'vi  ‘  V,\„  . 

Ihc  above  cv|uations  are  derived  under  the  following 
assumptions,  the  presence  of  the  SiC  particles  changes 
both  the  elastic  properties  of  the  composite  and  the 
texture  of  the  matrix  the  Si(  particles  are  randomly 
distributed  in  the  specimen  without  any  preferred 
orientation  I  he  ultrasonic  wavelength  is  much  larger 
than  the  average  particle  dimension  (  2  4  utn  I  so  (hat  no 
dispersion  eilecl'  mat  occur. 

Relationship  between  texture  and  Young  's  modulus 
Ihc  direction  dependence  of  Young's  modulus  I  in 
materials  w  it h  cubic  crystal  structure  containing  textures 
with  orthorhombic  svmmctrv  can  be  described  by 

/  i  o i  ,  ■  (  .  cos  2d  •  (  ,  cos  40  I S | 


where  0  is  the  angle  to  the  deformation  direction  trolling 
or  extrusion )  coinciding  with  the  I -direction  of  the  sample 
fixed  coordinate  system  t  Figure  1 1.  The  quantities  c,.  p, 
and  e  ,  in  equation  ( S )  are  constants  determined  from 
equation  t  10). 

According  to  Bunge  '  \  equation  I X )  can  also  be 
expressed  as 

EMI  I  =  E,  +  +  E,  (  4J  cos  20  +  T,C'0  cos  40 

(9) 

where  E,  is  the  Voigt  average  of  Young  s  modulus  for 
random  orientation  distribution  and  Ex.  E ,  and  E  ,  are 
combinations  ol  the  elastic  constants  of  the  single  crystals, 
liquation  I  8)  will  then  provide  the  expressions 


e,  =  I  4[  Ell)  i  +  2E{4 5  l  f  /it 90  |]  =  Em  I  l()a) 

e,  -  I  2[  /.(<)  )  -  El 90  |J  =  /.x  |  |()b) 

e,  =  1  4f  El  0  I  -  2/2(45  )  4-  Em  I]  =  I  2AE  I  Kiel 
and  from  equation  (9)  we  can  obtain 

=  /  ,( -I1  +  E,  Mia) 

/2v  =  /2,r;:  (Mb) 

A/2  =  /•.,(]'  (lie) 


Experimental 

T  he  metal  matrix  composites  examined  in  this  investiga¬ 
tion  are  the  aluminium  alloys  8001  and  7064  containing 
up  to  20" „  volume  fraction  of  silicon  carbide  particles. 
The  chemical  compositions  of  these  alloys  and  the  volume 
percentages  of  the  SiC  reinforcement  arc  shown  in 
table  I.  The  composites  were  received  as  extruded  rods 
of  25  mm  in  diameter.  Prismatic  specimens  15  to  18  mm 
wide  and  30  mm  long  were  machined  out  of  these  rods 
such  that  the  length  direction  is  parallel  to  the  extrusion 
direction  Opposite  laces  are  machined  Hat  and  parallel 
to  w  ithin  *  0.025  mm  on  these  specimens. 

In  order  to  determine  the  ODf  expansion  coeHieients. 
ultrasonic  velocity  measurements  are  performed  using  the 
pulse-echo-overlap  method,  which  is  described  in  detail 
bv  Papadakis1"  The  system  used  in  this  investigation 
consists  of  an  ultrasound  apparatus  which  generates 
pulses  of  approximate!)  I  ;/s  duration  of  a  variable 
repetition  rate.  T  hese  pulses  are  impressed  on  a  commercial 
transducer  of  a  fundamental  frequency  of  5  MM/  which 
is  acoustically  bonded  to  the  specimen  I  hc  reflected  r.f. 
echoes  are  received  by  the  same  transducer,  amplified 
and  displayed  on  a  screen  of  an  oscilloscope.  Two  of  the 
displayed  echoes  are  then  chosen  and  exactly  overlapped 
by  critically  adjusting  the  frequency  of  a  (  V\  oscillator 
and  the  vliv  isjon  factor  on  a  decade  div  ider  The  frequency. 
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I.  accurately  determined  by  an  electronic  counter  is 
employed  to  compute  the  ultrasonic  velocity  using  the 
relationship  I  =  sf,  where  s  is  the  path  length  of  the 
ultrasonic  wave.  The  system  is  capable  of  measuring 
changes  in  the  ultrasonic  velocity  to  an  accuracy  of  better 
than  one  part  in  10'.  while  the  inaccuracy  of  measuring 
both  the  longitudinal  and  the  shear  velocities  is  estimated 
to  he  0.3".. 

Results  and  discussion 

The  geometry  of  the  specimens  used  in  this  investigation 
allow  .  the  measurements  of  all  nine  ultrasonic  wave 
velocities  1  three  different  waves  in  three  orthogonal 
directions).  The  velocities  measured  are  listed  in  Table  2. 
In  these  measurements  the  3-direction  is  taken  to  be 
parallel  to  the  extrusion  direction  of  the  sample  (see 
Figure  I).  The  relationship  between  texture  and  Young's 
modulus,  described  by  equation  (  I  I  >.  is  valid  for  a  cubic 
crystal  structure  containing  textures  with  orthorhombic 
symmetry,  such  as  those  generated  by  rolling.  In  this 
arrangement  0  is  the  angle  to  the  rolling  direction,  which 
is  taken  to  be  parallel  to  the  1 -direction. 

Because  the  extruded  samples  under  examination  also 
exhibit  textures  with  orthorhombic  symmetry'  .  equation 
(III  can  be  applied  if  the  extrusion  direction  is  parallel 
to  the  1 -direction.  Rearranging  the  sample-fixed  coordinate 
system  in  this  way  leads  to  a  renaming  of  the  composite- 
velocities  f;""'1'  (  /able  2)  by  changing  the  subscripts  i 
and  /  such  that  I  -*2.  2 -->3.  3-*  I. 

The  velocities  which  should  be  equal  because  of  the 
rotational  symmetry  around  the  extrusion  axis  are  then 
averaged  to  provide  the  velocities  ( rV>)mp)a»..  (T'vTp)a>- 
(  I  and  I  I  V:"r,.i.  The  expansion  coefficients  Cj’ 

are  calculated  using  these  velocities  and  equations  (2)  (4). 
Their  values  are  plotted  versus  the  volume  percentage  of 
SiC  m  Heims  2.  2  and  4.  In  these  calculations  the 
following  Al-single  crvstal  constants  are  used":  c, ,  =  1 1 1 . 1 
GPa.  =  52  5  G Pa.  <44  =  32.5  GPa  for  AI-8091  [8  J; 
r,,  107.3  GPa.  .  =  60.3  GPa.  <44  =  28.(>  GPa  for 

■\  1-7064.  The  densitv  values  used  are1"  2524  kg  m  '  for 
AI-X09I  and  2864  kg  m  '  for  AI-7064. 

The  Young's  moduli  in  the  0  .  45  and  90  -directions 
are  determined  using  the  relationship 

,  ;d  d>(  3/.(  II)  +  2/d  d)) 

hUh  =■■  (12) 

/HO  +  /iid) 

where  the  I  ame  constants  /.(d)  and  fifth  are  determined 
by 

hi d t  ,> I s- 1 ii i  /(d)  t  2/i l d )  l  j  (13) 


Figure  1  Sample  fixed  coordinate  system  P  and  crystal  fixed 
system  K  for  characterization  of  orientation 


Figure  2  Expansion  coefficients  c”  of  the  MMC  specimens 
plotted  versus  volume  percentage  of  SiC  ((  direction  parallel  to 
extrusion  direction ) 


for  the  angles  d  =  0  .  45  .  90  .  EfO  )  and  62(90  I  are 
computed  using  the  velocity  data  shown  in  Table  2.  The 
shear  velocities  at  45  to  the  extrusion  direction  have 
also  been  measured  and  their  values  are  used  to  determine 
/.1 45  ).  The  values  of  /2<0  ).  62(45  )  and  62(90  )  as  well 
as  those  of  62m.  62 v  and  A62  for  the  AI-8091  and  AI-7064 
composites  are  included  in  Table  3.  The  plots  of  62m  versus 
C’i1-  62 A  versus  C42  and  A62  versus  are  also  shown 
in  Figures  5  to  7  and  indicate  linear  correlations  as  given 
by  equation  (11).  The  regression  equations  for  these 
relationships  can  then  be  written  as 
AI-8091 

6.m  =(2.238  (  4 1  +  78.17)  GPa  r  =  0.992 
6.\  =  (0.451  Gi2  0.003 1  GPa  r=  I 

A 62  =10  1 43  (  { '  0. 1 1 5 )  G  Pa  r  =  0.995 


Table  2 

Ultrasonic 

velocities  of  the  A! 

8091  end  A! 

7064  MMC  specimens 

Velocities  are  shown  in  units  of 

m  s 

8091 

8091 

8091 

7064 

7064 

7064 

6061 

6061 

•  0  • 

-  10% 

•  15% 

-  0% 

-  15% 

-  20% 

•  15% 

-  20 

SiC 

SiC 

SiC 

SiC 

SiC 

SiC 

SiC 

SiC 

V,  , 

661  2 

6922 

7050 

6238 

6731 

691 1 

6817 

7165 

V  •  ■ 

6606 

6913 

7011 

6232 

6734 

6937 

6780 

7103 

v\\ 

6564 

7043 

7168 

6194 

6920 

7108 

6869 

7207 

v,: 

3484 

25/08 

3816 

3059 

3435 

3584 

3413 

3695 

v.\ 

3477 

3702 

3805 

3066 

3440 

3567 

3418 

3701 

346  7 

3  728 

3834 

3082 

3479 

3641 

3438 

3731 

v‘  , 

3480 

3736 

3842 

3057 

3469 

3636 

3442 

3733 

v,\ 

34  72 

3  733 

3848 

3051 

3470 

3617 

3446 

3740 

l/.,' 

3482 

3736 

3844 

3062 

3469 

3640 

3435 

3774 
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Figure  3  Expansion  coefficient  C['  of  the  MMC  specimens  plotted 
versus  volume  percentage  of  SiC  (1  -direction  parallel  to  extrusion 

direction ) 


Figure  4  Expansion  coefficients  Ci*  of  the  MMC  specimens 
plotted  versus  volume  percentage  of  SiC  ( 1  direction  parallel  to 
extrusion  direction ) 
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!  he  lact  that  the  regression  lines  for  /•.  v  and  A/.  do  not 
cross  the  ordinate  at  the  orient  is  due  to  measuring  errors. 
NeuTtheless.  the  high  correlation  coeliicients  confirm  the 
xahdilx  of  the  relationships  between  texture  and  Young's 
modulus  (equations  ill)!  for  the  Al  SiC  composites 


Figure  5  Elastic  parameter  £m  plotted  versus  the  expansion 
coefficient  Cl1  fot  the  AI-8091  and  AI-7064  MMC  specimens 
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C4I2 

Figure  6  Elastic  anisotropy  parameter  fA  plotted  versus  the 
expansion  coefficient  C\ 2  for  the  Al  8091  and  AI-7064  MMC 
specimens 


examined.  It  is  important  to  note  that  in  these  relationships 
C|v  are  the  texture  coefficients  determined  for  the  AI- 
matrix.  whereas  F.m.  EA  and  AE  characterize  the  elastic 
anisotropy  of  the  whole  composite.  From  this  it  can  be 
concluded  that  the  elastic  anisotropy  of  the  composites 
is  influenced  only  b\  the  texture  of  the  Al  matrix,  which 
is  in  tu  ti  influenced  by  the  extrusion  process  and  the 
presence  of  the  SiC  particles.  Also  from  the  linear 
correlation  of  E  „  with  C{:.  shown  in  Figure  6.  it  is  seen, 
that  E\.  which  represents  the  difference  of  the  Young's 
modulus  in  the  extrusion  direction  and  the  Young's 


Table  3  Young  s  moduli  fiO  )  L  (45  i  and  6(90  )  and  the  resulting  parameters  6.,,  and  \F  for  the  Al  8091  and  Al  7064  MMC 
s|t»ir  intone  The-  values  ar»‘  shown  -n  units  of  GP<1 


8091 

•  0  1  - 

s  c 

8091 
•  10"-. 

SC 

8091 
-  lb".. 

SiC 

7064 
-  0% 

SiC 

7064 
•  15% 

SiC 

7064 
•  20% 
SiC 

t  1  t  >  i 

'  9  8 

91  1 

95  9 

72  1 
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99  1 

/  ■  1C  - 
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71  9 

90  0 
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f  i  no  , 
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90  0 

94  7 

72  0 
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72  0 

90  2 

97  5 
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0  66 

0  60 

0  06 

0  7b 

1  30 
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it 

0  66 

1  00 

0  1  6 

0  3b 

0  60 
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Figure  7  Elastic  anisotropy  parameter  \E  plotted  versus  the 
expansion  icienl  C.;  tor  the  Al  8091  and  Al  7064  MMC 

specimens 


modulus  perpendicular  to  that  direction,  is  equal  to  zero 
when  (  \'~  goes  to  /ero  Since  the  addition  of  SiC 
reinforcement  modifies  textures  during  the  extrusion 
process  |l>" .I  Sil  l',  the  texture  of  the  composite  can  he 
changed  b\  modifying  manufacturing  parameters  such  as 
extrusion  velocity  and  temperature  -  this  will  then  provide 
a  means  to  control  the  anisotropy  ol  the  composite,  and 
its  mechanical  properties. 
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IMTBODUCTION 

Metal  matrix  composites  hold  high  promises  as  engineering 
materials.  In  order  to  take  full  advantage  of  their  promising 
properties,  the  complex  nature  of  the  composites  must  be  understood. 

Some  questions  thus  arising;  hou  do  different  manufacturing  processes 
influence  the  microstructure  and  hou  can  the  mechanical  properties  of 
the  composites  be  explained  and  predicted  from  knowledge  of  their 
microstructure. 

The  elastic  properties  of  a  composite  material  depend  on  many 
parameters:  volume  fraction,  geometrical  shape,  sire  distribution, 
orientation  and  distribution  of  the  reinforcement  and  the  properties  of 
the  matrix.  Mathematical  models  have  been  developed  for  some  specific 
shapes  and  distributions  of  the  reinforcement.  Christensen  (1)  derived 
the  effective  shear  and  bulk  moduli  for  a  dilute  suspension  of  elastic 
spherical  particles  In  a  continuous  phase  of  another  elastic  material. 
Ledbetter  et  al.  (21  used  a  scattering  theory  to  explain  the  elastic 
behavior  of  a  particle  reinforced  composite  In  which  randomly  oriented 
ellipsoidal  particles  were  nonliomogeneously  distributed.  Experimental 
work  by  Lee  et  al.  (3|  showed  that  in  particle-reinforced  composites  the 
second-order  elastic  constants  increase  linearly  with  the  particle 
content.  Their  results  suggested  that  the  dominant  factor  In  Influencing 
the  anisotropy  Is  the  content  of  reinforcement.  Spies  C  Salama  (4| 
investigated  the  Influence  of  the  reinforcing  phase  on  the  texture  and 
found  that  the  fourth-order  expansion  coefficients  change  linearly  with 
the  particle  content.  Their  results  indicate  that  the  presence  of 
particles  in  the  composites  leads  to  considerable  changes  in  the  texture 
of  the  aluminum  matrix. 

The  objective  of  this  study  is  to  obtain  information  a£>out  what 
features  in  the  microstructure  of  the  composites  are  causing  their 
anisotropic  behavior.  This  was  accomplished  by  comparing  the  particle 
sire  distribution  in  three  orthogonal  directions,  with  the  ultrasonic 
velocities  measured  along  the  same  directions.  Three  series  of  SiC- 
particle  reinforced  composites  were  examined.  Two  series  comprised  of 
extruded  samples  while  the  third  consisted  of  pressed  specimens.  Each 
series  included  samples  of  different  particle  content.  Also,  three 
squeeze-cast  samples  with  different  particle/fiber  content  were 
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examined.  In  the  extruded  composites,  the  properties  along  the  extrusion 
direction  were  found  to  be  different  from  those  in  the  directions 
perpendicular  to  the  extrusion  direction.  This  behavior  is  explained  in 
terms  of  the  presence  of  texture  in  the  matrix.  The  squeeze-cast 
specimens  also  showed  an  anisotropic  behavior  which  is  induced  by 
preferred  orientation  of  the  fiber  reinforcement. 

MEASUREMENTS 

£c.scinigg.a 

The  metal  matrix  composites  (MMC)  used  in  this  investigation 
comprised  of  aluminum  alloys  as  the  matrix  material  and  either  SiC- 
particles  or  alumina  fibers  and  particles  as  the  reinforcement.  The 
specimens  which  have  A1-S061  as  matrix  were  received  as  pressed  plates, 
whereas  the  composites  containing  Al-7064  and  Al-8091  were  obtained  as 
extruded  rods.  The  alumina  reinforced  specimens  were  all  squeeze-cast. 
The  volume  fractions  of  reinforcement  in  the  composites  used  are  shown 
in  table  1. 

The  coordinate  systems  were  chosen  such  that  in  pressed  samples  the 
xi  and  xj-axes  are  in  the  plate  at  right  angles  to  each  other  and  the 
xj-axis  is  alon'-  the  compression  direction  (normal  to  the  plate)  .  In  the 
extruded  .«ie  Lmens  the  xj  and  xj-axes  are  perpendicular  to  the  extrusion 
direct!  ^..  ..J  orthogonal  to  each  other.  The  X3~axis  is  along  the 

extru  io>.  direction.  In  the  squeeze-cast  samples  the  xj  and  X2~ 
dirucf  .jns  are  in  the  fiber  rich  layers  at  right  angles  to  each  other, 
and  the  xj-direction  is  perpendicular  to  the  fiber  rich  layers. 

tlicxnatxucmxg 

The  particle  size  distribution  and  the  area  fraction  covered  by  the 
reinforcement  in  each  plane  were  estimated  in  each  specimen.  This  was 
accomplished  by  scanning  the  faces  of  the  specimen  under  an  optical 
microscope  and  taking  photographs  at  several  "representative"  locations 
along  the  three  chosen  directions.  The  particle  size  distributions  were 
obtained  from  these  micrographs  by  counting  the  particles  and  estimating 
their  size.  The  size  of  a  particle  was  estimated  using  its  projected 
area  on  the  face  which,  in  turn,  is  equal  to  the  area  of  a  circle  having 
the  diameter  d.  The  particle  was  considered  to  belong  to  the  size  range 


Table  1.  Metal-Matrix  Composites  used  in  investigations 


Manufacturing  method 

%  Reinforcement 

in  Compos 

ites 

Pressed  plates 

Al-6061 

0% 

SiC 

Al-6061 

251 

Sic 

Al-6061 

+ 

40% 

SiC 

Extruded  rods 

Al-7064 

+ 

0% 

SiC 

Al-7064 

+ 

15% 

SiC 

Al-7064 

+ 

20% 

SiC 

Al-8091 

0% 

SiC 

Al-8091 

+ 

10% 

SiC 

Al-8091 

+ 

15% 

SiC 

Squeeze-cast 

Al-Si-Cu-Ni-Mg 

20% 

A1203 

(Mat 

A) 

Al-Si 

+ 

20% 

ai2o3 

(Mat 

B) 

Huttenaluminum 

+ 

20% 

Al203 

(Mat 

C) 
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where  the  inequalities  n  <  d  <  n+l  are  satisfied,  where  n  and  n+1  are 
the  lower  and  upper  limits  of  the  size  range.  Also  from  the  micrographs, 
features  like  orientation,  shape  and  distribution  of  the  reinforcement 
were  examined. 

Ultrasonic. velocities 

Measurements  of  the  ultrasonic  velocities  were  performed  using  the 
pulse-echo-overlap  method,  which  is  described  in  detail  elsewhere  (51.  A 
pulse  of  approximately  l-(i3ec  duration  of  variable  pulse-repetition  rate 
is  generated  and  impressed  on  a  transducer  which  is  acoustically  bound 
to  the  specimen.  The  reflected  echoes  are  received  by  the  same 
transducer,  amplified,  and  displayed  on  an  oscilloscope.  Two  of  the 
displayed  echoes  are  then  chosen  and  exactly  overlapped  by  critically 
adjusting  the  frequency  of  the  cw  oscillator.  This  frequency  f  is 
employed  to  compute  the  ultrasonic  velocity  using  the  relation  V-21f, 
where  1  is  the  thickness  of  the  specimen.  X-  and  Y-cut  transducers  of  10 
and  2.25  MHz  were  used  for  the  generation  of  the  longitudinal  and 
transverse  waves  respectively. 


RESULTS  AND  DISCUSSION 

The  particle  size  distribution  and  the  fiber/particle  content, 
estimated  from  optical  micrographs,  ate  shown  in  figures  1-4.  Table  2 
and  3  contain  the  area  fractions  covered  by  the  reinforcement.  The  data 
are  accurate  to  within  101  of  the  nominal  values.  Table  4,5  and  6  give 
the  ultrasonic  longitudinal  and  shear  velocities  measured.  The 
velocities  are  denoted  where  i  and  j  are  the  directions  of 

propagation  and  polarization  respectively.  The  velocities  are  found  to 
be  reproducible  to  within  0.51. 

From  figures  1-3  one  finds  that  in  the  pressed  samples,  the 
particle  size  distributions  as  well  as  the  area  fractions  covered  by  the 
reinforcement  are  the  same  in  the  three  directions  within  the  accuracy 
of  the  measurement.  On  the  micrographs  the  reinforcement  showed  no 
features  explaining  the  differencies  in  ultrasonic  velocities  measured. 
This  suggests  that  the  anisotropy  is  due  to  the  te  ture  in  the  aluminum 
matrix.  The  similar  anisotropic  behavior  in  the  specimen  without 
reinforcement  also  confirms  this  statement.  However,  the  presence  of 
SiC-particles  is  found  to  vastly  enhance  the  anisotropy. 


Table  2. 

Area  fractions  covered  by  the 

SiC-particle 

reinforcement . 

Direction 

6061  + 

6061  + 

7064  + 

7064  + 

8091  + 

8091  + 

251  SIC 

1511  SIC 

101  sic 

15%  SiC 

1 

26.3 

37.7 

20.6 

23.9 

12.7 

17.5 

2 

22.6 

36.3 

18.7 

23.7 

12.1 

16.5 

3 

24.0 

31.7 

16.5 

21.0 

11.8 

15.3 

Table  3.  Area  fractions  covered  by  the  alumina  fiber  and  particle 
reinforcement  (1) . 


Direction  Material  A  Material  B  Material  C 

- Particle  Fiber  Total  Part  tel eFlber  Total  Particle  Fiber  Total 


1 

12.0 

11.9 

23.9 

8.5 

12.6 

21.1 

1.0 

13.2 

14.2 

2 

20.4 

2.5 

22.9 

15.2 

2.7 

17.9 

12.5 

3.2 

15.7 

3 

18 . 7 

2.9 

21.6 

15.7 

2.0 

17.7 

13.7 

2.7 

16.4 
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Table  4.  Ultrasonic  velocities  in  SiC  particle-Al  pressed 
composites  (m/s) 


Velocity 

6061  + 

6061  + 

6061  + 

01  SiC 

25%  SiC 

—  1  1  ■  II  ■ 

Vu 

6407 

7224 

6058 

V22 

6422 

7287 

8084 

6358 

6979 

7841 

V)J 

v2l 

3098 

3785 

4507 

3103 

3798 

4505 

V23 

3200 

3676 

4398 

V32 

3195 

3692 

4412 

Table  5.  Ultrasonic  velocities  In  SIC-particle-Al  extruded 
composites  (m/s) 


Velocity 

7064  + 

0%  SIC 

7064  + 
15%  SiC 

7064  + 

20%  SiC 

8091  + 

0%  sic 

8091  + 
10%  SIC 

8091  + 

Vll 

6251 

6728 

6935 

6617 

6890 

7002 

Vi3 

6233 

6733 

6875 

6611 

6903 

7000 

VJ3 

6250 

6902 

7071 

6626 

7025 

7141 

V)3 

3077 

3457 

3591 

3511 

3723 

3827 

Vn 

3069 

3485 

3667 

3498 

3738 

3848 

V21 

3075 

3448 

3593 

3507 

3719 

3818 

3095 

3477 

3625 

3496 

3731 

3B34 

v3l 

3056 

3490 

3666 

3492 

3745 

3856 

Vj2 

3090 

3469 

3595 

3494 

3741 

3852 

Table  6.  Ultrasonic  velocities  in  alumina  particle  and  fiber-Al 
squeeze-cast  composites  (m/s) 

Velocity 

Composite  A 

Composite  B 

Composite  C 

Vll 

7032 

6967 

6670 

7033 

6969 

6712 

VjJ 

6904 

6835 

6524 

V„ 

3694 

3656 

3419 

Vll 

3669 

3578 

3336 

Vjl 

3681 

3657 

3439 

V23 

3683 

3594 

3345 

V31 

3644 

3560 

3325 

V32 

3675 

3606 

3329 

In  the  two  series  containing  extruded  specimens  (Al-7064  and  Al- 
0091)  the  reinforcement  showed  the  same  features  in  the  three  directions 
(compare  figures  2  and  3) .  Nevertheless,  in  the  reinforced  samples  the 
longitudinal  velocities  are  higher  in  the  extrusion  direction,  whereas 
the  velocities  are  the  same  in  the  samples  without  reinforcement.  This 
behavior  further  Indicates  that  the  anisotropy  is  caused  by  texture  in 
the  aluminum  matrix  but  also  indicates  that  the  presence  of  the  SiC- 
particles  enhance  the  formation  of  texture. 

In  the  squeeze-cast  specimens  (composites  A,B,C)  the  ultrasonic 
longitudinal  velocities,  given  in  table  6,  are  considerably  higher  in 
the  plane  of  the  fibers  than  in  the  directions  perpendicular  to  that 


1444 


Number  of  particles  Number  of  particles  Number  of  particles 


300 


Direction  1 


1-2  2-3  3-4  4-5  5-6 

Particle  sire  (microns) 


Direction  1 


1-2  2-3  3-4  4-5  5-6 

Particle  size  (microns) 


Particle  size  (microns)  Particle  size  (microns) 


Fig.i  Particle  size  distribution  Fig. 2  Particle  size  distribution 

in  25%  SiC-6061  A1  in  201  SiC-7064  Al 
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Fig. 3  Particle  sire  distribution  Fig. 4  Area  fractions  covered  by  the 
in  151  SiC-8091  A1  reinforcement  in  composites 

A,B  and  C 


plane  and  the  higher  the  fiber  content  the  higher  the  difference  in 
velocity.  Ho  difference  in  longitudinal  velocity  was,  however,  found  for 
waves  propagating  in  the  plane  of  the  fibers  having  different 
polarization  directions.  Furthermore,  the  velocities  of  the  shear  waves 
propagating  normal  to  the  plane  of  the  fibers  but  polarized  in  different 
directions  are  the  same  whereas  the  velocities  of  the  shear  waves 
propagating  in  the  plane  of  the  fibers  are  higher  when  the  waves  are 
polarized  in  the  plane  than  when  the  waves  are  polarized  normal  to  the 
plane.  This  means  that  the  squeeze-cast  specimens  show  a  transversely 
Isotropic  behavior  'iich  is  in  agreement  with  their  microstructure. 
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Abstract.  A  model  has  been  developed  to  predict  the  elastic  moduli  in  composites 
reinforced  with  both  particles  and  fibers.  In  the  model  the  matrix  material  and  the 
particles,  which  are  assumed  to  be  homogeneously  distributed,  form  an  effective  ma¬ 
trix.  The  characteristics  of  this  effective  matrix  is  calculated  using  a  theory  formulated 
by  Ledbetter  and  Datta.  The  effective  matrix  is  then  considered  to  be  reinforced  with 
fibers  lying  in  one  plane  but  randomly  oriented  in  that  plane.  The  effect  of  the  2- 
dimensionally  random  orientation  of  the  fibers  on  the  elastic  moduli  of  the  composites  is 
determined  in  two  steps.  First  the  composite  cylinders  model  by  Hashin  and  Rosen  for 
an  aligned  fiber  system  is  employed,  and  then  a  geometric  averaging  procedure  sug¬ 
gested  by  Christensen  and  Waals  is  performed.  Using  this  model,  the  Young’s  and  shear 
moduli  were  calculated  for  three  samples  with  different  aluminum  matrices  and  volume 
fractions  of  particles  (9,  13,  and  17%)  but  the  same  fiber  content  (6%).  The  same  elastic 
moduli  were  also  determined  using  ultrasonic  velocity  measurements.  The  agreement 
between  calculated  and  measured  elastic  moduli  is  found  to  be  very  good.  Also,  the 
elastic  anisotropies  between  directions  of  the  fiber  rich  plane  and  that  normal  to  the 
plane  could  be  predicted  by  the  model. 

Introduction 

Many  models  have  been  developed  to  determine  the  effective  elastic  moduli  of 
composite  materials  [1-7].  Most  of  these  models  deal  with  reinforcements  in 
the  form  of  spherical  particles  [1,  2],  ellipsoidal  inclusions  [3,  4]  or  infinitely 
long  fibers  [5,  6].  However,  in  industrial  applications  the  composites  used  are 
often  of  a  more  complex  nature,  where  a  mixture  of  particles  and  fibers  is 
used  as  a  second  phase.  To  our  knowledge  no  models  that  describe  the  elastic 
properties  of  these  composites  are  available. 

Ledbetter  and  Datta  [1]  used  a  multiple  scattering  theory  to  predict  the 
elastic  behavior  of  composites  with  a  nonhomogeneous  particle  distribution.  In 
the  model  they  assume  that  the  particles  together  with  the  matrix  form  an 
enriched  “sea”  that  surrounds  “islands”  of  pure  matrix  material.  These  non- 
spherical  islands  are  aligned  and  produce  anisotropy.  The  elastic  constants 
predicted  by  the  model  are  found  to  agree  with  those  determined  by  ultrasonic 
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velocity  measurements.  Also,  the  elastic  anisotropies  in  these  composites 
could  be  explained  in  terms  of  the  nonhomogeneous  distribution  of  the 
particles. 

On  the  other  hand,  in  order  to  determine  the  effective  elastic  properties  of 
fiber  reinforced  materials,  Hashin  and  Rosen  [5]  introduced  the  composite 
cylinders  model.  In  this  model,  the  composite  is  considered  to  be  comprised  of 
infinitely  long  circular  cylinders  embedded  in  a  continuous  matrix  phase.  Each 
fiber  has  a  radius,  a ,  which  is  surrounded  by  an  annulus  of  matrix  material  of 
radius  b ,  and  the  ratio  a/o  is  considered  to  be  constant  for  ail  composite 
cylinders.  In  order  to  obtain  a  volume  filling  configuration,  the  absolute  size  of 
the  cylinders  must  vary  considerably.  Hence,  the  model  is  expected  to  provide 
reasonable  agreement  with  experimental  measurements  only  if  the  size  distri¬ 
bution  is  wide  or  the  fiber  concentration  is  low. 

The  present  study  is  concerned  with  composites  reinforced  with  both  parti¬ 
cles  and  fibers.  In  these  Composites,  the  presence  of  the  homogeneously  dis¬ 
tributed  particles  improve  the  overall  mechanical  and  thermal  properties, 
whereas  the  fibers,  randomly  oriented  in  one  plane,  introduce  the  directionality 
often  desirable  in  some  applications.  Due  to  the  mixing  of  particles  and  fibers  a 
model  determining  the  elastic  moduli  of  the  composites  is  expected  to  be  ex¬ 
tremely  extensive  and  complicated.  In  the  composites  investigated,  however, 
the  volume  fraction  of  the  reinforcement  is  relatively  low,  and  hence  we  as¬ 
sumed  no  interaction  between  the  fibers  and  the  particles.  This  assumption 
simplified  the  calculations  significantly  since  the  effects  of  fibers  and  particles 
on  the  elastic  properties  could  be  evaluated  independently. 

In  this  approach  we  first  considered  the  matrix  material  and  the  particles  to 
form  an  effective  matrix.  Since  the  particles  are  homogeneously  distributed  in 
the  metal,  the  effective  matrix  is  considered  to  be  homogeneous.  The  effective 
matrix  was  then  considered  to  be  reinforced  with  the  fibers  which  are  randomly 
oriented  in  one  plane.  The  influence  of  the  fibers  on  the  elastic  moduli  of  the 
composites  was  then  determined  first  by  using  the  composite  cylinders  model 
for  an  aligned  fiber  system  [5]  and  second  by  performing  a  geometric  average 
procedure  which  takes  care  of  the  2-dimensionally  random  orientation  of  the 
fibers.  A  good  agreement  was  obtained  between  the  calculations  and  the  experi¬ 
ment. 


Theory 

Composite  with  Spherical  Particles 

First,  we  consider  a  system  that  consists  of  an  aluminum  matrix  and  alumina 
particles.  To  find  the  elastic  properties  of  this  system,  the  results  obtained  by 
Ledbetter  and  Datta  [1]  are  adopted.  Using  a  multiple  scattering  approach 
these  authors  calculated  the  effective  wave  speeds  of  plane  waves  in  a  compos¬ 
ite  with  randomly  distributed  particles.  In  their  calculations  they  made  the 
assumptions  that  the  wavelength  of  the  ultrasonic  waves  is  long  compared  to 
the  dimensions  of  the  particles,  all  particles  have  the  same  shape  and  size,  and 
the  concentration  of  particles  in  the  composite  is  dilute.  The  expressions  for  the 
effective  bulk,  K and  shear  moduli,  Gtff,  they  obtained  are 
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Keff 

Kp 


Geff  -  Gm 
Gp  -  Gm 


CP  2  Tnj 


1  _  Cp  ( 


KJ 


3  Km  +  4G„ 


2  Tm 


c  Gp  Gm  T,jij  2  Tujjj 


_6Cj,(Gp-Gm)(Km  +  2G„ 

e- 


25 


3  Km  +  4G„ 


(i) 


(2) 


where  the  subscripts  m  and  p  refer  to  matrix  and  particles  respectively,  and 
Cp  is  the  volume  fraction  of  particles.  The  expressions  for  Tnil  and  7„,y  are  given 
in  [II. 

The  system  consisting  of  the  aluminum  matrix  and  the  alumina  particles  is 
now  considered  as  the  effective  matrix.  In  order  to  obtain  the  elastic  properties 
of  the  composites  under  consideration,  the  characteristics  of  a  composite  com¬ 
prised  of  the  effective  matrix  and  fibers  randomly  oriented  in  a  plane  are 
needed.  The  first  step  in  finding  these  characteristics  is  to  calculate  the  proper¬ 
ties  of  a  composite  with  fibers  aligned  in  one  direction. 


Composite  with  Aligned  Fibers 

The  most  common  model  used  for  this  geometry  is  the  composite  cylinders 
model  introduced  by  Hashin  and  Rosen  [5],  In  this  model,  the  fibers  are  as¬ 
sumed  to  be  infinitely  long  circular  cylinders  embedded  in  a  continuous  matrix 
phase.  With  every  fiber  of  radius,  a,  there  is  anvanmilus  of  matrix  material  of 
radius,  b,  associated  with  it.  The  ratio  of  these  radii  alb  is  considered  to  be 
constant  for  all  composite  cylinders  but  the  absolute  values  of  the  radii  a  and  b 
vary  such  that  a  volume  filling  configuration  is  obtained.  The  effective  Young's 
modulus  £n  .  the  Poisson's  ratio  vt2.  the  plane  strain  bulk  modulus  K2 3  and  the 
shear  moduli  Gl2  and  G25  of  the  composite  can  then  be  expressed  as 


£11  -  CfEf  +  (1  -  Cf)Em  +  4C/(I  -  Cf)Gn 


(1  ~  Cf)  Gm  +  CfGm 


Kr  +  GJ  3  Km  +  Gm/3 


+  1 


Cf(  1  -  C,){vf  -  vm) 


v,2  =  CfVf  +  (1  -  Cf)v„  + 


K„  +  GJ 3  Kf  +  G/3 


K2  3  =  Km  +  ^  + 


( 1  -  Cf)  Gm  Cf  Gm 
Kf  +  Gfl  3  Km  +  GJ  3  + 

Cf _ 


1 


1  -  Cf 


Kf-  Km  +  (Gf-  Gm)l 3  Kn 
Gn  =  Gr(l  +  Cf)  +  GJ  1  -  Cf) 

Gm  Gf(\  -  Cf)  +  GJ\  +  Cf) 


4Gm/3 


=  1  + 


d- 


Gm  (K„  +  7GW/3)(1  -  Q) 
Gf  -  Gm  +  2 (Km  +  4Gm/3) 


(3) 
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where  the  subscripts  m  and  /  refer  to  matrix  and  fibers  respectively,  and  Cf  is 
the  volume  fraction  of  fibers.  From  these  relationships  the  five  independent 
components  of  the  second  order  elastic  constant  tensor  can  be  calculated  using 
the  relationships  [7] 


C M  =  £n  +  4i^2  K2 3 
C|2  =  2KziV\2 

Cn  =  <J23  +  A*23 
C23  =  -G23  +  K2] 

C66  =  Gn  (4) 

The  above  expressions  give  relationships  for  the  elastic  properties  of  a 
composite  reinforced  with  spherical  particles  and  aligned  fibers.  The  character¬ 
istics  of  a  composite  where  the  fibers  are  randomly  oriented  in  a  plane  are  then 
to  be  determined. 

Composite  with  2-D  Randomly  Oriented  Fibers 

A  schematic  representation  of  the  problem  is  shown  in  Fig.  1.  The  figure  shows 
that  all  the  fibers  are  lying  in  the  jr^-plane  but  randomly  oriented  in  that  plane. 
To  determine  the  elastic  properties  of  this  configuration,  the  approach  of  Chris¬ 
tensen  and  Waals  [7]  is  used.  In  this  approach,  the  effect  of  a  random  orienta¬ 
tion  of  fibers  on  the  stress  to  strain  ratios  cr'/ej  Is  analytically  equivalent  to 
finding  the  average  value  of  the  ratio  aiJej  when"  all  possible  orientations  are 
taken  relative  to  a  fixed  axis  as  shown  in  Fig.  2.  In  the  2-dimensional  case  this 
can  be  expressed  as 


random 


\_ 

77 


(5) 


where  the  stress  to  strain  ratio  in  the  integrand  refers  to  the  aligned  fiber  system 
shown  in  Fig.  2  and  #  is  the  angle  between  the  x\-  and  the  xt-axis.  The  stress  to 
strain  ratios  are  then  calculated  using  the  tensor  transformation  laws  for 
stresses  and  strains 
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Fig.  1.  Composite  with  fibers 
randomly  oriented  in  the 
x,x2-plane. 


Fig.  2.  Composite  with  fibers 
aligned  along  the  x,-direction. 
The  X,-  and  x2-axes  are  ro¬ 
tated  an  angle  6  with  respect 
to  the  Xr  and  xj-axes. 
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and  the  stress-strain  relationships  for  a  system  with  aligned  fibers  which  are 
given  by 


(T\  =  C||E|  +  C|2E2  +  C|2E3 

(X  2  =  C|2E|  +  C22E2  +  C  23E3 

O-3  =  C|2E|  +  C23E2  +  C22E3 
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(8) 


The  stress  to  strain  ratio  <T\le\  can  then  be  evaluated  by  imposing  a  strain  along 
the  1 '-direction  such  that  the  conditions 


£1^0 

ei  =  0.  i  *  1  (9) 

are  satisfied.  The  transformation  laws.  Eqs.  (6)  and  (7),  and  the  stress-strain 
relationships,  Eq.  (8).  then  yield 


a\  =  cr,  cos:0  +  <r2  sin:0  -  cr6  sin20 
— r  =  C 11  cos:0  +  Cr  sin20 

E| 

=  C|>  cos:0  +  C22  sin20 
E] 

77  =-(:<*  sin ze  (10) 

£  1 

Combining  Eqs.  (10)  results  in 


-  Cn  cos40  +  2Ci2  sin:0  cos;0  +  On  sinJ0  +  40*,  sin:0  cos:0  (11) 
Integration  according  to  Eq.  (5)  gives 

Q 1.  =  (77)  ,  =5  (3C„  3-  2C|i  +  30:  +  4C*)  (12) 

The  stress  to  strain  ratio  evaluated  in  Eq.  (12)  will  correspond  to  one  of  the 
components,  (2n  •  in  the  second  order  elastic  constant  tensor  for  the  composite 
in  Fig.  1 .  The  other  components  in  this  tensor  can  be  evaluated  similarly  and  are 
found  to  be 
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(13) 


For  the  composite  with  2-D  randomly  oriented  fibers,  these  five  components 
fully  describe  the  stress-strain  relationships  which  are  given  by 
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From  the  five  components  of  the  second  order  elastic  constant  tensor,  the 
effective  Young's  and  shear  moduli  of  the  composite  can  be  calculated  using 
the  relationships 
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Measurements 


(15) 


Specimens 

The  metal  matrix  composites  used  in  this  investigation  are  comprised  of  an 
aluminum  alloy  as  the  matrix  nr'  '  and  alumina  particles  and  fibers  as  the 
reinforcements.  The  three  aluminum  alloys  used  as  matrix  materials  consist  of 
Al-129c  Si.  Al-12%  Si-Cu-Ni-Mg,  and  commercially  pure  aluminum. 

The  composites  were  manufactured  using  the  squeeze  casting  method  and 
produced  in  the  form  of  bars.  From  these  bars,  blocks  measuring  10  x  12  x  20 
mm  were  cut.  The  blocks  were  machined  and  ground  such  that  opposite  faces 
were  flat  and  parallel  to  within  10  ju.m. 

The  coordinate  systems  for  the  specimens  were  chosen  such  that  the  xt- 
and  .r;-  directions  are  in  the  fiber  rich  layers  at  right  angles  to  each  other  and  the 
xi-direction  is  perpendicular  to  the  fiber  rich  layers.  The  volume  fractions  of 
particles  and  fibers  in  these  samples  were  estimated  using  optical  microscopy. 
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Ultrasonic  velocities 

Measurements  of  ultrasonic  velocities  were  performed  using  the  pulse-echo- 
overlap  method,  which  is  described  in  detail  elsewhere  [8].  A  pulse  of  approxi¬ 
mately  l-/xsec  duration  of  variable  pulse-repetition  rate  is  generated  and  im¬ 
pressed  on  a  transducer  that  is  acoustically  bound  to  the  specimen.  The 
reflected  echoes  are  received  by  the  same  transducer,  amplified,  and  displayed 
on  an  oscilloscope.  Two  of  the  displayed  echoes  are  then  chosen  and  exactly 
overlapped  by  critically  adjusting  the  frequency  of  the  cw  oscillator.  This  fre¬ 
quency  /is  employed  to  compute  the  ultrasonic  velocity  using  the  relation  V  = 
Ilf,  where  /  is  the  thickness  of  the  specimen.  X-  and  F-cut  transducers  of  10  and 
2.25  MHz  were  used  for  the  generation  of  the  longitudinal  and  transverse 
waves  respectively.  The  elastic  constants  were  calculated  using  the  relation¬ 
ship 


Cpp  =  C,J,J  =  pVl  p  =  1.2 . 6.  ij  =  1,2,3  (16) 

where  p  is  the  mass  density,  and  indices  ij  are  contracted  to  index  p  according 
to  Voigt's  notation. 


Microstructure 

The  microstructure  of  one  of  the  composites  used  in  this  study  is  shown  in  Fig. 
3.  The  micrographs  show  the  microstructure  in  the  three  directions  for  the 
composite  with  hiittenaluminum  as  the  matrix  metal.  The  dark,  almost  circular 
areas  are  the  particles,  whereas  the  fibers  appear  as  light  grey  rods,  20-50  pm 
in  length.  The  micrographs  clearly  illustrate  that  the  fibers  are  lying  in  one 
plane,  with  its  normal  in  the  3-direction,  but  they  are  randomly  oriented  in  that 
plane.  Furthermore,  the  1-  and  2-directions  exhibit  the  same  features.  They  both 
show  a  low  fraction  of  fibers  oriented  in  the  X|X2-plane  and  have  the  same 
particle  content.  These  observations  confirm  that  a  transversely  isotropic  be¬ 
havior  of  the  composites  is  to  be  expected. 

The  area  fractions  covered  by  the  fiber  and  particle  reinforcements  in  each 
plane  were  estimated  in  the  three  specimens.  This  was  accomplished  by  scan¬ 
ning  the  faces  of  the  specimen  under  an  optical  microscope  and  taking  photo¬ 
graphs  at  several  "representative"  locations  along  the  three  chosen  directions. 
The  area  fractions  were  obtained  from  these  micrographs  by  adding  the  pro¬ 
jected  area  on  the  face  of  all  particles  and  fibers  and  dividing  by  the  total  area. 


Results 

Table  I  contains  the  area  fractions  of  the  fiber  and  particle  reinforcements  in 
the  composites  used  in  this  investigation  as  estimated  from  the  optical  micros¬ 
copy  study.  Table  2  includes  values  of  the  elastic  constant  found  in  literature 
[9]  for  the  aluminum  alloys  used  as  matrices  for  the  composites.  It  also  includes 
the  measured  mass  densities  of  the  composites  and  the  volume  fractions  of  the 
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FI*.  3.  Microstructure  in  three  orthogonal 
directions  of  the  composite  with  Hiitten- 
aluminum  as  the  matrix  metal:  a  direction  1 : 
b  direction  2:  c  direction  3.  The  dark  almost 
circular  areas  are  alumina  particles  whereas 
the  light  grey  rods  of  20-50  /im  in  length  are 
alumina  fibers.  Bar  =  20  *tm. 


alumina  reinforcements.  The  volume  fractions  of  particles  and  fibers  were  ob¬ 
tained  by  taking  the  average  of  the  area  fractions  in  the  three  directions. 

Table  3  lists  the  ultrasonic  velocities.  V„  ,  measured,  where  i  and  j  denote 
the  directions  of  propagation  and  polarization  respectively.  In  Table  4  the 
calculated  and  measured  elastic  constants  are  listed.  Column  2  gives  the  elastic 
constants  of  the  effective  matrix,  consisting  of  the  aluminum  matrix  and  the 
alumina  particles.  These  constants  were  calculated  using  Eqs.  (1)  and  (2). 
Using  Eq.  (3)  the  elastic  constants  of  composites  with  aligned  fibers  were 


Table  1.  Area  percentage  of  alumina  particle  and  fiber  reinforcements  in  MMC 
specimens 


Direction 

Al-Si-Cu-N 

i-Mg 

Matrix 

Al-Si 

Hiittenaluminum 

Particles 

Fibers 

Particles 

Fibers 

Particles 

Fibers 

1 

18.7 

2.9 

15.7 

2.0 

13.7 

2.7 

20.4 

2.5 

15.2 

2.7 

12.5 

3.2 

3 

12  0 

11.9 

8.5 

12.6 

1.0 

13.2 
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Table  2.  Elastic  constants,  densities,  and  volume  fractions  of  MMC 


specimens  used  in  calculations 


Material 

E  (GPa) 

V 

P  (kg/m!) 

Cp 

Q 

Al-Si-Cu-Ni-Mg 

73.1 

0.33 

2825 

0.17 

0.06 

Al-Si 

71.0 

0.33 

2760 

0.13 

0.06 

Huttenaluminum 

68.9 

0.33 

2793 

0.09 

0.06 

Alumina 

372.0 

0.224 

Table  3.  Ultrasonic  velocities  measured  in  MMC  specimens 

Matrix 

Velocity  tm  x) 

Al-Si-Cu-Ni-Mg 

Al-Si 

Huttenaluminum 

v„ 

7032 

6967 

6670 

V;: 

7033 

6%9 

6712 

V„ 

6907 

6835 

6524 

V,; 

3694 

3656 

3419 

V„ 

3669 

3578 

3336 

V:, 

3681 

3657 

3439 

V;j 

3683 

3594 

3345 

V., 

3644 

3560 

3325 

V,; 

3675 

3606 

3329 

Table  4.  Calculated  and  measured  elastic  constants  gf -MMC  specimens  with  Al-Si- 
Cu-Ni-Mg.  Al-Si.  and  Huttenaluminum  as  the  matrix  metal 


Matrix 

Elastic 

constant 

Effective 

matrix 

Aligned 

fibers 

Random 

fibers 

Measured 

Al-Si-Cu-Ni-Mg 

En 

92.2 

109  1 

102.2 

100.4 

E;; 

92.2 

99.5 

102.2 

100.4 

Ell 

92.2 

99.5 

99.6 

98.7 

G,: 

34.8 

37.5 

38.8 

38.4 

Gn 

34.8 

37.5 

37.4 

37.8 

G;i 

34.8 

37.2 

37,4 

38.2 

Al-Si 

Eh 

85.0 

102.3 

94.9 

95.0 

85.0 

92.0 

94.9 

95.4 

E„ 

85.0 

92.0 

92.1 

92.9 

O', ; 

32,0 

34.6 

35.9 

36.9 

Op 

32.0 

34.6 

34.4 

35.1 

G. , 

32.0 

34.2 

34.4 

35.8 

Huttenaluminum 

£i  i 

78.1 

95.8 

88.0 

84.6 

e:: 

78.1 

84.9 

88.0 

85.4 

Eh 

78.1 

84.9 

85.0 

81.9 

O'i: 

29.6 

32.0 

33.5 

32.8 

G,i 

29.6 

32.0 

31.9 

31.0 

29.6 

31.7 

31.9 

31.0 
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Fig.  4.  Calculated  effects  of 
changing  the  alumina  fiber  con¬ 
tent  on  the  anisotropies  in  the 
Young's  and  the  shear  moduli  of 
a  composite  containing  the  Al-Si 
alloy  as  the  matrix  metal  and  6 Vc 
alumina  particles. 


computed  and  the  results  are  listed  in  column  3.  Column  4  lists  the  predicted 
elastic  constants  when  the  fibers  are  randomly  oriented  in  one  plane.  These 
results  were  obtained  by  employing  Eqs.  (4).  (12).  (13).  and  (15).  These  values 
are  to  be  compared  with  the  measured  values  given  in  column  5.  The  measured 
data  were  computed  using  the  densities  given  in  Table  2  and  the  ultrasonic 
velocities  listed  in  Table  3  and  are  found  to  be  reproducible  to  within  1%. 

Figure  4  displays  the  calculated  effects  of  changing  the  fiber  content  on  the 
anisotropies  in  the  Young's  and  the  shear  moduli  of  the  composites.  In  these 
calculations  the  effective  matrix  consists  of  the  Al-Si  alloy  as  the  matrix  metal 
and  6%  alumina  particles.  Figure  5  illustrates  the  calculated  effect  on  the  same 
moduli  upon  changing  the  particle  content  in  the  effective  matrix.  In  these 
calculations  the  Al-Si  alloy  was  also  used  as  the  matrix  metal  and  the  fiber 
content  was  6 %. 


Discussion 

From  the  metallurgical  observations  it  is  seen  that  the  area  fraction  of  fibers  in 
the  3-direction  is  vastly  higher  than  those  measured  in  the  two  directions  in  the 
plane  perpendicular  to  the  3-direction.  In  the  I-  and  2-directions  the  fiber  con- 
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Fig.  5.  Calculated  effects  of 
changing  the  alumina  particle 
content  on  the  anisotropies  in 
the  Young's  and  the  shear  mo¬ 
duli  of  a  composite  containing 
the  A!-Si  alloy  as  the  matrix 
metal  and  6%  alumina  fibers. 


tent  is  the  same  to  within  ±  0.4%.  Also  the  micrographs  show  that  the  particle 
content  in  the  3-direction  is  significantly  lower  than  those  in  the  other  two 
directions.  However,  the  projected  area  of  the  particles  is  almost  circular  in  all 
three  directions  indicating  spherical  particles.  Due  to  this  result  and  the  tact 
that  the  manufacturing  process  should  not  produce  any  anisotropy  due  to  the 
distribution  of  particles,  it  is  assumed  in  the  calculations  that  the  particles  are 
spherical  and  that  they  are  randomly  distributed. 

The  experimental  results  in  Table  3  show  that  the  composites  exhibit  a 
transversely  isotropic  behavior  which  is  in  agreement  with  the  microstructure 
shown  in  Fig.  I.  This  behavior  requires  that  £ji  =  £3;  and  Gn  =  G21  which  is  in 
agreement  with  observations  shown  in  Table  4.  to  within  2%.  Furthermore,  the 
elastic  constants  in  the  3-direction  differ  from  those  in  the  I  -  and  2-directions 
such  that  £u  <  £11  and  Gp  <  G12 .  This  also  agrees  with  the  microstructure 
observed  in  the  three  samples  since  the  fiber  rich  plane  is  expected  to  be 
elastically  stiffer  than  the  planes  normal  to  that  plane. 

The  general  behavior  of  the  predicted  elastic  moduli  as  a  function  of  the 
fiber  or  particle  content  is  shown  in  Figs.  4  and  5  respectively.  It  can  be  seen 
from  these  figures  that  the  four  elastic  moduli  increase  almost  linearly  with  the 
fiber  content  and  as  a  consequence  the  anisotropy  exhibit  the  same  behavior. 
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As  the  particle  content  is  increased  the  elastic  moduli  deviates  slightly  from  a 
linear  relationship.  Furthermore,  £3 3  and  Gn  increase  at  a  faster  rate  than  £M 
and  Gn  such  that  the  anisotrophy  produced  by  the  fibers  is  slightly  lowered  as 
the  particle  content  is  increased. 

Since  the  fiber  content  in  the  three  samples  was  the  same,  the  anisotropy 
was  expected  to  be  more  pronounced  in  the  composite  with  the  soft  matrix  than 
that  with  a  stiffer  matrix.  As  can  be  seen  in  Table  4,  this  behavior  was  experi¬ 
mentally  observed  in  both  the  Young’s  and  the  shear  moduli.  As  the  stiffness  of 
the  effective  matrix  is  increased,  the  difference  in  the  moduli  between  the 
directions  in  the  plane  and  normal  to  the  plane  of  the  fibers  becomes  less 
pronounced.  The  quantitative  values  of  the  anisotropies  are  also  well  predicted 
by  the  model.  For  the  material  with  the  lowest  particle  content  (9%),  the  values 
agree  almost  exactly.  As  the  particle  content  gets  higher,  the  observed  aniso¬ 
tropies  deviate  from  those  predicted  but  the  agreement  is  still  good.  This  trend 
is  expected  since  the  model  assumes  no  interaction  between  the  particles  and 
the  fibers,  and  accordingly  the  model  will  be  more  applicable  to  composites 
with  low  particle  content. 

The  deviation  from  predicted  anisotropies  can  also  be  explained  from  a 
microstructural  point  of  view.  When  the  particle  content  is  increased  it  be¬ 
comes  more  difficult  to  orient  all  the  fibers  in  a  specific  plane.  In  addition,  it 
seems  that  the  presence  of  particles  will  force  some  of  the  fibers  to  be  oriented 
such  that  they  make  a  small  angle  with  their  ideal  orientation,  and  the  true 
configuration  will  differ  from  that  assumed  by  the  model.  The  anisotropy  pro¬ 
duced  will  then  be  lower  than  that  predicted. 


Summary 

The  elastic  behavior  of  metal-matrix  composites  reinforced  with  homogene¬ 
ously  distributed  particles  and  fibers  randomly  oriented  in  two  dimensions,  has 
been  investigated.  Using  a  model  which  assumes  no  interaction  between  parti¬ 
cles  and  fibers,  the  elastic  moduli  as  well  as  the  elastic  anisotropies  were 
calculated  for  the  composites  containing  9.  13,  and  17%  particles  and  6%  fibers. 
The  values  computed  by  the  model  are  found  to  agree  very  well  with  those 
determined  from  ultrasonic  velocity  measurements. 

These  results  indicate  that  for  relatively  low  volume  fractions  of  reinforce¬ 
ment  the  interaction  between  particles  and  fibers  can  be  neglected  on  determin¬ 
ing  the  effective  elastic  properties  of  these  composites.  This  simplifies  the 
calculations  significantly  when  these  composites  are  to  be  used  in  applications 
since  the  effects  of  adding  particles  that  improve  the  overall  properties,  and 
fibers  that  generate  the  directionality,  can  be  estimated  independently.  This  is 
applicable  for  volume  fractions  that  are  not  too  high. 


Acknowledgment.  This  work  is  sponsored  by  the  Army  Research  Office  under  contract  No. 
DAAL03-88-K-00% 


% 


Grelsson  and  Salama 


References 

!.  H.M.  Ledbetter  and  S.K.  Datta.  J.  Acoust.  Soc.  Am.  79:239  (1986) 

2.  B.  Budiansky.  J.  Mech.  Phys.  Solids  13:223  (1965) 

3.  J.D.  Eshelby.  Proc.  Roy.  Soc.  ( London )  241A.376  (1957) 

4  T,S.  Chow.  J.  Appl.  Phys.  48:4072  (1977) 

5.  Z.  Hashin  and  B.W.  Rosen.  J.  Appl.  Mech.  31:223  (1964) 

6.  R.  Hill.  J.  Mech.  Phys.  Solids  12:199  (1964) 

7.  R.M.  Christensen  and  F.M.  Waals.  J.  Comp.  Mater.,  6,  518  (1972). 

8.  K.  Salama  and  C.K.  Ling.  J.  Appl.  Phys.  51:1505  (1?8C) 

9.  Metals  Handbook,  vol.  1.  8th  ed.  ASM  (1961) 


R»c.  QA )DE  Mol  10,  till  in  Press 


ACOUSTIC  NONLINEARITY  IN  METAL  -  MATRIX  COMPOSITES 


H.  Mohrbacher,  D.  Lee,  E.  Schneider  *  and  K.  Salama 

Department  of  Mechanical  Engineering 

University  of  Houston 

Houston,  TX  77204,  (713)  749-4455 

*  Fraunhofer  Institute  for  Nondestructive  Testing,  Germany 


INTRODUCTION 

The  elastic  behavior  of  a  solid  consists  of  linear  and  nonlinear  contributions.  The 
linear  part  is  represented  by  the  well  known  Hooke’s  law  which  is  given  in  tensorial 
notation  as 


a. 

>j 


“kl 


(1) 


where  ajj  and  £ki  are  the  stress  and  strain  tensors,  respectively,  and  Cijki  is  the  tensor  of 
the  second  order  elastic  constants  (SOEC).  This  relationship  is  sufficient  for  many 
engineering  calculations  since  deviations  from  a  purely  linear  elastic  behavior  are  small. 
Hooke's  law,  however,  is  not  sufficient  for  an  advanced  characterization  of  the  elastic 
behavior  of  materials.  This  is  due  to  the  fact  that  many  of  the  physical  and  mechanical 
properties  of  materials  are  of  nonlinear  nature.  The  nonlinear  elastic  behavior  can  be 
investigated  using  ultrasonic  techniques  because  of  their  high  sensitivity  for  small  nonlinear 
effects.  Among  the  nonlinear  effects  are  the  stress  and  the  temperature  dependences  of 
ultrasonic  velocities  in  the  solid.  These  effects  have  gained  considerable  interest  in  the  last 
decade,  particularly  for  the  nondestructive  evaluation  of  applied  and  residual  stresses  [1], 
and  also  for  the  microstructural  characterization  of  materials  [2],  Another  physical 
manifestation  of  the  nonlinear  elastic  behavior  of  solids  is  the  acoustic  nonlinearity 
parameter.  This  parameter  can  be  determined  from  measurements  of  the  amplitudes  of 
fundamental  and  second  harmonic  when  an  originally  sinusoidal  wave  gets  distorted  while 
propagating  through  the  solid.  The  nonlinearity  parameter  can  also  be  calculated  from  a 
combination  of  second  and  third  order  elastic  constants.  In  previous  studies  [3],  the 
nonlinearity  parameter  was  found  to  be  sensitive  to  microstructural  changes  in  aluminum 
alloys  and  in  particular  to  the  content  of  precipitates  of  the  second  phase. 


Metal  -  matrix  composites  are  a  new  class  of  materials  wich  contain  a  metallic 
matrix  and  a  metallic  or  ceramic  material  as  a  reinforcement.  The  bulk  properties  of  these 
composites  can  be  tailored  by  changing  the  volume  percentage,  geometry,  distribution,  and 
orientation  of  the  reinforcement.  Due  to  the  different  coefficients  of  thermal  expansion  for 
the  matrix  and  the  reinforcement  material  in  MMCs,  the  creation  of  thermal  stresses  during 
manufacturing  is  unavoidable.  Therefore,  the  nondestructive  characterization  of  these 
composites  is  necessary  in  order  to  monitor  their  mechanical  properties  and  to  guarantee 
their  quality. 


In  the  present  study  the  ef'cct  of  the  volume  content  of  reinforcement  in  metal  - 
matrix  composites  (MMC)  on  the  two  nonlinear  elastic  quantities,  namely  acoustoelastic 


constants  and  nonlinearity  parameter  has  been  investigated.  The  nonlinearity  parameter  is 
determined  using  two  different  methods  and  the  results  are  compared. 


THEORETICAL 


Basically,  all  nonlinear  elastic  effects  are  due  to  the  anharmonicity  of  the  interatomic 
potential.  Thus,  relationships  betweei.  quantities  describing  the  elastic  nonlinearity  are 
expected.  In  order  to  develop  quantitative  parameters  for  the  description  of  elastic 
nonlinearity,  it  is  convenient  to  use  the  thermodynamic  derivation  of  the  elastic  constants 
starting  from  the  elastic  potential  of  the  solid.  If  the  lattice  arrangement  of  a  solid  is 
disturbed  by  an  infinitesimal  strain  e  due  to  the  presence  of  an  elastic  wave,  the  energy  of 
deformation  per  unit  volume  0(e)  can  be  expanded  as  a  power  series  of  strains  such  that 


0(e)  =  0  +  *>.£,.  +  -22-  £  .  +  —22 - 

0  ae..  y  dc.de..  ,J  kl  ae.a&^E 
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According  to  Brugger  [4],  the  elastic  constants  of  the  order  n  are  defined  as  the  n-th  partial 
derivatives  of  the  elastic  potential  with  respect  to  strain  as 


c..  =  a(n)o/ae.. 

Ij...  Ij. 


(3) 


If  up  to  third  order  terms  in  e  are  considered  in  eq.(2),  the  stress-strain  relationship  can  be 
written  as 


°ij  Cijkl  ^1  +  Cijklmn  Vmn 

where  the  Cijklmn  are  the  third  order  elastic  constants  (TOEC)  which  need  to  be  added  to 
Hooke's  law  (1)  to  allow  for  nonlinear  deviations. 


For  isotropic  materials,  the  tensor  of  the  second  order  elastic  constants  reduces  to 
two  independent  second  order  elastic  constants,  known  as  Lame  constants  X  and  p.  These 
constants  can  be  determined  directly  from  ultrasonic  experiments  by  measuring  the  wave 
speeds  of  longitudinal  and  shear  waves  using  the  relationships  : 

2  2 
|i  =  p  vy  and  X  +  2|i  =  p  vL  (5) 

The  tensor  of  the  third  order  elastic  constants  reduces  to  three  independent  third  order 
elastic  constants  in  isotropic  materials.  These  are  called  the  Mumaghan  constants  1,  m  and 
n.  In  order  to  measure  these  constants  using  ultrasonic  methods  the  propagation  velocities 
of  three  different  wave  modes  have  to  be  determined  as  a  function  of  an  applied  uniaxial 
strain.  The  relative  change  in  the  ultrasonic  velocity  of  a  given  wave  mode  with  the  applied 
elastic  strain  normalized  by  the  velocity  of  the  strain-free  specimen  is  called  the 
acoustoelastic  constant  (AEC)  which  is  a  characteristic  of  the  material.  The  Mumaghan 
constants  can  be  evaluated  from  acoustoelastic  constants  using  the  relationships  [6] 
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wnere  v  is  the  Poisson’s  iatio,  X  and  p  are  the  Lame  constants,  vij  is  the  velocity  of  an 
ultrasonic  wave  with  propagation  direction  i  and  polarization  direction  j,  v0  is  the  velocity 


in  the  unstrained  specimen  and  e  is  a  uniaxial  strain  applied  in  a  direction  perpendicular  to 
the  propagation  direction. 


Another  nonlinear  quantity  which  can  be  obtained  if  eq.(4)  is  inserted  into  the 
equation  of  motion  for  particles  in  the  solid's  lattice,  leads  to  the  nonlinear  wave  equation 


du.  d2u. 
1  1 


(7) 


where  u  is  the  particle  displacement,  i  is  a  mode  index  depending  on  the  polarization  and 
the  propagation  direction  of  the  wave,  a  is  a  coordinate  along  the  direction  of  wave 
propagation  and  vi  is  the  wave  velocity  of  the  mode  i.  The  quantity  pi  is  the  modal  acoustic 
nonlinearity  parameter  of  the  solid  and  can  be  expressed  as  a  linear  combination  of  the 
second  and  third  order  elastic  constants  [5].  The  solution  of  the  nonlinear  wave  equation 
can  be  given  as 


u;  =  Ai  sin  (ka  -  cot)  +  A2  cos  2(ka  -  cot) 


(8) 


where  Ai  and  A2  are  the  amplitudes  of  the  fundamental  and  second  harmonic  waves, 
respectively.  The  acoustic  nonlinearity  parameter  is  related  to  these  amplitudes  as 


P  = 


8  \ 


2  2 
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where  w  is  the  fundamental  frequency,  k  is  the  propagation  constant  and  a  is  the  distance 
measured  from  the  generating  transducer  to  the  instantaneous  position  of  the  fundamental 
wave  in  the  solid.  Its  magnitude  determines  the  extend  of  the  distortion  of  the  fundamental 
wave.  If  one  considers  a  longitudinal  wave  mode  in  an  isotropic  solid,  the  acoustic 
nonlinearity  parameter  is  related  to  the  Lame  and  the  Mumaghan  constants  by 


21  +  4m 
P  =  3  + - 

\  +  2p 


(10) 


Using  this  relationship  the  acoustic  nonlinearity  parameter  car.  be  calculated  when  the  third 
order  elastic  constants  are  known. 


EXPERIMENTAL 

In  this  study,  two  different  sets  of  specimens  were  used.  The  matrix  was  an  AI  - 
8091  alloy  in  one  case  and  in  the  other  case  an  AI  -  7064  alloy.  The  chemical  compositions 
of  both  alloys  are  shown  in  table  1.  As  reinforcement  material  these  composites  contain  SiC 
-  particles.  The  particles  had  a  more  or  less  globular  shape,  ranging  from  1  -  5  pm  in  size. 
The  material  has  finally  been  extruded  to  rods.  Micrographs  taken  for  different  cuts  of  the 
specimens  revealed  a  planar  random  distribution  of  the  reinforcement  in  the  plane  normal  to 
the  extrusion  direction  and  a  particle  alignment  along  the  extrusion  direction. 

A  block  diagram  of  the  experimental  set  up  for  the  determination  of  the  nonlinearity 
parameter  is  shown  in  Fig.  1  .  A  lithium  niobate  transducer  attached  to  the  specimen  by  a 
solid  bond  is  used  to  generate  ultrasonic  pulses.  These  pulses  have  typically  a  center 
frequency  of  10  MHz  and  a  bandwidth  of  200  kHz.  The  ultrasonic  signal  propagates  along 
the  extrusion  direction  of  the  specimen  whose  surfaces  are  made  parallel  to  each  other  and 
are  lapped  optically  flat.  The  10  MHz  fundamental  as  well  as  the  20  MHz  harmonic  signal 
cause  a  distortion  of  the  free  surface  of  the  specimen.  The  displacement  amplitudes  carried 
by  the  two  frequencies  are  measured  in  order  to  determine  the  nonlinearity  parameter  using 
eq.(9). 


Table  1:  Chemical  compositions  of  the  aluminum  alloys  in  weight  percent . 


Alloying  Elements 


Alloy 

Si 

Fc 

Cu 

Mg 

Zr 

Li 

Zn 

Cr 

Co 

Ai 

8091 

0.02 

0.01 

1.90 

0.80 

0.11 

2.70 

.... 

.... 

.... 

rem 

7064 

0.05 

0.10 

2.00 

2.30 

0.20 

7.10 

0.12 

0.22 

rem 

For  measurements  of  the  absolute  amplitudes,  the  capacitive  detector  technique 
described  in  [8]  has  been  used.  This  technique  allows  the  detection  of  displacement 
amplitudes  of  a  free  surface  with  a  sensitivity  of  103  A. 


In  order  to  verify  the  behavior  predicted  by  eq.(9),  the  amplitudes  of  the 
fundamental  and  second  harmonic  have  been  measured  as  a  function  of  increasing  source 
voltage.  A  plot  of  the  harmonic  amplitude  A2  vs.  the  square  of  the  fundamental  a  aptitude 
Aj  yields  a  linear  relationship  over  the  whole  range  of  the  driving  voltages  used  shown 
in  Fig.2  .  The  slope  of  a  linear  fit  through  these  data  is  used  to  calculate  the  no  a  linearity 
parameter.  In  many  of  the  engineering  materials  the  attenuation  of  ultrasonic  waves 
depends  on  their  frequencies.  Unless  the  attenuation  of  the  second  harmonic  is  twice  that  of 
the  fundamental,  a  correction  has  to  be  made  according  to  relationship  [9] 


P  = 
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“me  as. 
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where  oq  and  aj  are  the  attenuation  coefficients  of  the  fundamental  and  the  second 
harmonic  waves  respectively.  The  attenuation  coefficients  for  10  and  20  MHz  longitudinal 
waves  in  the  MMC-specimens  have  been  used  to  perform  this  correction.  Another 
correction  due  to  diffraction  effects  has  been  neglected  because  of  the  large  lateral 
dimensions  of  the  specimens. 


Fig.  1.  Block  diagram  of  the  capacitive  detector  system. 


Fig.  2.  Measured  displacement  amplitudes  for  Fig.  3.  Velocity  designations  and 
fundamental  and  second  harmonic  signals.  direction  of  applied  stress 

Acoustoelastic  constants  are  determined  by  applying  a  uniaxial  compressive  stress 
while  ultrasonic  waves  are  propagating  transversely  to  the  load  direction.  The  relative 
changes  in  velocity  for  three  different  wave  modes  (Fig.3)  are  measured  using  the  pulse 
echo  overlap  technique.  This  technique  is  described  in  detail  in  [10]  and  can  resolve 
velocity  changes  up  to  one  pan  in  106.  The  changes  caused  by  the  acoustoelastic  effect  are 
typically  in  the  order  of  one  part  in  104  and  depend  on  the  polarization  of  the  wave.  The 
slope  of  the  velocity  data  vs.  the  elastic  strain  is  used  to  determine  the  acoustoelastic 
constants  and  then,  with  eqs.(6),  to  calculate  the  third  order  elastic  constants.  The  average 
values  resulting  from  measurements  in  two  different  propagation  directions  within  the 
quasi-isotropic  plane  are  used  to  represent  the  material  parameters. 


RESULTS  AND  DISCUSSION 

From  measurements  of  densities  and  of  ultrasonic  velocities  the  second  order  elastic 
constants  of  the  metal  matrix  composites  have  been  evaluated  by  Lee  et  al.  [13]  using 
eq.(5).  Their  results  are  listed  in  Tab. 3  .  From  this  table  it  is  seen  that  both,  the  Young’s 
and  shear  moduli  increase  linearly  with  the  SiC  content.  The  Young's  modulus  follows 
closely  the  isostress  condition  represented  by 

Ecomp  =  [EAI*ESiC]/[ESiC*fAl  +  EAffSiC]  ( 1 2) 

where  f  is  the  volume  fraction  of  the  indicated  phase. 

The  values  of  the  acoustoelastic  constants  of  the  metal-matrix  composites, 
determined  from  the  relative  changes  in  the  ultrasonic  velocity  as  a  function  of  elastic  strain, 
are  the  average  values  of  the  data  obtained  by  Lee  et  al.  [13]  and  are  listed  in  Tab.2  .  The 
AECs  are  reproducible  within  2%.  From  the  table,  one  can  see  that  the  values  of  the  AECs 
of  the  7064  alloy  are  larger  than  those  of  the  8091  alloy.  Also  the  AECs  for  both  sets  of 
MMCs  decrease  as  the  content  of  SiC  increases.  Smaller  absolute  values  of  AECs  indicate  a 
smaller  change  in  the  ultrasonic  velocity  as  a  function  of  elastic  strain  which  means  a 
smaller  deviation  from  the  ideal  Hookean  behavior  of  the  material. 


Table  2  :  Averaged  acoustoelastic  constants  of  the  examined  MMCs. 


Material 


Acousloelastic  Constants 


AEC  22 

A% 

AEC  21 

A% 

AEC  23 

A% 

A1  -  8091 

1.12 

0 

-2.59 

0 

0.87 

0 

+  10%  SiC 

1.08 

-3.7 

-2.17 

-19.4 

0.82 

-5.7 

+  15%  SiC 

0.92 

-17.9 

-1.91 

-35.6 

0.74 

-14.9 

A1  -  7064 

1.44 

0 

-3.10 

0 

1.93 

0 

+  15%  SiC 

1.20 

-16.7 

-2.35 

-24.2 

1.27 

-34.2 

+20%  SiC 

0.88 

-38.9 

-1.85 

-40.3 

0.92 

-52.3 

From  the  acoustoelastic  constants  and  the  second  order  elastic  constants,  the  third 
order  elastic  constants  were  determined  using  eq.(6).  From  an  analysis  of  error  propagation 
the  inaccuracy  is  found  to  be  highest  in  the  Mumaghan  constant  1  and  is  equal  to  15%.  The 
inaccuracies  of  the  constants  m  and  n  are  5%  and  3%  respectively.  In  general,  the  values 
of  the  Mumaghan  constants  m  and  n  are  not  significantly  influenced  by  the  content  of  SiC 
within  the  experimental  error.  The  constant  m,  however,  has  a  tendency  towards  less 
negative  values,  whereas  the  constant  n  is  showing  no  trend  at  all.  The  Mumaghan  constant 
1  changes  considerably  and  even  becomes  positive  when  SiC  is  added.  Compared  to  the 
elastic  moduli  (i  and  E,  the  third  order  elastic  constants  do  not  show  a  clear  relationship 
with  the  content  of  reinforcement,  whereas  those  of  the  elastic  moduli  increase  in  a 
predictable  manner. 

The  Mumaghan  constants  1  and  m  are  used  to  calculate  the  nonlinearity  parameter 
according  to  eq.(8) .  The  values  of  the  calculated  nonlinearity  parameter  are  listed  in  Tab.4 
as  [kale.  Also  included  in  Tab.4  are  the  values  of  the  directly  measured  nonlinearity 
parameter  as  (3meas.  The  experimental  error  is  10%  for  the  calculated  and  15%  for  the 
measured  nonlinearity  parameter.  For  both  Peak,  as  well  as  pmeas.  the  values  decrease 
considerably  with  increasing  volume  percentages  of  SiC. 


As  can  be  seen  from  Fig.4,  the  values  of  the  calculated  nonlinearity  parameter  for 
the  8091  as  well  as  for  the  7064  alloys  change  linearly  as  a  function  of  second  phase 
content.  The  values  of  the  composites  with  the  A1  -  8091  matrix  appear  to  be  smaller  than 
those  of  the  composites  with  the  A1  -  7064  matrix. 


Table  3  :  Second  and  third  order  elastic  constants  of  MMCs. 


Material  SOEC  [GPa]  TOEC  [GPa] 


\ 

P 

E 

1 

m 

n 

A1  -  8091 

42.0 

30.1 

77.7 

-34 

-320 

438 

+  10%  SiC 

42.8 

35.4 

91.0 

34 

-313 

466 

+  15%  SiC 

42.4 

37.6 

95.7 

33 

-288 

454 

A!  -  7064 

54.1 

26.9 

71.4 

-33 

-359 

-515 

+  15%  SiC 

57.4 

35.1 

91.5 

43 

-343 

-516 

+20%  SiC 

54.3 

38.0 

99.2 

24 

-309 

486 

Table  4  :  Calculated  and  directly  determined  nonlinearity  parameter. 


Material 

[ic  ale. 

A% 

fimeas. 

A% 

At  -  8091 

10.2 

0 

10.6 

0 

+  10%  SiC 

7.4 

-27.5 

8.5 

-19.8 

+  15%  SiC 

6.2 

-39.2 

6.9 

-34.9 

A1  -  7064 

10.9 

0 

8.7 

0 

+  15%  SiC 

7.1 

-34.9 

6.6 

-24.1 

+20%  SiC 

6.2 

-43.1 

5.8 

-33.3 

Fig.  4.  Calculated  nonlinearity  parameter  Fig.  5.  Measured  nonlinearity  parameter 
as  a  funtion  of  particle  content  as  a  funtion  of  particle  content 

A  plot  of  Pmcas.  vs.  the  volume  fraction  of  reinforcement  is  displayed  in  Fig. 5  .  The 
values  of  both  composites  are  close  to  a  linear  relationship  between  [Was.  and  the  volume 
fraction  of  SiC.  One  value,  namely  that  of  the  unreinforced  7064  specimen,  is  clearly 
deviating  from  this  behavior.  Interferences  in  the  backwall  echo  sequence  obtained  on  this 
specimen  indicate  a  strong  texture  in  the  extrusion  direction.  Since  die  nonlinearity 
parameter  vanes  significantly  in  different  lattice  directions  of  single  crystalline  materials,  a 
texture  is  likely  to  change  the  value  of  the  nonlinearity  parameter. 


The  lower  values  for  the  nonlinearity  parameter  of  SiC  reinforced  aluminum  alloys 
can  be  understood  from  the  fact,  that  ceramic  materials  with  low  nonlinearity  will  lower  the 
nonlineanty  parameter  of  aluminum  alloys  when  a  composite  is  formed.  This  does  not 
mean  that  a  law  of  mixture  is  applicable  to  model  the  bulk  nonlinearity  of  a  metal-matrix 
composite.  Influences  from  the  interfacial  region  between  matrix  and  reinforcement  are 
expected  to  contribute  to  the  nonlinearity  of  the  composite. 


Because  there  is  no  significant  difference  between  the  calculated  nonlinearity 
parameter  uetermined  in  the  isotropic  plane  and  the  one  directly  measured  along  the 
extrusion  direction,  it  can  be  assumed  that  the  nonlinearity  parameter  depends  primarily  on 
the  overall  content  of  reinforcement.  Influences  from  the  particle  alignment  in  the  extrusion 
direction  could  not  be  detected  within  the  accuracy  of  the  measurements. 


Theoretical  investigations  by  Cantrell  [  1 2)  have  shown  that  the  nonlinearity 
parameter  is  related  to  the  coefficient  of  thermal  expansion  (CTE)  which  is  also  a  nonlinear 


quantity.  In  previous  studies  [11,  14],  the  coefficient  of  thermal  expansion  in  MMCs  has 
been  found  to  decrease  with  increasing  amounts  of  reinforcement.  Measurements  of  the 
CTEs  of  the  composites  investigated  in  this  paper  are  in  progress. 


CONCLUSIONS 

The  results  of  this  study  show  that  the  acoustoelastic  constants  and  the  acoustic 
nonlinearity  parameter  are  influenced  by  the  amount  of  reinforcement  in  metal-matrix 
composites.  Therefore,  they  are  promising  candidates  to  characterize  the  mechanical 
behavior  of  MMCs  nondestructively.  Also,  the  two  quantities  clearly  indicate  a  decreasing 
elastic  nonlinearity  of  the  composite  with  the  increasing  content  of  SiC.  The  nonlinearity 
parameter  changes  linearly  as  a  function  of  second  phase  content 

The  absolute  values  of  the  calculated  as  well  as  of  the  directly  measured  acoustic 
nonlinearity  parameter  are  in  good  agreement  within  the  accuracy  of  the  measurements. 
This  shows  that  both  techniques,  measurements  of  absolute  amplitudes  using  the  capacitive 
gap  receiver  and  measurements  of  the  acoustoelastic  effect,  are  suitable  methods  for  the 
determination  of  the  nonlinearity  parameter. 


The  direct  measurement  of  the  nonlinearity  parameter  using  the  capacitive  gap 
receiver  requires  a  careful  preparation  of  the  specimen  surfaces.  The  measurement  of  the 
acoustoelastic  effect  is  restricted  to  simple  specimen  geometries  since  it  requires  the 
application  of  external  stresses.  The  selection  between  the  two  methods  depends  on  the 
geometry  condition  of  the  sample. 
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Acoustoelastic  characterisation  of  materials 


J.  H.  Cantrell  and  K.  Salama 


The  increasing  demand  for  new,  more  reliable 
materials,  which  are  often  used  in  hostile 
environments,  has  led  to  the  necessity  of 
establishing  equally  reliable,  quantitative 
techniques  for  the  non-destructive  evaluation 
(NDE)  and  characterisation  of  such  materials. 
Non-destructive  evaluation  methods  are 
commonly  used  in  applications  ranging  from 
materials  processing  and  control  to  monitoring 
the  effects  of  environmental  degradation  and 
the  estimation  of  remaining  useful  life  of 
materials.  Although  linear  ultrasonic  methods 
have  long  been  among  the  most  popular  and 
useful  of  NOE  methodologies,  this  review  is 
concerned  with  the  considerable  effort  that  has 
oeen  expended  recently  on  understanding, 
developing,  and  applying  non-linear 
acoustoelastic  techniques.  Although 
applications  to  complex  materials  are  usually 
correlative,  recent  progress  in  the  quantitative 
modelling  of  the  acoustoelastic  properties  of 
multiphase  alloys  in  terms  of  material 
composition  is  also  reviewed.  Considerable 
emphasis  is  placed  on  understanding  the 
relationship  between  non-linear  acoustoelastic 
properties  and  the  fundamental  atomic 
structure  of  simple  materials.  Such 
considerations  lead  to  an  enhanced 
understanding  of  the  effect  of  residual  and 
applied  stresses  on  the  acoustoelastic  and 
thermoelastic  measurements  of  metallic  alloys. 
Similar  considerations  of  magnetic  domain 
structure  provide  an  explanation  of  the  effect 
of  stress  on  the  magnetoeiastic  properties  of 
ferromagnetic  materials.  Implications  of  these 
advances  to  the  non-destructive 
characterisation  and  evaluation  of  materials 
are  discussed.  IMR/223 
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List  of  symbols 


Lagrangian  (material)  coordinates 
Huang  coefficients 
crack  half  length 

components  of  Boussinesq  stress  tensor 

Brugger  elastic  constants 

depletion  constants 

volume  fraction  of  spherical  inclusions 

total  volume  fraction  of  second  phase 

components 

inverse  of  wave  propagation  time 
acoustoelastic  constants 
magnetic  held  strength 
effectre  stress  intensity  factor 
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bulk  modulus  _  co<{ 

bulk  modulus  of  inclusion*’"*— £-****** — ~~  ^ 

bulk  modulus  of  matrix 

Salama-Ling  proportionality  constants 

compressibility 

wave  propagation  distance 

Murnaghan  constants 

propagation  matrix 

components  of  uniaxial  stress  direction 
components  of  unit  propagation 
direction 

magnitude  of  uniaxial  or  hydrostatic 

stress 

entropy 

compliance  coefficients 
time 

thermodynamic  tension 

temperature 

initial  stress  tensor 

components  of  particle  displacement 

vector 

displacement  gradient 
displacement  gradient  referred  to 
deformed  state 

acoustic  radiation  induced  static 
displacement 

internal  energy  per  unit  mass 
components  of  wave  polarisation  vector 
true  sound  velocity 
longitudinal  wave  true  sound  velocity 
shear  wave  true  sound  velocity 
natural  sound  velocity 
natural  sound  velocity  at  zero  stress 
components  of  particle  position  vector 

direction  cosines 

transformation  coefficients 

thermal  strain  tensor 

thermal  expansivity 

acoustic  non-linearity  parameter 

acoustic  non-linearity  parameter  for 

solid  s  tion 

Kronecker  delta 

direction  cosines 

linear  elastic  strain 

magnetoeiastic  strain 

Lagrangian  strains 

dilatation 

wave  propagation  vector 

Lame  constants 

thermal  stress  tensor 

saturation  magnetostriction  along  (abc: 

shear  modulus 

shear  modulus  of  inclusion 

shear  modulus  of  matrix 

Poisson's  ratio  of  matrix 

amplitudes  of  fundamental  and  second 

harmonic  waves 
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p  average  mass  uensuy  ot  matrix  anu 
inclusion 

p,  mass  density  of  inclusion 
pm  mass  density  of  matrix 
Po  mass  density  of  material  in  unperturbed 
state 

a;,  stress  tensor 
at  angular  frequency 

Introduction 

Ultrasonic  methods  have  long  been  used  to  charac¬ 
terise  various  properties  of  materials  ranging  from 
the  fundamental  to  the  most  practical  applications. 
Ultrasonic  measurements  of  fundamental  material 
properties,  for  example,  have  provided  some  of  the 
first  experimental  evidence  of  the  double  well 
interatomic  potential  in  vitreous  silica.1  They  have 
also  been  used  to  quantify  the  superconducting 
state  of  materials,2  and  have  provided  well  estab¬ 
lished  methods  of  testing  lattice  dynamical  theories 
of  the  solid  state.3-5  Several  recent  review  articles 
and  monographs.6"0  highlighting  the  practical 
applications  of  ultrasonic  methods,  are  testimonials 
to  the  flexibility  and  utility  of  ultrasound  in  assess¬ 
ing  properties  of  materials  from  simple  crack 
detection  to  the  characterisatiorf  of  material  micro¬ 
structure  and  defect  induced  changes  in  elastic 
properties.  Indeed,  ultrasonic  techniques  have 
become  an  indispensable  tool  in  the  non-destructive 
characterisation,  evaluation,  and  testing  of 
materials. 

The  purpose  of  the  present  review  is  to  sum¬ 
marise  recent  developments  in  non-linear  acousto¬ 
elasticity  that  bear  directly  on  the  thermoelastic 
properties,  arid  in  appropriate  cases  on  magneto¬ 
elastic  properties,  of  multiphase  materials.  Non¬ 
linear  acoustoelasticity  entails  consideration  of  the 
variation  in  the  sound  velocity  as  a  function  of  the 
state  of  stress  in  the  material  as  well  as  the 
consideration  of  acoustic  harmonic  generation  and 
radiation  induced  static  stresses  and  strains  in  the 
material  generated  from  an  initially  pure  sinusoidal 
source.  In  order  to  emphasise  the  common  basis  of 
these  related,  but  quite  distinct,  non-linear  phe¬ 
nomena  the  theoretical  treatment  of  each  phenom¬ 
enon  is  initiated  from  the  same  set  of  generalised 
wave  equations.  While  it  will  become  apparent  that 
acoustoelastic  measurements  can  be  used  to 
determine  the  effective  second  and  higher  order 
elastic  constants  as  well  as  various  mechanical 
moduli  of  materials,  the  emphasis  here  is  on  the 
parameters  associated  with  and  defined  from  the 
measurement  process  itself.  These  measurement 
defined  acoustoelastic  constants  and  non-linearitv 
parameters  offer  the  most  direct  route  to  the 
characterisation  of  many  important  material  prop¬ 
erties  from  interatomic  potentials  to  stress  fields 
and  the  state  of  fatigue.  The  connection  with 
fundamental  dynamical  properties  is  emphasised 
throughout  this  review  in  an  effort  to  underscore 
the  role  of  material  anharmonicity  in  acoustoelas¬ 
ticity  3nd  thus  to  provide  some  link  between 
macroscopic  acoustic  measurements  and  the  atomic 
properties  of  elemental  or  single  phase  material. 
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i  nc  complexities  mtrouuceu  in  measurements  ui 
multiphase  materials  are  addressed  from  a  consider¬ 
ation  of  models  giving  the  effective  parameters  as  a 
function  of  total  volume  fraction  of  second  phase 
constituents. 

A  considerable  portion  of  the  review  is  also 
devoted  to  the  effects  of  residual  and  applied 
stresses  on  ultrasonic  measurements.  Indeed,  one 
of  the  most  important  problems  in  the  non¬ 
destructive  characterisation  of  materials  is  the 
determination  of  residual  and  applied  stresses.  Such 
stresses  play  a  fundamental  role,  for  example,  in 
the  corrosion  or  embrittlement  of  metals  and  in 
crack  initiation.  Acoustoelastic  techniques  as  a 
means  of  characterising  such  fields  have  been 
extensively  reviewed  elsewhere,6  and  so  will  not  be 
repeated  here.  Rather,  in  keeping  with  the  spirit  of 
this  review,  the  role  of  stress  fields  as  they  affect 
the  thermoelastic  and,  in  ferromagnetic  solids, 
magnetoelastic  properties  of  materials  is  addressed. 
In  doing  so  recent  advances  in  stress  field  assess¬ 
ment  not  previously  considered  are  presented.  In 
order  to  clarify  the  meaning  of  the  various  elastic 
and  mechanical  moduli  from  a  measurement  per¬ 
spective.  the  connection  between  static  or  low 
frequency  dynamic  moduli  measurements  generally 
performed  in  mechanical  testing  and  the  second 
order  elastic  constants  obtained  in  bulk  ultrasonic 
measurements  is  given  in  Appendix  1. 

The  microstructural  variation  in  acoustic  prop¬ 
erties  has  been  the  key  to  the  development  in  the 
past  decade  of  a  number  of  new  microscopes  which 
utilise  acoustic  waves  in  some  fashion  for  obtaining 
the  microscopic  images."  12  The  images  so  attained 
provide  information  quite  different  from  that 
obtained  in  conventional  light  or  electron  micro¬ 
scopes.  Different  types  of  acoustic  microscopes  use 
different  energy  sources,  such  as  electron  beams, 
lasers,  and  piezoelectric  transducers,  for  generating 
the  acoustic  images  and  are  descriptively  named, 
for  example,  the  scanning  electron  acoustic  micro¬ 
scope.  the  thermoacoustic  microscope,  and  the 
scanning  acoustic  microscope.  Although  the  image 
contrast  and  spatial  resolution  attainable  in  such 
microscopes  are  in  large  measure  dictated  by  the 
details  of  energy  source-material  interactions,  the 
dependency  on  the  local  variations  in  the  material 
elastic  or  thermoelastic  properties,  including 
residual  stresses,  is  a  necessary  consideration  in 
such  interactions.  Indeed,  in  some  cases  features 
directly  associated  with  the  residual  stress  fields  can 
be  imaged.1-’  revealing  even  the  subtlest  influence 
of  thermoelasticity  on  the  contrast  mechanisms. 
The  concepts  summarised  in  this  review,  when 
combined  with  an  appropriate  model  of  the  micro¬ 
scope  under  consideration,  are  pertinent  to  an 
understanding  of  many  of  the  contrast  mechanisms 
in  the  various  acoustic  microscopes,  particularly 
those  involving  acoustic  bulk  wave  propagation. 

General  equations  of  elastic 
wave  motion 

Consider  first  the  adiabatic  propagation  of  an 
elastic,  bulk  travelling  wave  in  a  lossless  solid  of 


arbitrary  crystalline  symmetry.  The  equations  of 
motion  (Newton's  law)  in  Lagrangian  (material) 
coordinates  a ,  are14  (Einstein  summation  of 
repeated  indices  assumed  throughout  text) 


Pu  — 
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00(1 
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where  x,  are  the  components  of  the  particle  position 
vector.  Oi,  are  components  of  the  stress  tensor,  r  is 
time,  and  p0  is  the  mass  density  of  the  solid  in  the 
unperturbed  (natural)  state.  If  the  stress  tensor  is 
defined  in  terms  of  Lagrangian  strain  derivatives  of 
the  internal  energy  per  unit  mass  f/(uk,r;ks,S)  where 
rjks  are  the  Lagrangian  strains  and  S  is  the  entropy, 
the  stress  tensor  is  commonly  called  the  first 
Piola-Kirchhoff  tensor.  Formally,  the  first  Piola- 
Kirchhoff  tensor  is  defined  bv15 


3U 

=  Poar-.k  — —  (2) 

where  cyik  are  the  transformation  coefficients 
defined  by 

on, 

-  d,k  A  —  (3) 

ouk 


In  equation  (3),  (3lk  are  Kronecker  deltas  and  ti-,  = 
a,  -  a,  are  components  of  the  displacement  vector. 
The  Lagrangian  strains  r/k)  are  defined  by1" 


>h\-  =  -  <3jk)/2 


(4) 


Many  researchers  (Brillouin17  and  Born  and 
Huangis  among  others)  have  preferred  to  use  the 
Boussinesq  stress  tensor  B„.  which  is  defined  in 
terms  of  derivatives  of  the  internal  energy  per  unit 
mass  with  respect  to  the  displacement  gradient, 
( SitJSa j)  =  it,,.  Using  equations  (2)-(4).  it  is 
straightforward  to  show  the  equivalence  of  the  two 
stress  tensors  by  writing 


B„  =  Pn 
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Another  stress  tensor  of  relevance  here  is  the 
thermodynamic  tensions  r,,  defined  by 

SU 

U  =  Pn  — —  . (6) 

oth, 

Huang19  expanded  the  internal  energy  per  unit 
mass  in  terms  of  the  displacement  gradients  as 

Put  —  A. ,n,,  —  -  A,,k|U,,uki 

|  -4 1  ,k!mnUijUkiUn,n  ^ . (0 

where  A,,...  are  referred  to  as  the  Huang  (or 
propagation)  coefficients.  Equation  (7)  serves  to 
define  the  nth  order  Huang  coefficients  as 

i  5nU  \ 

=  Pn{- — ; - I  ...  (8) 
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1  he  relationship  between  the  Huang  coefficients 


and  the  elastic  constants  referred  to  the  Lagrangian 
strain  measure  Cijki...  was  found  by  Wallace14  to  be 

A,j  =  C,j  =  7), 
r^ijkl  =  T'jttS.k  A  C„kl 

Hijklmn  C,|mn(5ik  A  U)kn|(5im  A  U,jn|(5km  A  Oijkimn 
. (9) 

The  C.jki ...  are  the  nth  order  elastic  coefficients  of 
Brugger20  defined  for  adiabatic  conditions  by  the 
nth  order  derivative  of  the  internal  energy  per  unit 
mass  with  respect  to  the  Lagrangian  strains 

C,jki...  =  Po{~j- — )  (10) 

Mr/ijdrjki.. 'p=o 

It  is  clear  from  equations  (6)  and  (10)  that  C, j  =  T,t 
is  the  thermodynamic  tension  evaluated  at  T]  =  0. 
From  equations  (5)  and  (3)  we  find  that  a,,  ~  r(j  = 
Cjj  =  T, j  when  oulda,  hence  p,  is  zero.  The 
coefficient  C,j  is  thus  identified  as  the  initial  stress 
7*jj  in  the  solid. 


Acoustoelasticity 

Traditionally,  acoustoelasticity  has  been  viewed  as 
that  aspect  of  material  anharmonicity  giving  rise  to 
the  variation  in  sound  velocity  as  a  function  of  the 
state  of  stress  in  the  material.  Although  acousto- 
elastic  measurements  are  a  popular  means  of  assess¬ 
ing  stress  fields,  the  emphasis  here  is  to  explore  the 
relationships  between  acoustoelasticity  and  certain 
aspects  of  material  microstructure.  In  particular, 
the  use  of  a  measurement  defined  acoustoelastic 
constant  as  a  material  characterisation  parameter  is 
considered.  Applications  of  the  concepts  developed 
here  to  two  phase  materials  are  also  emphasised. 


General  theory 


It  is  necessary  to  obtain  the  elastic  wave  velocities 
as  a  function  of  applied  or  residual  stress.  It  is 
expedient  then  to  expand  the  stress  ctj,  in  equation 
(1)  about  the  homogeneously  deformed  initial  state 
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where  uk  =  »k  —  uk.  the  overbar  represents  the 
value  of  the  parameter  referred  to  the  deformed 
state,  and  the  expansion  coefficients  are  evaluated 
at  the  deformed  state  X.  From  equations  (1),  (5), 
(7).  rnd  (II)  we  thus  obtain 
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Assuming  wave  propagation  of  the  form 

<<v  =  U\  exp  i(wt  -  k  ■  5)  ....  (14) 
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where  /V  =  k!\k\  is  the  unit  propagation  direction 
and  W  =  oj/\k\  is  designated2  as  the  natural  velo¬ 
city,  since  it  is  the  velocity  referred  to  the  natural 
state.  Substituting  equation  (14)  into  equation  (12), 
the  set  of  linear  equations  for  the  unit  polarisation 
vector  U  is  obtained 

(L.juiV^V,  -  p0W26^)U^  =  0  ....  (15) 

Equation  (15)  is  the  basic  equation  of  acoustoelasti¬ 
city.  It  provides  the  variation  in  the  sound  velocity 
(here  the  natural  velocity)  as  a  function  of  applied 
or  residual  strain  and  wave  polarisation. 

Consider  now  the  rate  of  change  of  the  square  of 
the  natural  velocity  with  respect  to  the  applied 
stress  p  (uniaxial  or  hydrostatic)  evaluated  at  the 
natural  configuration  (i.e.  where  p  =  0). 


lowest  crystalline  symmetry  (triclinic  crystal  system) 
this  number  reduces  to  21  independent  second 
order  elastic  constants  and  56  independent  third 
order  elastic  constants.  For  higher  symmetries  the 
number  of  independent  constants  are  even  further 
reduced.  For  example,  the  independent  elastic 
constants  reduce  to  3  second  order  and  6  third 
order  constants  for  cubic  crystals  of  highest  symme¬ 
try  and  for  isotropic  solids  (the  highest  attainable 
symmetry)  the  numbers  are  2  second  order  and  3 
third  order  elastic  constants. 

For  isotropic  solids  the  second  order  constants 
often  used  in  the  literature  are  the  Lame  constants 
A  and  p.  and  the  third  order  constants  are  often  the 
Murnaghan  constants  /,  m ,  and  n.  These  constants 
are  related  to  the  Brugger  elastic  constants  Cijki.. 
for  isotropic  solids  ass 

A  -r  2p.  =  Cu  =  C22  =  C33 

A  C|2  —  C13  -  C23  —  C21  =  C31  =  C32 

/JL  —  C44  =  C55  =  C(£ 

21  +  4 m  ~  Cm 
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21  —  2m  +  n  —  C123 
m  —  nl 2  =  Cu4 
/  =  Cll2/2 


where  Srs3b  =  (<5r/rs/5f3b)p=o  are  the  isothermal 
compliance  coefficients.  For  applied  hydrostatic 
pressure  or  uniaxial  compression  of  magnitude  p 
in  the  direction  \t 


£4b) 

op  p=0 


-(5jb  for  hydrostatic  pressure 

-A/3A/b  for  uniaxial 
compression 


(17) 


Substituting  equations  (9),  (13).  and  (15)  (for  the 
natural  state)  into  equation  (16)  gives  the  equations 
of  Thurston  and  Brugger1 
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where 
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Equations  (18)  and  (19),  derived  for  solids  of 
arbitrary  crystalline  symmetry,  specifically  show  the 
dependence  of  the  natural  velocity  and  the  change 
of  natural  velocity  with  stress  in  terms  of  the  second 
and  third  order  elastic  constants  of  the  solid.  From 
the  experimental  measurements  of  natural  velocity 
as  a  function  of  stress  the  equations  suggest  that 
one  can  determine  the  second  and  third  order 
elastic  constants  of  the  material.  Indeed,  calculation 
of  the  elastic  constants  from  such  measurements 
have  become  quite  well  established.3^-21  There  are. 
in  general.  81  second  order  elastic  constants  and 
729  third  order  elastic  constants.  For  solids  of 


m  =  Ci  55 

n  =  4C456  . (20) 

w-here  all  other  second  and  third  order  elastic 
constants  are  zero.  Voigt  notation  is  used  in  the 
above  equations  for  the  Brugger  elastic  constants 
whereby  the  single  subscript  k  replaces  the  sub¬ 
scripted  pair  (/,/)  according  to  the  following 
scheme: 

/,  j  =  11  22  33  32  or  23  31  or  13  21  or  12 

k  =  1  2  3  4  5  6 

In  a  typical  ultrasonic  experiment,  one  generally 
measures  as  a  function  of  applied  stress  the 
variation  in  time  required  for  a  generated  plane 
w-ave  front  to  propagate  between  parallel  sample 
surfaces.  An  inverse  measure  of  that  propagation 
time  is  the  parameter  F  defined  by 

F=l~lv . (21) 

where  v  is  the  true  sound  velocity  and  /  is  the 
propagation  distance  in  the  sample  in  the  homo¬ 
geneously  deformed  state.  Thurston  and  Brugger21 
have  shown  that  the  true  velocity  is  related  to  the 
velocity  referred  to  the  natural  state  W  by 

v  =  /u'/W  . (22) 

where  /«  is  the  propagation  distance  in  the  undefor- 
med  sample.  From  equations  (18)  and  (22)  and  an 
expansion  of  (/q1/)  in  terms  of  the  Lagrangian 
strains, one  can  obtain  the  expressions  of  Hughes 
and  Kelly23  derived  from  isotropic  solids  which 
relate  the  true  sound  velocity  to  the  applied  stress 
a.  Their  expressions,  for  example,  for  axial  stress 
applied  in  a  direction  perpendicular  to  the  wave 
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propagation  direction  are 

,  ,  _  a  [  2A 

p i iV'7  —  A  +  2 p. - 21  —  —  (m  +  A  4-  2/a) 

3A'0  L  p 

,  a  [  hi 

W2  -  M  -  - —  /n  +  —  4-  /.  -f  2/x 

3A'o  c  4/i 

(7  A  +  /A  I 

p,A'j  =  M  -  — —  m - /i  -  2A  I  .  .  (23) 

->A:0  c  2/a  J 

where  v/  is  the  longitudinal  velocity  and  v:  and  v3 
are  the  velocities  of  shear  waves  when  the  polari¬ 
sation  direction  is  parallel  to  and  perpendicular  to 
that  of  stress,  respectively. 

There  is  advantage,  however,  in  measuring 
changes  in  the  natural  velocity  as  a  function  of 
applied  stress,  since  such  measurements  are 
obtained  trom  time  measurements  in  the  deformed 
state  but  the  length  measurements  are  referred  to 
the  unaeformed  length  of  the  solid.  This  is  seen  bv 
substituting  equation  (21)  into  equation  (22)  to  get 

F=tfW  . (24) 

From  equation  (24) 


_1 _ o(pVV2) 

F  5p  \V  Sp  2pu\V2  Sp 


(25) 


When  evaluated  at  p  =  0,  it  has  become  popular  to 
refer  to  H  Uj  the  'stress  acoustic  constant'  or  the 
'acoustoelastic’  constant.  It  is  clearly  a  material 
parameter,  since  it  depends  on  the  second  and  third 
order  e last  1  _•  constants  of  the  material.  Cantrell" 
has  determined  the  general  relationship  between 
the  fractional  change  in  the  natural  velocity  with 
respect  to  stress  for  a  wave  of  polarisation  0 
propagating  along  ;V,  i.e.  the  generalised  acousto¬ 
elastic  constant  /7UV(  t/,:V),  and  the  fractional 
change  in  the  true  velocitv  v  to  be 


=  f-v—  -  A',.V$)sr„v  .  .  (26) 

where  prs  are  the  Lagrangian  strains  resulting  from 
the  imposed  stress,  and  5rsuv  are  the  compliance 
coetficients.  For  acoustic  wave  propagation  perpen¬ 
dicular  to  the  direction  of  applied  stress  in  isotropic 
solids,  equation  (26)  shows  that  acoustoelastic 
constants  are  also  equal  to  the  fractional  change  in 
the  true  velocity  with  respect  to  the  applied  stress. 
For  acoustic  longitudinal  wave  propagation  along 
the  direction  of  applied  stress  in  such  solids,  the 
fractional  variations  in  the  true  and  natural  veloci¬ 
ties  differ  by  an  additive  constant. 

The  exact  meaning  of  the  experimental 
parameter  F  in  equation  (25)  depends  on  the 
particular  experimental  technique.  For  example,  if 
one  uses  a  resonance  or  resonance  derived  tech¬ 
nique.  then  F  is  the  acoustic  standing-wave  reson¬ 
ance  frequency  f.  If  one  uses  a  pulse  coincidence 
technique,  then  F  is  the  inverse  pulse  repetition 
rate.  It  is  important  to  emphasise  that  all  the  above 
equations  specifically  refer  to  acoustic  bulk  wave 
propagation  in  materials.  The  elastic  moduli  or 
constants  used  to  quantify  this  wave  propagation 


are  constants  referred  to  unbounded  or  bulk  sotios 
and  accurately  represent  acoustic  wave  propagation 
when  the  dimensions  of  the  sample  are  generally 
large  compared  with  an  acoustic  wavelength.  Such 
elastic  constants,  however,  are  quantitatively 
related  to  the  elastic  or  mechanical  moduli  typically 
obtained  from  static  engineering  measurements  or 
low  frequency  dynamic  measurements  as  shown  in 
Appendix  1. 

In  general,  for  a  given  comoination  of  stress, 
propagation,  and  polarisation  directions,  the 
relationship  between  the  natural  sound  velocity  and 
stress  may  be  obtained  by  integrating  equation  (25) 
to  get 

VV(o)  =  +  W<yHo  . (27) 

where  W0  -  14(0),  the  natural  velocity  at  rero 
stress.  As  suggested  by  equation  (27)  unknown 
stresses  can  be  determined  when  both  the  natural 
velocity  in  the  absence  of  stress  and  the  acoustvj- 
elastic  constant  are  known  independently.  This 
approach  was  first  used  by  Crecraft,24  who  showed 
that  acoustoelasticity  could  be  used  as  a  practical 
method  for  determining  applied  and  residual  stres¬ 
ses  in  engineering  materials.  Theoretically,  one 
could  calculate  the  acoustoelastic  constants  from 
equations  (18),  (19),  and  (25)  if  the  second  and 
third  order  elastic  constants  are  known.  Unfortu¬ 
nately,  measurements  of  the  third  order  elastic 
constants  have  been  performed  mostly  on  pure 
materials  and  some  simple  alloys.  Indeed,  the 
technique  mostly  used  for  such  measurements  have 
been  that  of  determining  the  stress  dependence  of 
the  natural  velocity,  i.e.  the  acoustoelastic  constant 
of  the  materia).  It  is  important  to  point  out  here 
that  the  measured  velocity  and  stress  dependent 
velocity  changes  in  engineering  materials  depend 
strongly  on  microstructural  features  which  make  it 
necessary  to  perform  a  calibration  between  velocity 
and  stress  (i.e.  determine  the  acoustoelastic  con¬ 
stant)  for  each  material  in  order  for  the  method  to 
be  used  in  the  determination  of  unknown  stresses. 
Microhomogeneity,  texture,  and  weak  anisotropy 
of  the  material,  which  are  usually  neglected  in 
engineering  applications  of  the  theory  of  elasticity, 
cannot  be  neglected  in  the  applications  of  acousto¬ 
elasticity.  Results  by  Smith  et  al.2i  show  that  third 
order  elastic  constants  for  polycrvstalline  materials 
can  differ  widely,  even  for  alloys  having  the  same 
composition.  For  structural  aluminium  of  slightly 
different  composition  differences  as  large  as  50% 
are  observed. 

Models  of  two  phase  systems 

Most  experimental  measurements  of  the  acousto¬ 
elastic  constants  have  been  performed  either  on 
simple,  single  phase  materials  or  on  complex, 
multiphase  materials  that  have  been  treated  as 
single  phase  materials.  It  is  important  to  recognise 
that  in  quantitatively  characterising  complex 
materials  one  must  also  understand  the  effects  of 
material  composition  on  the  measured  acoustic 
parameters.  In  order  to  calculate  the  acoustoeiastic 
constant  in  a  two  phase  alloy,  Salama  et  alF6 
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assume  that  the  precipitates  in  the  alloy  can  be 
represented  by  a  dilute  suspension  of  spherical 
particles  in  an  infinite  solid  solution  matrix.  The 
physical  meaning  of  this  idealisation  is  that  the 
precipitates  are  so  small  and  so  far  apart  that  all  the 
interactions  between  them  can  be  neglected,  no 
matter  what  the  size  of  the  representative  volume 
may  be.  The  effective  elastic  constants  of  such 
materials  have  been  calculated  by  Christensen2- 
using  the  energy  methods  approach  and  Eshelby's 
formula28  for  the  calculation  of  strain  energy  in 
systems  containing  inhomogeneities.  For  a  dilute 
concentration  of  inclusion  i,  Christensen  found  that 
the  effective  shear  and  bulk  moduli  p  and  k  are 
related  to  those  of  the  matrix  as 

M  ;  15(1  -  vm)[l  -  (pjp<n)\f 

Mm  (7  ~  5vm)  4-  2(4  -  5vm)(pJ  pm) 

and 


[1  ~  (fcj/fcm))/ 

1  -r  [(k,  -  km)!(km  -r  4  pj 3) 


(29) 


where  the  c,ibscripts  tn  and  /,  respectively,  denote 
matrix  and  inclusion,  /  is  the  volume  fraction  of  the 
spherical  inclusions  under  dilute  conditions,  and  vm 
is  Poisson's  ratio  of  the  matrix.  The  ultrasonic 
longitudinal  velocity  vt  is  re lated  to  the  elastic 
moduli  p  and  k  as'J 


Pl'T 


k  ~  4/x/3 


(30) 


where  p  is  the  average  mass  density  of  the  medium 
(matrix  plus  inclusions). 


Case  1 

For  an  inclusion  with  elastic  moduli  not  too 
different  from  those  of  the  matrix,  such  that 


M,  ~  Mm  and  k,  ~  km.  but  p,  -  Mm  -  0  and 
k,  -  km  --  0,  equations  (28)  and  (29),  respectively, 
become 

Mi  =  Mm  1-  (Mi  -  Mm)/  . (31) 

and 

k  ~  km  4-  (k,  -  km)f . (32) 


Equations  (31)  and  (32)  indicate  that  the  elastic 
moduli  of  a  two  phase  material  is  a  function  of  the 
inclusion  concentration  as  well  as  the  difference 
between  the  corresponding  constants  of  the  inclu¬ 
sion  and  matrix.  Substituting  equations  (31)  and 
(32)  into  equation  (20)  gives  for  the  longitudinal 
velocities 

pV]  Pm(C)rn  “^[pm(^l)ni  “  Pi(^l)i]  fOj) 


For  alloys  containing  dilute  concentrations  of  pre¬ 
cipitates,  p  —  pm  and  equation  (33)  can  be  written 

as 


v7  —  (vT)m 


( t'T  )  m 


(34) 


Fiquation  (34)  indicates  that  the  longitudinal  velo¬ 
city  in  the  two  phase  alloy  increases  or  decreases 
from  that  of  the  solid  solution  phase  according  to 
whether  tiie  density  of  the  precipitates  is 
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respectively  higher  or  lower  than  that  of  the  solid 
solution  phase,  since  the  difference  between  the 
longitudinal  velocities  in  the  precipitates  and  the 
solid  solution  is  usually  small. 

Consider  the  application  of  a  uniaxial  stress  in  a 
direction  perpendicular  to  that  of  the  longitudinal 
wave  propagation  direction.  Differentiating  equa¬ 
tion  (34)  with  respect  to  stress  a  and  assuming  that 
in  the  dilute  approximation  (i’i),,,  =  v,  and  that, 
tenuously,  the  densities  of  the  two  phases  present 
remain  roughly  unchanged  with  stress  variations, 
the  fractional  change  in  ultrasonic  velocity  in  the 
two  phase  material  is  estimated  as 

2  (22  =  1  (21s) 

IV)  ''do '  v’i  '  da ' 


From  the  definition  of  the  acoustoelastic  constant, 
equation  (25),  it  is  concluded  from  equation  (35) 
that  the  change  in  the  effective  longitudinal  wave 
acoustoelastic  constant  for  the  two  phase  system 
caused  by  the  presence  of  precipitates  is  linearly 
proportional  to  the  volume  concentration  of  the 
inclusion  and  depends  on  the  elastic  moduli  of  the 
precipitates  and  the  solid  solution  phases,  as  well  as 
their  stress  dependences.  These  latter  quantities  are 
generally  unavailable  and  make  the  calculations 
of  the  effective  acoustoelastic  constants  rather 
difficult. 


Case  2 

When  the  elastic  moduli  of  the  precipitate  are  much 
greater  than  those  of  the  matrix,  i.e.  p,  »  pm 
and  k,  »  km.  equations  (2S)  and  (29)  become, 
respectively 

15(*  "  Vm)/ 

M  Mm  Mm  _  . (-16) 

2(4  -  Dv,n) 


and 


k  =  k„  +  {km  +  4pJ3)  f  . (37) 

These  equations  indicate  that  the  effective  sheai 
and  bulk  moduli  of  the  two  phase  alloy  system  do 
not  depend  on  those  of  the  precipitates:  however, 
they  van/  linearly  with  their  concentrations.  Substi¬ 
tuting  equa.ions  (36)  and  (37)  into  equation  (30) 
and  considering  v*  =  i’:(2(i  —  v))1'2,  where  is  the 
shear  wave  velocity,  we  obtain  in  the  dilute 
approximation  that 


( 1  -  f)v\  -  ( v  1  )  m  + 


3(3  -  5vm) 
(4  -  5vm) 


f(v f)m 


(38) 


Again  in  order  to  obtain  an  expression  for  the 
effective  acoustoeiastic  constant  in  alloys  containing 
precipitates  with  elastic  moduli  much  higher  than 
those  of  the  solid  solution,  equation  (38)  is  differen¬ 
tiated  with  respect  to  the  stress  a.  again  under  the 
condition  that  the  stress  is  applied  in  a  direction 
perpendicular  to  the  wave  propagation  direction. 


Under  such  conditions  and  considering  the 
relationship  between  the  longitudinal  and  shear 
wave  velocities,  17=  vy[2(l  —  v)/(  1  -  2v)|,  the 
relative  change  in  the  effective  longitudinal  wave 
acoustoelastic  constants  H\  is  given  as 

■  (-'9) 
where 


15(1  -  vj(l  -  2vm) 
(4  -  5vm)2 


and  where  (f/;)m  is  the  shear  wave  acoustoeiastic 
constant  and  (H ,)m  is  the  longitudinal  wave  acous¬ 
toeiastic  constant  of  the  matrix  material.  Equation 
(39)  indicates  that  the  relative  change  in  the 
effective  longitudinal  acoustoeiastic  constant  H,  of  a 
two  phase  material  is  linearly  proportional  to  the 
volume  friction  of  the  dilute  second  phase.  The 
constant  of  proportionality  depends  on  the  Pois¬ 
son's  ratio  of  the  solid  solution  phase  and  the 
relative  difference  between  longitudinal  and  shear 
wave  acoustoeiastic  constants.  Comparison  of  the 
theoretical  models  with  experimental  remits  is 
presented  in  the  section  on  'Experimental  results’ 
below,  but  first  techniques  for  measuring  the 
acoustoeiastic  constants  are  briefly  discussed. 


Measurement  techniques 

The  ultrasonic  measurement  techniques  generally 
used  in  the  determination  of  acoustoeiastic  con¬ 
stants  arc  based  the  propagation  of  continuous 
waves  (C'.V),  radiotrequency  (rf)  bursts  (tonebur- 
sts),  or  broadband  pulses  in  the  medium  of  interest. 
Comprehensive  descriptions  of  these  techniques 
have  been  given  in  several  review  articles  and 
monographs,  e.g.  McSkimin.29  Truell  et  a!..A)  Papa- 
dakisT'1  Breazeale  et  al..  ':  and  Ratcliff.-’-'  In  view  of 
the  extensive  review  iiterature  on  acoustoeiastic 
measurement  methodology  details  are  not  given 
here,  but  it  is  pointed  out  that  the  techniques 
generally  are  based  on  measured  parameters  that 
are  ultimately  referred  to  a  time  standard.  The 
specific  technique  may  utilise  a  standing-wave  pat¬ 
tern.  a  pulse  coincidence  or  pulse  overlap  scheme  in 


which  a  certain  phase  relationship  among  the  pulses 
is  maintained,  or  even  a  hybrid  of  the  above.  In  any 
case  the  measurement  parameter  of  interest  is 
generally  a  frequency  or  the  inverse  of  a  pulse 
repetition  rate  which  is  measured  as  a  func'ion  of 
the  impressed  strain  or  stress  in  the  material. 

Experimental  results 

In  order  to  study  the  sensitivity  of  acoustoeiastic 
constants  to  changes  in  microstructure,  Heyman  et 
al .-U  measured  the  acoustoeiastic  constants  in  the 
four  carbon  steels  AISI  1020,  A1SI  1045,  AISI 
1095,  and  ASTM  533B.  In  steel  alloys,  carbon  in 
excess  of  the  solubility  limit  (0  02  wt-%)  forms  a 
second  phase,  namely  iron  carbide  (cementite), 
which  precipitates  from  the  solid  or-ferrite. 
Compared  with  ferr,>e,  cementite  is  very  hard  and 
its  presence  within  the  turnte  increases  the  strength 
of  the  steel.  Table  1  lists  the  values  of  the 
acoustoeiastic  constants  and  the  longitudinal  wave 
velocities  obtained  for  these  steels.  Also  included  in 
the  table  are  the  volume  fraction  of  second  phase 
precipitates  in  these  alloys  and  the  percentage 
changes  in  the  acoustoeiastic  constants  with  respect 
to  those  of  the  100%  solid  solution.  The  latter  are 
obtained  bv  extrapolating  the  experimental  values 
of  the  acoustoeiastic  constants  to  100%  solid 
solution.  The  relative  variations  in  the  acoustoeias¬ 
tic  constants  as  a  function  of  volume  fraction  of 
precipitates  are  plotted  in  Fig.  1,  which  shows  that 
the  acoustoeiastic  constant  decreases  linearly  as  the 
amount  of  carbide  (cementite)  phase  is  increased  in 
the  alloy.  In  this  study  the  amount  of  carbide  phase 
is  calculated  using  the  lever  rule  and  the  nominal 
carbon  content  in  the  alloy.  These  results  were  also 
confirmed  by  Allison  et  al.2*  where  the  acoustoeias¬ 
tic  constants  in  the  four  steels  AISI  1016,  AISI 
1045.  AISI  1095,  and  AISI  8620  are  found  to 
decrease  linearly  with  the  increase  of  the  carbide 
phase. 

More  recently  Schneider  et  ai.-'b  studied  the 
relationship  between  the  acoustoeiastic  constants 
and  percentage  of  precipitates  in  the  five  aluminium 
aliovs  1100.  3003.  5052,  6061,  and  2024.  Aluminium 
alloys  that  contain  small  percentages  of  other 
elements  such  as  copper,  magnesium,  and  silicon 


Table  1  Variations  of  acoustoeiastic  constants  and  longitudinal  ultrasonic  veloci¬ 
ties  with  volume  fraction  of  compound  precipitates  in  carbon  steels  and 
heat  treatable  aluminium  alloys 


Alloy 

Volume 
fraction  of 
precipitates.  % 

Velocity, 
m  s_1 

Acoustoeiastic 

constant. 

G?3 

Change  in 
ascoustoelastic 
constant.  % 

Steel 

100%  ferrite 

0 

440* 

0 

AISI  1020 

3  1 

5889 

418 

5-0 

ASTM  5338 

3  9 

$380 

410 

6  8 

AiSI  1C45 

7-3 

5383 

375 

14-8 

AISI  1095 

14-7 

5910 

325 

26-1 

Aluminium 

100ao  solid  solution 

0 

83* 

0 

6061-T5  (Al-Mg-Sij 

3  0 

6112 

79 

4  8 

2CM-T351  tAI-Cu! 

5-9 

6 1 53 

74-7 

10-0 

*  bro-n  extraeoia'ion  of  acoustoeiashc  constants  v.  volume  fraction 
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%  Carbide  Phase 


1  Graph  of  acoustoelastic  constant  of  carbon  steel 
as  function  of  volume  percent  carbide  (cementite) 
phase 

form  intermetallic  precipitates  such  as  CuAF, 
CuMgAU,  and  Mg;Si.  These  precipitates  are  very 
hard  and  cause  the  strengthening  of  the  aluminium 
alloy.  Some  values  of  the  acoustoelastic  constants 
and  the  longitudinal  wave  velocities  are  also  given 
in  Table  1.  The  results  of  this  study,  showm  in  Fig. 
2,  reveal  that  the  acoustoelastic  constants  in  the 
heat  reatable  alloys  6061  and  2024  decrease  as  a 
function  of  the  second  phase  precipitates.  The 
behaviour  in  the  work  turdenable  alloys  1100. 
3003.  and  5052,  however,  is  opposite,  i.e.  the 
acoustoelastic  constants  increase  with  the  increase 
of  percentage  of  precipitates. 

A  plot  of  the  acoustoelastic  constant  as  function 
of  aging  time  in  aluminium  alloys  is  showm  in  Fig.  3. 
The  results  indicate  that  the  acoustoelastic  constant 
does  not  change  significantly  with  aging  time, 
although  the  average  size  of  the  precipitate  particles 
are  expected  to  change.37  It  is  thus  inferred  that  the 
acoustoelastic  constant  is  not  significantly  influ¬ 
enced  by  changes  in  the  size  and  distribution  of 
second  phase  precipitate  particles,  at  least  of  the 
size  represented  in  the  figure. 

Equations  (35)  and  (39)  obtained  above  repre¬ 
sent  me  changes  in  the  acoustoelastic  constants 
when  a  second  phase  precipitate  is  in  equilibrium 
with  a  solid  phase.  By  considering  that  the  changes 
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%  Second  Phase  Precipitate 

2  Graph  of  acoustoelastic  constant  of  heat  treat¬ 
able  Al  alloy  6061  (Al-Mg-Si)  and  2024  (Al-Cu)  as 
function  of  volume  fraction  of  second  phase 
precipitates 
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3  Graph  of  acoustoelastic  constant  of  aluminium 
alloys  as  function  of  aging  time 

in  the  acoustoelastic  constants  of  carbon  steels 
measured  bv  Heyman  et  al.:'~  and  by  Allison  et 
al.J'6  and  of  the  aluminium  alloys  6061-T6  and 
2024-T351  measured  by  Schneider  et  al/6  are  due 
only  to  the  presence  of  precipitates  in  the  solid 
solution  phase,  a  comparison  between  theory  and 
experiment  can  be  made.  Further,  since  the  precipi¬ 
tates  in  both  systems  investigated  are  of  higher 
rigidity  than  those  of  the  solid  solution  phase, 
equation  (39)  is  considered  to  be  a  better  represen¬ 
tation  of  the  experimental  data. 

Figures  1  and  2  indicate  that  the  relative  changes 
in  the  acoustoelastic  constants  for  the  longitudinal 
waves  vary  linearly  with  the  volume  fraction  of 
second  phase,  as  predicted  by  equation  (39).  The 
figures  also  show  that  the  constant  of  proportion¬ 
ality  between  the  relative  changes  of  the  acousto¬ 
elastic  constants  and  the  volume  fraction  of  second 
phase  in  the  steel  and  aluminium  systems  are  the 
same  and  approximately  equal  to  TS.  This  means 
that  the  acoustoeiastic  constants  in  these  alloys  can 
be  calculated  empirically  using  the  values  of  this 
quantity  in  the  100%  solid  solution  material  and  the 
volume  fraction  of  second  phase  in  the  alloy.  This 
acoustoeiastic  constant  is  the  proper  vaiue  to  be 
used  in  the  residual  stress  determination  of  the 
alloy. 

Considerable  research  effort  has  been  devoted  in 
recent  years  to  the  development  of  metal  matrix 
composites.  The  non-destructive  characterisation  of 
the  elastic  and  anelastic  behaviour  of  such  materials 
is  essential  to  their  manufacturing  processes.  The 
effects  on  the  measured  elastic  properties  of 
variations  in  tne  volume  fraction  of  second  phase 
SiC  particles  in  a  number  of  aluminium  alloy 
matrices  have  been  reported. ’s~i0  Of  particular 
interest  here  are  the  measured  values  of  the 
acoustoeiastic  constants  as  a  function  of  the  volume 
percent  SiC  in  the  matrix.  A  typical  data  set  is 
hown  in  Fig.  4  for  SiC  particles  in  a  7064  alu¬ 
minium  alloy  matrix.  The  relationship  between  the 
acoustoeiastic  constant  and  the  percent  SiC  is  seen 
to  be  approximately  linear  up  to  roughly  1C% 
volume  fraction.  This  behaviour  ii  a  nominal 
agreement  with  equation  (39)  and  thus  indicates 


4  Graph  of  acoustoelastic  constant  as  function  of 
volume  percent  SiC  particles  in  7064  (Al-Zn-Mg) 


that  up  to  10%  volume  fraction,  the  SiC  particles 
are  considered  to  he  of  dilute  concentration  in  the 
aluminium  matrix.  Thereafter,  the  curve  rapidly 
becomes  non-linear. 


Non-linear  wave  propagation 

A  further  aspect  of  acoustoelasticity  is  the  distor¬ 
tion  of  an  acoustic  waveform  as  it  propagates 
through  a  material.  Such  distortion  is  another 
manifestation  of  the  same  material  interatomic 
anharmonicity  that  governs  the  stress  dependence 
of  the  sound  velocity,  but  it  results  from  a  rather 
different  mechanism  involving  self-induced  stress 
fields.  The  quantitative  measure  of  such  non-linear 
wave  propagation  is  the  acoustic  non-linearity 
parameter  /3.  It.  like  the  acoustoelastic  constant,  is 
found  to  be  a  very  useful  parameter  for  assessing 
material  properties. 


General  theory 

Consider  the  solid  to  be  in  the  natural  state,  i.e. 
r,,  -  0.  Substituting  equations  (5)  and  (7).  into 
equation  (I)  and  retaining  first  order  non-linear 
terms  gives 

,  ~~  'h  _  |  ,  .  I  -  /11V 

P«  _  .  |  A,, Cl  ■  >iiklrnn  ~Z  „  \  1 1  / 

or  -  i 


The  self-resonant  or  mutual  resonant  solution  to 
equation  (41)  assuming  the  boundary  condition 

u,  =  e Lficoswr  at  a  =  0  ....  (42) 

is  to  a  first  order  in  the  non-linearity14 

u,  -  (//,)  -  c;C,cos(eof  -  k:  ■  a) 

-i-  /jij-A.-/6',cos(2t:tf  —  2k  ■  J)/S  +  ...  (43) 

where  /  is  the  distance  of  wave  propagation  and  />  is 
•fie  non-lineantv  parameter  defined  by41 
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(44) 


Just  as  in  the  case  of  the  acoustoelastic  coefficient, 
the  non-linear. rv  parameter  is  seen  from  equations 


\'J }  ailu  (  i  iw  uCpUia  uu  wiv,  oCwuiiu  auu  timu 

order  elastic  constants  of  the  material,  and  the 
directions  of  wave  propagation  and  polarisation. 

The  term  (it,)  in  equation  (43)  is  a  constant  that 
represents  a  static  or  dc  displacement  arising  from 
the  non-linearity  of  the  material.  The  magnitude  of 
the  static  displacement  has  been  the  source  of  some 
controversy  which  has  been  recently  resolved  by 
Cantrell  ei  al. 4:43  They  showed  that  the  static 
displacement  is  intrinsically  linked  with  the  exist¬ 
ence  of  a  non-zero  acoustic  radiation  stress  in  the 
solid  which  when  properly  included  in  the  pertur¬ 
bation  solution  to  equation  (41)  yields  the  value 

<«.>  =  /3c,V/f/,/S  . (45) 

The  sign  of  the  static  displacement  is  governed  by 
the  sign  of  the  non-linearity  parameter.  For  pure 
mode  acoustic  wave  propagation  directions  the  sign 
of  the  non-linearity  parameter  is  typically  positive 
and  according  to  equation  (45)  leads  to  a  radiation 
induced  static  dilation  in  the  solid  within  the  spatial 
extent  of  the  acoustic  wave.  A  negative  non¬ 
linearity  parameter,  on  the  other  hand,  gives  rise  to 
a  radiation  induced  static  contraction.  This  phe¬ 
nomenon  has  been  confirmed  in  single  crystal 
silicon.4’-44  germanium.4^  and  vitreous  silica. 44 


Measurement  technique 

Equation  (43)  indicates  that  the  acoustic  non¬ 
linearity  parameters  may  be  obtained  either  from 
measurements  of  the  absolute  amplitudes  of  the 
fundamental  acoustic  wave  and  generated  second 
harmonic  wave  or  from  measurements  of  the 
fundamental  and  static  displacement  signals.  Spe¬ 
cifically.  if  gi  and  are  'he  measured  amplitudes  of 
the  fundamental  and  second  harmonic  signals,  then 
the  non-linearity  parameter  /3  is  determined  by 


A  block  diagram  of  the  equipment  arrangement 
typically  used  in  non-linearity  parameter  meas¬ 
urements  is  shown  in  Fig.  5.  The  diagram  shows  the 
setup  for  both  the  harmonic  generation  and  the 
static  displacement  techniques.  A  rf  signal  (typically 
in  the  range  5-30  MHz)  from  a  CW  oscillator  is 
combined  in  a  mixer  with  a  gating  pulse  from  a 
logic  and  timing  generator.  The  rf  pulse  is  amplified 
by  a  broadband  linear  rf  power  amplifier  and  is 
used  to  drive  a  narrow  band  lithium  niobate 
transducer  bonded  to  one  end  of  a  cylindrical 
sample.  The  transducer  generates  a  g3ted  ultrasonic 
signal  or  toneburst  which  propagates  through  the 
sample.  The  toneburst  is  received  by  a  broadband 
air-gap  capacitive  transducer  at  the  opposite  end  of 
the  sample. 

The  capacitive  transducer  is  a  parallel  plate 
arrangement  in  which  the  sample  end  surface 
functions  as  the  ground  plate.  The  other  plate  is  an 
optically  flat  electrode  which  is  recessed  approxi¬ 
mately  7  urn  from  the  sample  surface.  A  dc  bias 
(tvpicallv  50-150  V)  is  applied  through  a  1  MQ 
resistor  to  the  electrode.  Gauster  3nd  Breazeale4* 
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5  Block  diagram  of  typical  equipment  arrangement  for  measurements  of  acoustic  non-linearity  parameters 


have  shown  that  the  measurement  of  the  output 
voltage  from  such  a  capacitive  arrangement  may  be 
converted  to  a  displacement  measurement  of  the 
sample  surface  using  the  relation 

V'j  =  V'b  («/5)  . (47) 

where  Vd  is  the  output  voltage,  V’b  the  dc  bias 
voltage,  «  the  displacement  of  the  sample  surface, 
and  s  the  gap  spacing  between  the  capacitor  plates. 
Yost  and  Breazeale47  have  shown  that  the 
capacitive  transducer  is  capable  of  measuring 
displacement  amplitudes  of  the  order  10~16  m, 
which  is  quite  adequate  for  measurements  m  the 
range  10-lj-l0-,°  in  representing  typical  ampli¬ 
tudes  of  second  harmonic  and  static  displacement 
signals. 

The  displacement  amplitudes  of  the  fundamental 
and  harmonically  generated  toneburst  are  obtained 
by  switching  the  output  of  the  capacitive  transducer 
into  the  appropriate  (preferably  narrow  band  if) 
preamplifier.  The  output  of  the  preamplifier  is  sent 
into  a  rectifier  and  filter  assembly  which  converts 
the  toneburst  into  a  detected  signal.  The  output  of 
the  sample  and  hold  is  measured  by  a  voltmeter 
and  recorded.  The  capacitive  transducer  is  now 
switched  out  and  a  Thevenin  equivalent  network  is 
switched  into  the  circuit.  A  substitutional  cali¬ 
bration  signal  corresponding  to  either  the 
fundamental  or  harmonically  generated  signal,  as 
the  case  may  be.  is  now  switched  into  the  Thevenin 
equivalent  network  and  adjusted  in  amplitude  until 
the  voltmeter  connected  to  the  sample  and  hold 
reads  the  same  value  as  that  of  the  ultrasonic 
toneburst  measurement.  Since  the  amplitude  of  the 
calibration  signal  is  measured  at  the  input  to  the 
Thevenin  equivalent  network,  the  measured  value 
is  equal  to  the  signal  produced  by  the  capacitive 
transducer. 

In  order  to  display  the  acoustic  radiation  induced 
static  displacement  pulse  the  output  of  the 


capacitive  transducer  is  connected  to  a  low  fre¬ 
quency  (of  order  3-300  kHz)  preamplifier  by  a 
short,  low  capacitance  lead.  The  output  of  the 
preamplifier  is  sent  to  an  oscilloscope  where  the 
static  displacement  signal  is  displayed.  A  template 
is  placed  over  the  screen  to  record  the  static 
displacement  shape.  The  shape  of  the  static 
displacement  signal  is  that  of  a  right  angled  triangle 
whose  slope  is  equal  to  the  right-hand  side  of 
equation  (45)  with  /  =  I.  Calculation  of  the  non¬ 
linearity  parameter  is  obtained  from  measurements 
of  the  slope  and  the  fundamental  amplitude  e‘i . 

Acoustic  non-linearity  parameters 
and  material  properties 

The  acoustic  non-linearity  parameters  p  play  a 
central  role  in  determining  the  thermoelastic  prop¬ 
erties  of  crystalline  solids,  since  they  directly 
quantify  the  anharmonic  character  of  the  lattice 
modes.  Cantrell'18  has  shown  that  if  one  considers  a 
crystalline  solid  to  consist  of  a  large  number  of 
incoherent  non-linear  acoustic  radiation  sources 
identified  with  the  vibrating  particles  of  the  lattice, 
then  randomisation  of  the  resulting  acoustic  field 
together  with  the  assumption  of  a  stochastically 
independent,  fluctuating  radiation  field  at  the  abso¬ 
lute  zero  of  temperature  leads  to  an  exp^ssion  of 
the  temperature  dependent  rauiation  field  in  terms 
of  the  zero  point  field.  This  leads  directly  to 
expressions  of  the  thermodynamic  state  functions 
that  intrinsically  include  the  non-linearity 
parameters  as  a  direct  measure  of  structure  depen¬ 
dent  modal  anharmonicity.  These  'non-linear'  state 
functions  in  turn  lead  to  calculations  of  the  tem¬ 
perature  dependence  of  elastic  moduli49  as  well  as 
to  calculations  of  the  thermal  expansion  coeffi¬ 
cients48  directly  in  terms  of  the  non-linearity 
parameters. 

The  relationship  between  the  non-linearity 
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Table  2  Comparison  ot  structure,  uuiiumy,  anu 
.  acoustic  non-linearity  parameters  along 

[TOO]  direction  of  cubic  crystals 


Structure 

Bonding 

Range  of  I) 

NaC! 

Ionic 

>4  6 

1 4  0-  1 5-4 

bcc 

Metallic 

8  2 

7-4-8  8 

fee  'men.  gas) 

van  der  Waals 

5-4 

5-8-7  0 

fee 

Metallic 

5-6 

4-0-7  0 

Fluorite 

Ionic 

3-8 

3  4-4-6 

Zinc  blende 

Covalent 

2-2 

1-8-3  0 

parameters  and  fundamental  lattice  properties  is 
intrinsic.  Cantrell50  ha<:  shown  from  a  survey  of 
crystals  of  cubic  symmetry  that  for  a  given  acoustic 
mode  (i.e.  direction  of  acoustic  wave  propagation 
and  polarisation)  the  non-linearity  parameters  are 
found  to  be  ordered  according  to  the  crystalline 
structure.  Results  for  wave  propagation  along  the 
[100]  direction  in  cubic  crystals  are  given  in  Table  2. 
Listed  in  the  table  are  the  structure  of  the  crystal, 
the  type  of  atomic  bonding,  the  range  of  values  of  [3 
for  all  crystals  having  a  given  structure,  and  the 
average  values  of  the  j3s  in  that  range.  The  acoustic 
non-linearity  parameters  are  seen  to  be  stronglv 
ordered  according  to  the  type  of  crystalline  struc¬ 
ture.  The  range  of  values  of  (3  for  a  given  structure 
is  distinct;  overlap  of  ranges  occur  ■,  nly  slightly  for 
fee  and  fluorite  structures.  The  influence  of  the  type 
of  atomic  bonding  may  be  inferred  from  a  compari¬ 
son  of  the  fee  structured  crystals  The  /3s  for  the  fee 
metallic  bonded  crystals  and  the  fee  van  der  Waals 
bonded  crystals  are  approximately  equal  even 
though  the  difference  in  strength  of  these  bonds  is 
very  large.  It  is  inferred  that  the  influence  of  the 
bonding  on  the  value  of  the  .ion-linearitv  parameter 
is  small  compared  with  that  of  the  geometrical 
arrangement  of  atoms  defining  the  structure. 

In  order  to  explain  these  results  Cantrell."1 
•■allowing  an  approach  suggested  by  Ghate52  and  by 
.  liki  ana  Granaro,"  proposed  a  model  based  in 
t  st  approximation  on  a  short  range,  two  bodv. 
o  atral  force  potential  of  the  Born-Mayer  type.  He 
sh  wed  that  the  non-linearity  parameters  depend 
on  v  on  the  atomic  arrangement  of  the  crystal  and 
the  'hardness'  parameter  of  the  Born- Mayer 
port  itial.  He  also  found  that  the  hardness 
para  teter  together  with  the  atomic  nearest  neigh¬ 
bour  eparation  distance  determines  the  shape  or 
curvature  of  the  potential  curve.  The  dependence 
on  the  shape,  rather  than  the  magnitude,  of  the 
potenti  1  curve  provides  a  general  explanation  of 
the  resi.'ts  of  Table  2  including  the  insensitivity  of 
the  acoi  Stic  non-linearity  parameters  to  the  bond 
strength. 

While  the  model  successfully  explains  the  results 
of  Table  2  for  crystalline  solids,  it  does  not  readiiv 
predict  the  sign  and  magnitude  of  the  non-linearity 
parameter  ar  vitreous  silica,  an  amorphous  solid. 
The  /3  parar  ieter  for  vitreous  silica  is  negative54  in 
contrast  to  he  positive  values  of  the  parameter 
typically  men  ured  for  the  pure  mode  propagation 
directions  in  rystalline  solids.  In  addition.  Yost 
and  Cantrell4-  have  observed  what  thev  believe  to 
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induced  static  displacement  pulse,  (u,)  in  equation 
(43).  into  a  series  of  bulk  solitons  in  samples  of 
vitreous  silica.  The  large  negative  non-lineariiv 
parameter  gives  rise  to  a  contractive  static  pulse  in 
the  material  which  when  combined  with  the  large 
velocity  dispersion  in  the  medium  provides  condi¬ 
tions  suitable  for  the  generation  of  solitary  waves  or 
solitons.  The  observations  of  bulk  solitons  and 
acoustic  radiation  induced  static  contraction  in 
vitreous  silica  may  be  of  some  significance  in  the 
statistical  mechanical  treatment  of  amorphous 
structures.5"  a  subject  of  considerable  current  inter¬ 
est.  Vitreous  silica  is  known  to  have  a  negative 
thermal  expansion  coefficient  at  low  temperatures 
where  long  wavelength  vibrational  modes  dominate 
the  dynamical  properties.  At  higher  temperatures, 
where  the  expansivity  is  positive,  the  short 
wavelength  vibrational  modes  become  more  popu¬ 
lated  and  the  lattice  dynamics  is  dominated  by  a 
local  quartz-like  structure  having  positive  non¬ 
linearity  parameters  along  the  pure  mode  propa¬ 
gation  directions.  The  sign  of  the  thermal  expans¬ 
ivity  at  high  and  low  temperatures  is  thus  reflected 
in  the  sign  of  the  non-linearity  parameters  appro¬ 
priate  to  the  atomic  structure  ’seen’  by  the  domin¬ 
ate  lattice  vibrational  mode  at  that  temperature. 

The  significance  of  atomic  structure  and  its  effect 
on  the  non-linearity  parameter  is  not  diminished 
in  considerations  of  complex  materials.  Material 
microstructure  is  governed  fundamentally  by  even 
more  complicated  arrangements  of  atoms,  including 
random  or  amorphous-like  configurations  in  many 
cases,  which  also  affect  the  non-linearity  parameter. 
Of  particular  concern  here,  however,  is  the  fact  that 
the  mechanical  properties  of  many  engineering 
materials  are  derived,  at  least  in  part,  from  the 
presence  of  secondary  phases  in  the  solid  solution 
matrix.  The  presence  of  the  second  phase,  for 
example,  raises  the  flow  stress;  and  the  extent  of 
strengthening  depends  to  first  order  on  the  volume 
fraction,  size,  and  characteristics  of  the  second 
phase  precipitates  which  form  during  the  manufac¬ 
turing  process.  A  mathematical  model  has  been 
proposed  by  Cantrell  et  air 6  giving  the  effective 
non-linearity  parameter  of  an  alloy  as  a  function  of 
total  volume  fraction  of  second  phase  precipitates. 
Although  the  relationship  is  in  general  non-linear 
(see  Appendix  2  for  derivation)  the  equation  is 
approximated  to  within  experimental  error  for 
volume  fractions  of  second  phase  up  to  —10%  by 
the  linear  expression 

13  =  /5(1  -  Rfv)  . (48) 

where 


In  these  equations  /p  denotes  the  total  volume 
fraction  of  second  phase  precipitates;  the  CnS  are 
'quasi-isotropic'  second  order  elastic  constants 
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6  Graph  of  acoustic  non-linearity  parameter  of  7075 
(Af-Zn-Mg!  as  function  of  volume  fraction  of 
second  phase  precipitates 


(Voigt  notation).  The  overbar  denotes  the  values  of 
the  parameters  for  which  fp  =  0,  i.e.  pure  solid 
solution;  the  subscript  p  refer  to  values  for  the  total 
second  phase  precipitates  taken  collectively;  the 
subscript  i  refers  to  the  individual  constituents  of 
the  solid  solution;  the  superscript  c  refers  to  solid 
state  constituents;  and  e,  is  the  'depletion  constant' 
for  constituent  i  formed  from  solid  solution. 

Razvi  et  al.'1  provided  experimental  confirmation 
of  equation  (49)  for  measurements  of  fa  in  the  heat 
treatable  aluminium  alloy  7075  (A!  Zr,  -Mg).  The 
results  of  aeir  measurements  of  the  effective 
non-linearity  parameter  ji  as  a  function  of  volume 
fraction  of  second  phase  precipitates  /p  are  shown 
m  Fig.  h.  A  least  squares  fit  to  the  data  results  in  a 
linear  curve  with  a  1-0  correlation  coefficient  given 
by  the  sold  line  in  the  figure.  The  lack  of 
experiment.  I  measurements  of  the  non-linear  tv 
parameters  and  the  elastic  constants  of  individual 
second  ph  'Scs  of  Al  7075  prevents  an  exact  calcu¬ 
lation  of  the  effective  non-linearity  parameter  of 
the  materia'  directly  from  equation  (4$).  None  the 
less,  the  agreement  between  the  linear  form  of 
equation  (-5)  and  the  linearity  of  the  experimental 
data  is  evif  ent.  Measurements  of  the  intercept  and 
slope  of  the  curve  yield  a  value  of  0-29  for  the 
non-lineari  v  parameter  of  pure  solid  solution  of  Al 
7075  and  a  value  of  0T2  for  R. 

It  was  sumed  in  the  mathematical  model  that 
the  numbi  of  randomly  oriented  grains  contained 
within  a  p  thlemjth  of  the  propagating  sound  wave 
is  sufficien  iy  large  that  a  good  statistical  sampling 
of  quasi-is' tropic  behaviour  is  achieved.  It  was  thus 
expected  o  at  the  5  cm  pathlength  in  the  specimens 
allows  a  vide  variation  of  average  gr3in  size 
without  violating  the  quasi-isotropic  assumption.  A 
mamfestatit  n  of  proper  statistical  sampling  would 
be  the  invariance  of  the  non-linearity  parameter  as 
a  function  of  grain  size.  Figure  7  shows  the 
measured  non-linearity  parameters  in  Al  7075. 
corrected  for  attenuation,  as  a  function  of  aging  up 
to  237  h.  No  siemheant  variation  in  the  effective  fa 
occurs  although  the  average  gram  size  is  expected 
to  change.  It  is  inferred  from  these  results  and  trom 
Fig.  6  that  the  mathematical  model  qualitatively 
predicts  the  correct  variation  in  the  effective  non- 
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7  Graph  of  acoustic  non-linearity  parameter  of  7075 
(Al-Zn-Mg)  as  function  of  aging  time 


linearity  parameter  of  Al  7075  as  a  function  of 
volume  fraction  of  second  phase  precipitates. 

The  dependence  of  the  non-linearity  parameter 
on  the  second  phase  precipitates  for  alloys  and  the 
dependence  on  the  hardness  parameter  of  the 
Born-Mayer  potential  for  single  crystals  suggest  a 
possible  relationship  between  the  non-linearity 
parameter  and  the  engineering  hardness  number 
even  though  engineering  hardness  is  measured  by  a 
plastically  deforming  test  and  the  pertinent  atomic 
mechanisms  are  complicated  by  microstructural 
details.  None  the  less,  measurements  of  the  non¬ 
linearity  parameters  on  aluminium  alloys  as  a 
function  of  Rockwell  F  hardness  were  reported  by 
Li  el  al.-1  Their  results  are  shown  in  Fig.  8.  Similar 
measurements  on  iS%Ni  maraging  steel  as  a 
function  of  Rockwell  C  hardness  were  reported  by 
Yang  a  al.'*  and  are  shown  in  Fig.  9.  In  both  cases 
there  is  a  strong  correlation  between  the  non¬ 
linearity  parameters  and  the  hardness  number. 

Finally.  Elber59  reported  the  discovery  of  a  crack 
closure  phenomenon  that  occurs  with  metal  fatigue 
He  noted  that  closure  of  the  crack  planes  near  the 
crack  tip  can  occur  while  the  applied  stress  is  still 
tensile.  The  existence  of  a  crack  closure  stress 
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8  Plot  of  acoustic  non-linearity  parameter  of  Al 
alloys  v.  Rockwell  F  hardness 
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9  Plot  of  acoustic  non-linearity  parameter  of  18%  Ni 
maraging  steel  v.  Rockwell  C  hardness 


opens  the  way  tor  defining  an  effective  stress 
intensity  factor  la,  given  by 

la  =  (o-  ac,o,ure)(.T0)ia  . (50) 

where  a  is  the  applied  tensile  stress,  aC|0$urc  the 
crack  closure  stress,  and  b  the  crack  half  length. 
The  precision  of  l„  depends  on  how  well  one  can 
determine  aci0surc.  It  is  generally  difficult  to 
determine  aC|0,ure  experimentally,  since  conven¬ 
tional  crack  opening  determination  is  imprecise. 
However,  Yost  et  al T3  have  resolved  this  difficulty 
by  developing  a  technique  based  on  the  fact  that  as 
a  crack  in  a  compact  tension  specimen  is  opened 
acoustic  second  harmonics  are  generated  at  the 
(unbounded)  surfaces  of  the  crack  interface.  A 
typical  data  set  is  shown  in  Fig.  10  wiiere  the 


Tension  Load,  lb 

10  Graph  of  acoustic  second  harmonic  amplitude 
generated  at  unbounded  surface  of  crack 
interface  in  compact  tension  specimen  as  func¬ 
tion  of  tensile  load  (1  lb  e  0-454  kg) 


second  harmonic  generated  at  constant  fundamen¬ 
tal  amplitude  is  plotted  as  a  function  of  tension 
load.  Crack  opening  occurs  at  a  load  corresponding 
to  the  peak  of  the  generated  harmonic  amplitude. 


Stress-temperature  dependence  of 
sound  velocity 

The  variation  in  the  temperature  dependence  of  the 
sound  velocity  as  a  function  of  stress  is  the  basis  of 
a  new  method  of  characterising  residual  or  applied 
stress  fields  in  materials.  Consider  the  thermo¬ 
dynamic  tensions  (equation  (6))  to  be  a  function  of 
the  Lagrangian  strains  77^  and  temperature  T,  i.e. 
r, j  =  qj(rj, 7") .  The  differential  can  be  written 

dr,j  =  -  A.,d T  +  C,jlddr?k|  . (51) 

where  A,j  =  -(5^/57),,  is  defined  as  the  thermal 
stress  tensor.  For  constant  stress  (i.e.  dr^  =  0)  from 
equation  (51) 


where 


is  defined  as  the  thermal  strain  tensor.  For  isotropic 
solids  and  cubic  crystals 


. (54) 

where  rj  is  a  scalar  strain  parameter.  For  this  case 


h  j  =  ~P< 5ij 


(55) 


where  p  is  a  hydrostatic  pressure.  Thus 
equations  (52)— (55) 

dp  _  nrT 
ST  ~ 


from 


(56) 


where  K T  is  the  compressibility  and  aT  is  the 
thermal  expansion  coefficient.  From  equation  (56) 
it  can  be  inferred  that  a  temperature  dependent 
stress  is  induced  in  the  material  through  the  thermal 
expansion  coefficient.  For  relatively  small  changes 
in  the  temperature,  the  stress  p  is  linearly  depen¬ 
dent  on  temperature  and  from  equation  (IS)  is  seen 
to  give  rise  to  a  linear  variation  of  sound  velocity 
with  temperature. 

Salama  and  Ling61  considered  the  effect  of  stress 
on  the  temperature  dependence  of  the  ultrasonic 
velocity  of  aluminium  and  copper  alloys  by  experi¬ 
mentally  determining  the  relative  change  in  the 
temperature  derivative  of  the  natural  velocity 
(SWIST)  as  a  function  of  applied  stress  a.  A 
graphical  representation  of  their  results  is  shown  in 
Fig.  11.  They  find  that  the  linear  curves  of  Fig.  11 
3re  represented  quite  well  by  the  empirical 
equation 


where  (cWi3T)0  is  the  temperature  derivative  of 
sound  velocity  at  zero  stress,  ( c\VI3T)a  is  that  at 
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1 1  Plot  of  temperature  dependence  of  acoustic 
natural  velocity  as  function  of  applied  stress 


applied  stress,  o  is  the  applied  stress  (compressional 
in  their  experiments),  and  K  is  a  proportionality- 
constant  equal  to  2-4  x  10“’’  MPa"1  for  aluminium 
alloys  and  2-5  x  lO"4  MPa”1  for  copper  alloys. 
These  results  were  confirmed  in  steel  and  other 
aluminium  alleys  by  Salama  ei  al62'*5 

A  preliminary  calculation  of  the  temperature 
dependence  of  the  sound  velocity  at  zero  initial 
stress  was  reported  by  Salama.66  He  also  obtained  a 
comparison  between  the  calculated  stress  derivative 
of  the  temperature  dependence  at  zero  initial  stress 
and  experimental  values  of  the  stress  derivative. 
The  results  yielded  reasonable  values  for  the  stress 


UCpCllUClilCb  UWU  ICO  IVJ  COilviuJiO.i  •■»»>-*» 

temperature  dependence  of  ultrasonic  velocity  in 
metallic  alloys  varies  linearly  with  stress  as  ob¬ 
served  experimentally.  A  similar  calculation  and 
conclusion  were  obtained  independently  by  Chern 


More  recently,  Cantrell0*’  obtained  an  analytical 
expression  for  the  constant  K  by  considering  the 
explicit  temperature  dependence  of  the  elastic 
moduli  via  the  quasiharmonic-anisotropic- 
eontinuum  model  of  crystals.69  The  effect  of 
applied  or  residual  stress  was  made  explicit  by 
expanding  the  elastic  moduli  about  the  state  of  zero 
stress  in  a  manner  analogous  to  that  leading  to 
equation  (18).  The  generalised  proportionality  con¬ 
stant  Kpq,  analogous  to  K  in  equation  (57),  for 
arbitrary  applied  stress  crpq  is  found  to  be 
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■vhere  kB  is  Boltzmann’s  constant,  with 

rf"  =  -  [2 PuW'-U.U,! 

2  p0W- 

-f-  C0^rnunv  (7U  Lj  v)A  mA  n  ] 


—  =  2  yfyf 


(  (-'n/lydmn 


Srjyd  2po\V~ 

A~  d  a(5yc5munv  Uu  Uv  "h  ^Ccy^rnunv Tu  L,})  A/m  A' r, 
. (62) 

The  experimental  measurements  of  Salama  and  Ling 
used  a  longitudinal  ultrasonic  wave  propagating 
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in  a  direction  perpendicular  to  a  uniformly  applied 
cdmpressional  static  stress.  Their  work  was  per¬ 
formed  on  quasi-isotropic  polycrystalline  solids.  For 
such  an  experiment  situation  the  generalised 
expression  for  Kpq.  equation  (58).  reduces  to08 


. (63) 

Note  that  equations  (5S)  and  (60)  involve  the 
temperature  derivatives  of  third  order  elastic  con¬ 
stants.  From  equations  (61)  and  (62)  it  can  be  seen 
that  such  derivatives  may  be  written  explicitly  in 
terms  of  fourth  and  fifth  order  elastic  constants 
when  the  quasiharmonic-anisotropic-continuum 
model  is  used  Conversely,  the  experimental  meas¬ 
urements  of  K  provide  a  method  of  determining  the 
fifth  order  elastic  constants  of  the  material  in  this 
approximation  provided  the  fourth  order  constants 
can  be  determined  from  other  methods.  For  pure 
aluminium  and  copper  single  crystals  the  K  con¬ 


stants  calculated  for  uniaxial  stress  applied  along 
[001  ]  and  compressional  wave  propagation  along 
[100]  in  these  materials  are  1-8  x  10"  J  MPa-1  for 
aluminium  and  3  0-x  lfl'J  MPa-1  for  copper. 
These  theoretical  values  are  in  very  good 
agreement  with  the  experimental  values  obtained 
by  Salama  and  Ling  for  corresponding  alloys  of 
these  elemental  metals  and  suggest  that  the  base 
metal  of  a  given  alloy  system  dominates  the 
stress-temperature  of  the  sound  velocity  in  the 
alloys. 

Low  field  magnetoacoustics 

A  novel  technique  for  assessing  stresses  in  ferro¬ 
magnetic  materials  has  been  developed  by  Nam- 
kung  ei  al.‘°  The  technique  is  based  on  the 
variation  of  the  sound  velocity  measured  as  a 
function  of  low  field  magnetisation  of  the  material. 
The  normal  crystalline  structure  of  a  solid  is 
generally  altered  when  the  solid  is  ferromagneti- 
callv  ordered.  For  solids  such  as  iron  the  cubic  urn 
cells  are  spontaneously  deformed  into  tetragonal 
structures  with  the  longer  edges  aligned  along  the 
magnetisation  vector.  The  magnetisation  vectors  in 
turn  are  oriented  along  one  of  the  six  equivalent 
crystallographic  (100)  directions.  A  two  dimension¬ 
al  representation  of  the  net  domain  structure  for 
the  case  of  zero  net  magnetisation  and  zero  initial 
stress  is  shown  at  the  left  in  Fig.  12a.  The  effect  of  a 
residual  or  applied  stress  a  on  the  ferromagnetic 
state  is  to  alter  the  energy  density  of  a  given 


-•?  effect  of  :ero  initial  stress.  ^  effect  of  tensile  initial  stress:  c  effect  of  corror^ssive  initial  stress 

12  Two  dimensional  representation  of  net  magnetic  domain  structure 
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+  a^Qi/jv,)  . (64) 

where  £')'  is  the  first  anisotropy  constant,  A"^  the 
saturation  magnetostriction  along  the  (100)  axis. 
An i  that  along  the  (111)  axis,  and  the  ais  and  yiS 
are  direction  cosines  of  the  magnetisation  vector 
and  uniaxial  stress  axis,  respectively,  with  respect  to 
the  cube  axes. 

For  transdomain  stresses  the  magnetoeiastic 
energy  densities  are  generally  different  for  adjacent 
domains  and  give  rise  to  a  net  pressure  acting  on 
the  domain  walls  separating  them.  From  equation 
(64)  the  net  pressure  (i.e.  energy  density)  is 
non-zero  only  for  those  adjacent  domains  whose 
magnetisation  vectors  are  at  right  angles  to  each 
other  (i.e.  90°  domain  walls).  Under  tension  the  90° 
domain  wall  motion  is  such  that  the  net  volume  of 
domains  oriented  along  the  uniaxial  stress  axis  is 
increased  as  indicated  by  the  domain  representation 
at  the  left  of  Fig.  12 b.  Under  compression  the 
domains  orient  close  to  a  plane  perpendicular  to  the 
stress  axis  as  shown  at  the  left  of  Fig.  12c.  In  both 
cases  the  effect  of  stress  is  to  decrease  the  area  of 
90°  walls.  The  stress  induced  903  wall  motion, 
however,  is  generally  restricted  bv  the  interaction 
between  the  walls  and  various  lattice  defects. 
Hence,  complete  domain  alignment  under  stress 
does  not  generally  occur  in  the  case  of  impure 
ferromagnetic  materials  like  steels. 

When  an  external  magnetic  field  is  appued  to  the 
ferromagnetic  solid,  domains  with  low  Zeeman 
energy  become  seed  domains  and  their  volume 
begins  to  expand.  In  this  case  both  ISO3  and  90° 
wall  motion  occurs  such  that  as  the  applied  field 
increases  more  domains  become  oriented  along  the 
applied  field  direction.  The  case  for  zero  initial 
stress  in  the  material  ;c  illustrated  in  Fig.  12a.  The 
relative  change  in  net  domain  structure  as  a 
function  of  applied  magnetic  field  is  dependent  on 
the  domain  structure  initially  induced  by  the  stress 
at  zero  magnetic  field  conditions  as  illustrated  in 


13  Predicted  changes  in  sound  velocity  as  function 
of  magnetic  field  strength  for  various  states  of 
initial  stress 


Fig.  12/t  and  c.  The  applied  magnetic  field  strength 
is  shown  to  increase  from  zero  at  the  left  to  some 
maximum  value  indicated  at  the  right  in  the  figure. 

The  effect  of  net  domain  orientation  on  the 
sound  velocity  in  the  ferromagnetic  material  is 
governed  by  the  relative  amount  of  90°  domain 
walls  available  in  the  material  at  the  time  of  the 
velocity  measurement.  The  elastic  modulus  £  of 
ferromagnetic  materials  is  given  by 

£  =  — . .  . (65) 

■  £mc 

where  o  is  the  applied  stress,  aC|  the  linear  elastic 
strain,  and  cme  the  magnetoeiastic  strain.  The 
magnetoeiastic  strain  is  generally  assumed  to  be 
proportional  to  the  total  area  of  90°  domain  walls. 
From  equation  (65)  and  Fig.  12,  the  expected 
functional  change  in  sound  velocity  as  a  function  of 
applied  magnetic  field  strength  H  for  various  states 
of  initial  stress  is  qualitatively  shown  in  Fig.  13.  It  is 
of  interest  to  note  that  tensile  and  compressive 
stresses  give  rise  to  initial  slopes  of  opposm- 
polarity  and  thus  provide  a  means  of  testing  the 
sign  of  the  initial  stress. 


a  norvcompressive  stresses:  b  compressive  stresses 

14  Changes  in  sound  velocity  in  AISI  1020  steel  as 
function  of  magnetic  field  strength  for  various 
states  of  initial  stress;  PPM  parts  per  million 


International  Mater,ais  Reviews  1991  VoL  36  Mo.  4 


Namkung  et  alP'112  have  performed  an  exten¬ 
sive  series  of  natural  velocity  measurements  in 
various  steel  samples  using  the  pulsed  phase-locked 
loop  technique.-'4  Their  results7-'  for  AISI  1020  steel 
are  shown  in  Fig.  14.  The  curves  are  in-  nominal 
agreement  with  the  model  curves  in  Fig.  13.  Their 
results  °'7~  for  railroad  steel,  however,  show  some 
deviation  from  the  model.  Specifically,  the  curves  at 
tension  indicate  a  greater  modulus  or  velocity 
change  than  the  curve  at  zero  stress.  In  addition, 
the  unstressed  curve  is  relatively  flat  until  a  rather 
large  magnetisation  is  induced  in  the  sample. 
Namkung  et  alP  explain  that  the  deviation  results 
from  the  high  degree  of  local  lattice  strain  that 
significantly  impedes  the  motion  of  90’  walls  but  has 
little  influence  on  the  180°  wall  motion.  While  the 
180°  wall  motion  contributes  to  the  magnetisation 
of  the  material,  only  the  90°  wall  motion  affects  the 
elastic  modulus  and  velocity  changes.  The  unstres¬ 
sed  curve  is  thus  expected  to  remain  relatively  fiat 
until  high  field  strengths  are  achieved.  The  velocitv 
increase  at  higher  fields  is  caused  mainly  by  domain 
rotation  and  in  part  by  90°  wall  motion.  The  effect 
of  initial  tension  is  to  aid  the  magnetic  field  in 
moving  90°  walls  at  a  somewhat  reduced  field 
strength.  This  results  in  an  increasingly  upward  shift 
of  the  tension  curves  as  the  initial  tension  increases. 

Conclusions 

It  is  clear  that  non-linear  acoustic  methods  can  be 
used  to  obtain  useful  and  in  some  cases  unique 
information  about  materials.  Although  the  same 
basic  experimental  acoustic  techniques  are  applic¬ 
able  to  most  materials,  the  meaning  and  inter¬ 
pretation  of  the  experimental  results  are  limited  by 
the  validity  of  the  acoustic  models  used  to  describe 
the  material.  Most  measurements  of  non-linear 
thermoelastic  properties  have  been  performed  on 
materials  having  a  relatively  simple,  well  defined, 
single  phase  structure  in  which  elastic  wave  behav¬ 
iour  is  understood  from  the  generalised  wave 
equations  for  solids  of  arbitrary  crystalline  symmet¬ 
ry.  Indeed,  elastic  constant  data  (particularly  the 
higher  order  elastic  constants)  have  been  reported 
mostly  on  relatively  defect  free  single  crystals. 
Multiphase  and  composite  materials,  although  more 
complex  structurally,  have  generally  been  treated 
acoustically  as  quasi-isotropic  or  at  best  mildly 
textured  materials  without  regard  to  their  micro- 
structural  composition.  Although  acoustic  wave 
propagation  can  in  principle  be  treated  as  taking 
place  in  such  materials  with  effective  elastic  moduli 
that  reflect  whatever  degree  of  symmetry’  the 
researcher  feels  necessary  to  consider,  the  compo¬ 
sitional  character  of  complex  materials  in  itself  is 
important  and  deserves  specific  consideration. 

An  attempt  has  been  made  in  this  review  to 
redress  the  general  lack  of  compositional  consider¬ 
ation  by  emphasising  the  interrelationship  of 
material  composition,  material  properties,  and 
acoustic  measurements.  Although  the  acoustic 
models  presented  here  in  terms  of  volume  fraction 
of  material  constituents  are  in  some  cases  simplistic, 
the  essential  features  predicted  by  the  models  are 


borne  out  in  experiment.  The  models,  indeed,  do 
provide  a  more  comprehensive  understanding  of 
the  influence  of  material  composition  on  acoustic 
properties  and.  hopefully,  will  serve  as  catalysts  for 
the  development  of  better,  more  comprehensive, 
models  of  compositional  influence  on  wave  behav¬ 
iour.  The  success  of  many  non-destructive 
evaluation  efforts  and  the  evolution  of  new  non¬ 
destructive  methodologies  rests  in  large  part  on  the 
ability  to  model  comprehensively  such  microstruc- 
tural  influences  on  acoustic  properties.  The  con¬ 
sideration  of  compositional  influence,  including 
magnetic  or  other  domain  effects  where  appro¬ 
priate.  and  the  linking  of  acoustic  non-linearity  to 
the  basic  lattice  dynamical  behaviour  of  materials 
are  necessary  to  provide  the  proper  science  base  for 
advancement  and  innovation  in  non-destructivelv 
assessing,  characterising,  and  testing  the  thermo- 
elastic  and  mechanical  properties  of  materials. 
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APPENDIX  1 

Dynamic  moduli  and  elastic  constants 

The  mechanical  properties  of  isotropic  solids  are 
usually  obtained  from  static  or  low  frequency 
dynamic  testing  of  specimens  having  a  specified 
geometrical  shape.  Here  the  relationships  between 
the  elastic  moduli  measured  in  such  tests  and  the 
elastic  constants  (specifically  Lame  constants) 
obtained  from  ultrasonic  measurements  of  the  same 
material  are  summarised.  The  relationship  between 
the  Lame  constants  and  the  Brugger  second  order 
elastic  constants,  generally  used  in  the  text,  is  given 
in  equation  (TO). 


i  iiC  lineal  bllebb- ollaiu  i cianuiLtiup  i \_»vjcw_  :> 

law)  for  an  isotropic  solid  in  teims  of  the  Lame 
constants  a  and  p  can  be  written  in  the  form 

a,  =  2 -f  kOS,,  . (66) 

where  6  is  the  dilation  or  change  in  material  volume 
defined  by 

6  -  £U  +  E22  +  fc'33  . (67) 


it*  Cv.{UUU\Jii  ‘ - -  I  ■  i.is-  j_.  •  - - 

can  be  considered  effectively  unbounded  when  the 
ultrasonic  frequencies  are  sufficiently  large  that  the 
acoustic  wavelength  is  small  compared  with  the 
dimensions  of  the  sample. 

Poisson's  ratio  v  is  also  defined  from  static 
conditions  as  the  ratio  of  the  lateral  strain  to  the 
longitudinal  strain  and  is  obtained  from  equation 
(70)  as 


and  c,j  are  the  linear  strains  defined  in  terms  of  the 
displacement  gradients  (oujda^  as 


T,  =  T 

'  '  C7C1 


(68) 


Equation  (66)  is  obtained  directly  from  equations 
(5),  (7),  and  (20)  where  only  the  linear  terms  are 
retained  in  the  final  expression.  The  p  constant  is 
often  referred  to  as  the  shear  or  rigidity  modulus 
and  is  given  the  symbol  G  by  many  authors.  The 
Lame  constants  can  be  obtained  from  meas¬ 
urements  of  the  compressional  and  shear  bulk  wave 
velocities  of  the  unstressed  (i.e.  a  =  0)  sample  as 
indicated  by  equations  (23). 

The  bulk  modulus  k  is  an  important  engineering 
modulus  defined  from  a  static  measurement  of  the 
change  in  the  material  volume  resulting  from  an 
applied  hydrostatic  pressure  p  as 


The  right-hand  side  of  equation  (69)  is  obtained  by- 
writing  (7,,  =  -p<5,,  for  the  hvdrostatic  pressure  and 
using  this  expression  in  equation  (66).  As  indicated 
in  equation  (69),  the  bulk  modulus  can  also  be 
calculated  (torn  knowledge  of  the  Lame  constants 
determined  from  ultrasonic  bulk  wave  meas¬ 
urements. 

Consider  now  a  cylindrical,  isotropic  solid  and 
deform  the  solid  statically  along  the  cylindrical  axis, 
which  is  defined  as  the  x  or  1-axis.  The  surface  of 
the  cylinder  is  stress  free  and  so  from  equation  (66) 

<7i  |  =  2pft  |  +  kO  : 

0  -  2 pr::  +  kd  . .  .  (70) 

0  ==  2/XS33  -T  kd 

Young's  modulus  £  is  defined  as  the  ratio  of  the 
applied  axial  stress  aM  to  the  longitudinal  strain  fu 
and  from  equation  (70) 


Ou  _  Ml  3A  +  2 p.) 

£l  l  k  ~  fX 


(71) 


Although  equation  (71)  is  defined  from  static 
conditions.  Young's  modulus  can  be  obtained  from 
measurements  of  the  velocity  vE  of  low-  frequency- 
longitudinal  waves  propagating  along  the  cylinder 
axis  as 


E  -  pov~f 


(72) 


It  is  important  to  emphasise  that  the  velocity  vE 
measured  under  the  above  conditions  is  not  equal 
to  the  velocities  measured  in  unbounded  media,  as 


£  1 1  2(A  +  p) 

It  is  clear  from  the  above  relationships  that  the 
various  moduli  are  intrinsically  related.  For 
example,  it  is  easily  shown  that 


(1  +  v)(  1  -  2v) 


2(1  +  v) 


APPENDIX  2 


Non-linearity  parameters  of 
multiphase  alloys 

A  mathematical  model  was  developed  by  Cantrell 
er  a/.56  to  describe  the  effective  non-linearity 
parameters  of  multiphase  alloys.  The  model  is 
based  on  the  law  of  mixtures  for  the  non-linearity 
parameter  and  is  reproduced  here  to  underscore  the 
essential  features  and  assumptions  not  brought  out 
in  the  text. 

Consider  a  rotation  of  the  Cartesian  reference 
frame  such  that  the  direction  of  acoustic  wave 
propagation  is  always  along  the  x  or  1-axis  of  the 
rotated  coordinate  system.  Then  perform  a  second 
orthogonal  transformation  on  the  system  which 
diagonalis^s  the  non-linear  equations  o'"  motion, 
equation  (41),  such  that  the  equations  are  decoup¬ 
led  into  three  independent  equations  corresponding 
to  the  three  wave  polarisation  directions.  Perfor 
ming  these  same  transformations  on  the  general 
stress-strain  relationship  of  equation  (5)  and  keep¬ 
ing  only  terms  to  first  order  in  the  non-linearity  lead 
to  the  relationship  between  the  transformed  stress 
apq  and  the  transformed  displacement  gradients 
(SiijISa ,)  given  by11  (no  sum  on  /) 


CU; 

-a  i 


(76) 


where  p”  and  v^Cre  linear  combinations  of  second 
and  third  order  elastic  constants,  is  the  La- 
grangian  coordinate  along  the  x  or  l-axis,  and 
j  =  1.2,3  is  a  mode  index  representing  the  appro¬ 
priate  polarisation  direction.  Consideration  is 
restricted  to  longitudinal  waves  (j  =  1)  in  quasi 
isotropic  solids  (i.e.  solids  consisting  of  random-.' 
oriented  grains).  Thus  the  superscripts  arm  r  n 
scripts  are  dropped  and  the  longitudinal  stresses  a 
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written 


equality  of  stresses,  we  obtain 


o  ~  C,,(— j  -  -  Cn/il  )  .  .  .  ( / 7) 

'  3a  ■'  2  '  --a' 

where  C'n  is  the  longitudinal  second  order  elastic 
constant  written  in  Voigt  notation  and  />  is  the 
acoustic  non-linearity  parameter.  Solving  for 
( 3ui3a )  in  terms  of  a  we  obtain 

/  -’w  )  5  1  P  - 

W  C'H  2  Cl, 

Consider  now  the  solid  to  consist  of  any  number  of 
phases  S’.  Assume  that  for  a  given  phase  i,  the 
grain  orientations  are  perfectly  randor  (no  texture) 
and  that  the  number  such  grains  contained  within  a 
path  length  of  sound  is  sufficiently  large  to  provide 
a  good  statistical  sampling.  To  the  extent  that  such 
conditions  are  maintained  the  value  of  tire  effective 
non-linearity  parameter  is  expected  to  be  indepen¬ 
dent  of  grain  size. 

In  order  to  obtain  the  appropriate  mixing  law  for 
/>  detine  Vo  and  p0  to  be  the  initial  (unperturbed) 
volume  and  mass  density,  respectively,  of  the  solid. 
The  local  transformation  from  the  initial  state  to 
the  deformed  state  V  or  p  is  defined  through  the 
Jacobian 


The  volume  at  any  time  is  considered  to  consist  of 
a  number  .V  of  constituent  phases  i  such  that 


Then  t';om  equations  (79)  and  (SU) 

j  -  -7-Xp  =  rr  2m  =  2jj\ 


where  7,  (WWi)  is  the  Jacobian  for  phase  i  and 
/,  ---•  (KVV'n)  is  the  volume  fraction  of  phase  i. 
Expanding  the  Jacobian  in  terms  of  the 
displacement  gradients  (3uJ3a.)  =  u,t  and  keeping 
the  linear  terms  (small  strains!,  then  (Einstein 
summation ) 

J  1  i-  m  ,  . (82) 


Substituting  equation  (82)  into  equation  (81)  gives 


Consider  now  the  quasi-isotropic  solid  and  assume 
that  for  a  given  phase  i  the  grain  orientations  are 
sufficiently  random  and  of  sufficiently  large  number 
that  each  phase  responds  individually  as  an  isotro¬ 
pic  structure.  Under  such  conditions,  equation  (83) 
may  be  written 


C)  -  2(4“)/ 


1  1  p 

—  o  - - - —  a~ 

Cu  2  Ci, 


=  (s4-/.)» 


7(0 1): 


Equating  like  powers  of  a.  we  obtain 

±  =  2—f,  . 

c„  .  o. 


P=C2n2^~f i  . (87) 

In  general,  a  non-linear  relationship  is  found 
between  the  effective  non-linearity  parameter  p  and 
the  volume  fraction  f\  of  individual  phases  because 
of  the  appearance  of  Cn  in  equation  (87).  It  is  of 
interest  to  point  out  that  for  liquid  media  the 
second  order  elastic  constants  Cu  =  0  and 
On  =  0)2.  In  this  case  (Cii)-1  can  be  identified 
with  the  liquid  state  compressibilities  and  equation 
(S7)  becomes  identical  to  the  results  of  Apfel'4  for 
immiscible  liquid  mixtures. 

Assume  now  that  the  solid  consists  of  any 
number  of  distinct  second  phase  precipitates  and 
that  the  relative  volume  fractions  of  constituent 
second  phase  precipitates  for  a  given  alloy  are 
constant  -  only  the  total  volume  fraction  of  second 
phase  precipitates  is  considered  to  change.  Hence, 
the  effective  non-linearity  parameter  pp  and  the 
effective  (Ci;)p  of  the  second  phase  precipitates 
taken  collectively  remain  unchanged.  The  inva¬ 
riance  of  the  relative  volume  fractions  of  second 
phase  precipitates  must  necessarily  come  at  the 
expense  of  the  solid  solution  constituents.  It  is 
assumed  that  the  depletion  of  the  constituents  of 
solid  solution  occurs  linearly  as 

n  =  fl  -  e.fp  . (SS; 


in  the  notation  of  equation  (77).  From  equations 
(78)  and  (79),  together  with  the  assumption  of  local 


where  f;  is  the  present  volume  fraction  of  solid 
solution  constituent  /,  /p  is  the  total  volume  fraction 
of  second  phase  precipitates,  /f  is  the  voiurm 
fraction  of  constituent  i  in  pure  solid  solution  (i.e. 
when  /p  =  0),  and  e,  is  the  'depletion'  constant  for 
constituent  i. 

From  the  above  considerations  equation  (86)  can 
be  written 

—  =  2  — —fi  +  — —  /p  ....  (89) 

Cu  •  (Ci,),'  (Cu)p  P 

where  in  equation  (89)  and  in  all  following  equa¬ 
tions  the  summation  is  taken  over  solid  solution 
constituents  o.’.ly,  the  superscript  c  refers  to  solid 
solution  constituents,  and  the  subscript  p  refers  to 
the  second  phase  constituents  taken  collectively.  It 
has  been  assumed,  somewhat  tenuously,  in  writing 
equat'on  (89)  that  the  solid  oolution  constituents 
behave  as  an  immiscible  mixture.  While  this 
assumption  is  not  strictly  true,  the  summation  term 
in  the  equation  is.  none  the  less,  representative  of 
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lions  (88)  and  (S')) 

i  _  1  .  r  i  _  v  i 

til  C  I  I  -(  C|  [);,  I  (C, :  ),  - 


K'lii  Cv.|ua- 


solution.  For  typical  values  of  the  Cu*  and  /is  the 
coefticients  of  the  terms  containing  fp  in  equation 
(92)  are  estimated  to  be  of  order  unity.  Expanding 
equation  (92)  in  a  power  series  for  small  values  of 
/.,  and  keeping  only  the  linear  terms,  gives 

(i  -  />’(  1  -  Rfv)  . (94) 


where  the  constant 


is  the  solid  solution  contribution. 

Similarly,  from  equations  (87),  (S8),  and  (90)  it  is 
found  that  the  effective  non-linearity  parameter  /3 
of  the  solid  in  terms  of  total  volume  fraction  of 
second  phase  precipitates  fp  is  given  by 


(92) 


where 


Equation  (95)  is  a  linear  approximation  to  equation 
(92).  A  survey  of  typical  values  of  (Cu)i  and  $ 
indicates  that  equation  (95)  should  be  accurate  for 
most  materials  to  within  typical  experimental 
uncertainty  for  volume  fractions  as  high  as  approxi¬ 
mately  10%. 


(93) 
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ABSTRACT 

Conventionally,  metal  matrix  composites  (MMC)  are  reinforced  with  either  particles  or  fibers. 
Recently,  a  new  class  of  composites  where  a  mixture  of  panicles  and  fibers  is  used  as 
reinforcement  has  emerged.  The  particles  improve  the  isotropic  mechanical  and  thermal 
properties,  whereas  the  fibers  introduce  directionally  favorable  properties  for  specific 
applications  of  the  material.  The  elastic  properties  of  three  different  matrix  alloys  containing  6% 
alumina  fibers  and  varying  alumina  particle  volume  fractions  of  9,  13,  and  17%  have  been 
determined  using  ultrasonic  velocity  measurements.  The  results  show  that  the  elastic  moduli 
increase  with  the  particle  content  and  the  composites  have  the  highest  elastic  stiffness  in  the 
directions  of  the  fiber  plane.  A  model  is  developed  to  explain  the  observed  elastic  moduli  of  this 
type  of  composites.  The  model  uses  results  of  the  theories  by  Ledbetter  and  Datta  for  spherical 
inclusions  and  Hashin  and  Rosen  for  aligned  fibers.  Furthermore,  it  includes  an  averaging 
procedure  suggested  by  Christensen  and  Waals.  The  agreement  between  measured  and 
calculated  elastic  moduli  is  found  to  be  good.  In  a  second  series  of  measurements,  the  elastic 
moduli  in  two  sets  of  extruded  MMCs  and  one  set  of  pressed  MMCs  are  determined.  These 
composites  are  reinforced  with  silicon  carbide  particles.  Their  elastic  moduli  as  well  as  their 
elastic  anisotropies  are  explair  ,d  using  the  theories  discussed  earlier. 


INTRODUCTION 

In  order  to  increase  structural  efficiencies  in  modem  design,  materials  possessing  high  stiffness 
and  high  strength  are  required.  One  class  of  engineering  materials  fulfilling  these  requirements 
are  metal -matrix  composites  (MMCs).  In  these  composites,  properties  of  the  material  can  be 
tailored  by  the  appropriate  selection  of  matrix  and  reinforcement  materials  and  by  their  mutual 
arrangement  in  order  to  meet  specific  needs  of  the  designed  component.  The  matrix  and  the 
reinforcement  are  to  be  selected  so  that  they  combine  their  different  mechanical  and  elastic 
properties  in  a  synergistic  way.  Also,  micro  structures  resulting  from  different  fabrication 
processes  are  found  to  influence  properties  of  these  composites  and  provide  valuable 
information  for  their  further  development. 

Many  models  have  been  developed  to  determine  the  effective  elastic  moduli  of  composite 
materials.  Most  of  these  models  deal  with  reinforcement  in  the  form  of  spherical  particles 
(Ledbetter  and  Datta,  1986,  Budiansky,  1965),  ellipsoidal  inclusions  (Eshelby,  1957,  Chow, 
1977)  or  infinitely  long  fibers  (Hashin  and  Rosen,  1964,  Hill,  1964).  Ledbetter  and  Datta  used 
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a  multiple  scattering  theory  to  predict  the  elastic  behavior  of  composites  with  a 
nonhomogeneous  panicle  distribution.  In  the  model  they  assume  that  the  particles  together  with 
the  matrix  form  an  enriched  "sea"  that  surrounds  "islands"  of  pure  matrix  material.  These  non- 
spherical  islands  are  aligned  and  produce  anisotropy.  On  the  other  hand,  in  order  to  determine 
the  effective  elastic  properties  of  fiber  reinforced  materials,  Hashin  and  Rosen  introduced  the 
composite  cylinders  model.  In  this  model,  the  composite  is  considered  to  be  comprised  of 
infinitely  long  circular  cylinders  embedded  in  a  continuous  matrix  phase.  Each  fiber  has  a 
radius,  a,  which  is  surrounded  by  an  annulus  of  matrix  material  of  radius  b,  and  the  ratio  alb  is 
considered  to  be  constant  for  all  composite  cylinders.  In  order  to  obtain  a  volume  filling 
configuration,  the  absolute  size  of  the  cylinders  must  vary  considerably. 

In  some  industrial  applications  composites  used  are  of  a  more  complex  nature,  where  a  mixture 
of  particles  and  fibers  is  used  as  a  second  phase.  To  our  knowledge  no  models  that  describe  the 
elastic  properties  of  these  composites  are  available.  The  present  study  is  concerned  with 
composites  reinforced  with  low  values  of  both  particles  and  fibers,  and  hence,  no  interaction 
between  fibers  and  particles  is  assumed.  In  modeling  these  composites,  we  first  consider  the 
matrix  material  and  the  particles  to  form  an  effective  matrix.  Since  the  particles  are 
homogeneously  distributed  in  the  metal,  the  effective  matrix  is  considered  to  be  homogeneous 
(Ledbetter  and  Datta,  1986).  The  effective  matrix  is  then  considered  to  be  reinforced  with  fibers 
which  are  randomly  oriented  in  one  plane.  The  influence  of  the  fibers  on  the  elastic  moduli  of 
the  composites  is  then  determined  first  by  using  the  composite  cylinders  model  for  an  aligned 
fiber  system  (Hashin  and  Rosen,  1964)  and  second  by  preforming  a  geometric  average 
procedure  (Christensen  and  Waals,  1972)  which  takes  care  of  the  2-dimensionally  random 
orientation  of  the  fibers.  The  results  obtained  show  a  good  agreement  between  calculations  and 
measurements  where  details  are  given  elsewhere  (Grelsson  and  Salama,  1990). 


EXPERIMENTAL 

The  metal  matrix  composites  used  in  this  investigation  consist  of  aluminum  alloys  as  the  matrix 
material  and  either  SiC-particles  or  alumina  fibers  and  particles  as  the  reinforcement.  The 
composites  based  on  AI  7064  and  A1  8091  were  obtained  as  extruded  rods  of  25  mm  in 
diameter.  The  specimens  based  on  Al  6061  were  received  as  pressed  plates  of  the  dimension  6  x 
25  x  60  mm.  The  alumina  reinforced  specimens  were  manufactured  by  squeeze  casting  and 
received  as  bars  of  the  dimension  12  x  12  x  50  mm.  The  volume  percentage  of  reinforcement  in 
the  different  sets  of  MMCs  is  shown  in  Table  1.  Specimens  are  cut  from  the  as-received 

Table  1 .  Manufacturing  method,  matrix  alloy  and  volume  fraction  of 
reinforcement  of  the  MMCs  investigated. 


Manufacturing  method 

Matrix  alloy  and  volume  fraction 
of  reinforcement 

Extruded  rods 

Al  7064 

+ 

0% 

SiC 

Al  7064 

+ 

15% 

SiC 

Al  7064 

+ 

20% 

SiC 

Al  8091 

+ 

0% 

SiC 

Al  8091 

+ 

10% 

SiC 

Al  8091 

+ 

15% 

SiC 

Pressed  plates 

Al  6061 

+ 

0% 

SiC 

Al  6061 

+ 

25% 

SiC 

Al  6061 

+ 

40% 

SiC 

Squeeze  cast  bars  A 132  +  23%  AI2O3 

A13  +  19%  AI2O3 

A 1100  +  15%  AI2O3 
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Extrusion 


Extruded  A1 7064  and 
Al  8091  MMC  rods 


Pressed  A1  606 1  Squeeze  cast 

MMC  plates  MMC  bars 


Fig.  1.  Specimen  geometry  and  designated  coordinate  systems. 


composites  and  the  faces  of  each  specimen  are  machined  flat  and  parallel  to  within  ±25  Jim.  The 
coordinate  systems  are  chosen  such  that  in  extruded  samples  the  X)-  and  X2-axes  are 
perpendicular  to  the  extrusion  direction.  In  the  pressed  specimens  the  Xi-  and  X2-axes  are  in 
the  plate  perpendicular  to  each  other  and  the  X3*axis  is  oriented  parallel  to  the  compression 
direction.  In  the  squeeze  cast  samples  the  Xi-  and  X2 -directions  are  in  the  fiber  rich  plane  at 
right  angles  to  each  other  and  the  X3-direction  is  normal  to  the  fiber  rich  plane.  Fig.  1  displays 
schematically  the  cut  of  the  specimens,  the  manufacturing  parameters,  and  the  orientation  of  the 
coordinate  systems.  For  all  specimens  the  microstructure  is  examined  in  the  three  orthogonal 
planes  in  terms  of  particle  size  distribution  and  area  fraction  covered  by  the  reinforcement. 
Measurements  of  ultrasonic  velocities  were  performed  using  the  pulse-echo-overlap  method 
with  is  described  in  details  elsewhere  (Salama  and  Ling,  1980).  X-  and  Y-cut  transducers  of  10 
and  2.25  MHz  were  used  for  the  generation  of  the  longitudinal  and  transverse  waves, 
respectively.  The  elastic  constant  are  calculated  using  the  relationship  Qj  =  p  Vjj2  where  p  is  the 
mass  density. 


RESULTS  AND  DISCUSSIONS 
Extruded  MMCs 

The  measured  Young's  and  shear  moduli  along  the  three  principal  axes  in  the  A1  7064  MMCs 
are  listed  in  Table  2  as  a  function  of  the  volume  fraction  of  SiC.  The  same  moduli  for  the  A1 
8091  MMCs  are  shown  in  Table  3.  Also  included  in  these  tables  are  the  elastic  moduli  predicted 
using  a  model  in  which  the  composites  are  assumed  to  consist  of  an  aluminum  matrix  and 

Table  2.  Calculated  and  measured  Young’s  and  shear  moduli  of  the  A1 7064 
MMC  specimens.  Moduli  are  shown  in  units  of  [GPa]. 


MMC 

Modulus 

Al  7Q64  +_ 
calc. 

rniic 

meas. 

Al  7064  + 

calc. 

15%  SiC 

meas. 

Ell 

90.0 

90.3 

96.8 

99.9 

E  2  2 

90.0 

90.5 

96.8 

97.2 

E33 

90.9 

91.4 

98.1 

100.5 

G  1  2 

34.0 

34.9 

36.6 

38.0 

G 1  3 

34.1 

35.6 

36.8 

39.6 

G23 

34.1 

35.3 

36.8 

38.4 

Table  3.  Calculated  and  measured  Young's  and  shear  moduli  of  the  A!  8091 
MMC  specimens.  Moduli  are  shown  in  units  of  [GPa], 


MMC 

Modulus 

M  809  L 

calc. 

t.-lfl2i_SiC 

mcas. 

A1  8091 
calc. 

+  15%  SiC 

mcas. 

Ell 

92.6 

91.3 

99.3 

96.8 

E22 

92.6 

91.8 

99.3 

96.7 

E33 

93.2 

92.6 

100.2 

98.3 

G 1 2 

35.7 

35.8 

38.4 

38.1 

G 1 3 

35.8 

36.2 

38.5 

38.7 

G23 

35.8 

36.1 

38.5 

38.5 

homogeneously  distributed  particles  of  spherical  shape.  The  panicles,  however,  are  arranged  so 
that  they  form  unidirectional  fibers  of  particle-rich  and  particle-free  aluminum.  This  arrangement 
is  chosen  to  model  the  microstructure  observed  in  the  extruded  composites  containing  spherical 
particles  where  the  particles  cluster  together  and  form  areas  with  high  particle  concentration  and 
islands  of  pure  matrix  material.  Due  to  the  extrusion,  these  areas  align  in  a  rod-like  shape  along 
the  extrusion  direction. 

From  Table  2  and  Table  3  one  can  see  that  the  anisotropy  for  predicted  Young's  moduli  between 
the  extrusion  direction  and  the  two  transverse  directions  is  of  the  same  order  (-1%)  as  that  of 
the  measured  moduli.  The  anisotropy  predicted  for  the  shear  moduli  is  negligible  and  the 
measured  shear  moduli  is  within  the  experimental  error.  However,  there  is  a  tendency  in  all 
measurements  indicating  that  the  shear  stiffness  is  lower  in  the  XiX2-plane  than  in  the  other  two 
planes.  This  agrees  with  the  difference  observed  in  the  calculations  of  Young's  moduli  where 
the  extrusion  direction  is  found  to  be  stiffer  than  the  two  other  directions. 

Pressed  MMCs 

The  measured  Young's  and  shear  moduli  along  the  principal  axes  in  the  A1  6061  MMCs 
together  with  the  elastic  moduli  predicted  by  the  model  assuming  dilute  concentrations  of 
homogeneously  distributed  particles  of  spherical  shape  are  plotted  in  Figs.  2  and  3  as  a  function 
of  the  volume  fraction  of  SiC-particles.  The  plots  clearly  shows  that  the  model  predicts  the 
overall  increase  of  the  elastic  moduli  relatively  well  for  specimens  containing  25%  SiC. 
However,  it  fails  for  the  40%  SiC-reinforced  specimen.  In  this  specimen,  the  values  predicted 
by  the  model  are  much  smaller  than  the  measured  ones.  The  disagreement  between  predicted 
and  measured  moduli  is  due  to  the  fact  that  the  model  assumes  a  dilute  concentration  of 
particles.  The  plots  in  Figs.  2  and  3  also  indicate  that  the  increase  of  the  moduli  deviates  from  a 
linear  relationship  towards  higher  values  of  volume  fractions  of  reinforcement.  This  behavior 
may  be  attributed  to  particle  interactions  when  the  mean-free-pathlength  between  them  is 
reduced  as  their  volume  fraction  is  increased. 

The  measured  elastic  anisotropies  found  between  the  directions  in  the  plate  and  that  normal  to 
the  plate  can  be  explained  by  the  metallurgical  investigations  where  it  is  observed  that  the  area  of 
the  particles  in  the  compression  direction  is  slightly  larger  in  the  compression  direction  than  that 
in  the  plane  of  the  plate.  This  suggests  that  the  particles  have  the  shape  of  spheroids  and  that 
they  tend  to  align  during  the  compression  such  that  their  c-axes  coincide  with  the  compression 
direction.  The  alignment  of  the  particles  results  in  a  higher  elastic  stiffness  in  the  plane  of  the 
plate  which  agrees  with  the  measured  elastic  properties.  The  metallurgical  investigations  also 
show  that  the  aspect  ratio  is  only  slightly  less  than  one.  According  to  calculations  by  Ledbetter 
and  Datta,  such  a  low  aspect  ratio  does  not  result  in  elastic  anisotropies  as  large  as  they  are 
found  in  the  investigated  MMCs  which  are  approximately  5%.  Calculations  assuming  an  aspect 
ratio  of  0.8  predict  an  anisotropy  of  the  order  of  1%  for  the  A1  6061  matrix  reinforced  with  25% 
SiC. 


Volume  percentage  of  SiC 


b 


Tabic  4.  Calculated  and  measure*]  Young's  and  shear  moduli  of  the  squeeze  cast  MMC 
specimens.  Moduli  arc  shown  in  units  of  IGPa], 


MMC 

Modulus 

A1  123  + 
calc. 

23%  AbCh 

meas. 

A!  13  + 
calc. 

19%  AbQa 
mcas. 

Al  1100  + 
calc. 

15%  AhCh 

mcas. 

Ell 

102.2 

101.9 

94.9 

97.1 

88.0 

86.4 

E  2  2 

102.2 

102.0 

94.9 

97.2 

88.0 

88.2 

E33 

99.6 

100.7 

92.1 

94.3 

85.0 

83.9 

Gl2 

38.8 

38.4 

35.9 

36.9 

33.5 

32.8 

Gl3 

37.4 

37.8 

34.4 

35.2 

31.9 

31.0 

G23 

37.4 

38.2 

34.4 

35.8 

31.9 

31.1 

Measurements  in  the  unreinforced  material  have  shown  that  the  texture  present  influences  the 
elastic  properties  significantly.  Here,  the  texture  accounts  for  a  4%  higher  Young’s  modulus  in 
the  compression  direction.  Thus,  it  seems  reasonable  to  consider  that  the  texture  is  responsible 
for  the  major  pan  of  the  anisotropy  observed  in  the  reinforced  specimens. 

Squeeze  Cast  MMCs 

In  the  squeeze  cast  specimens,  the  alumina  fibers  are  randomly  oriented  in  one  plane  whereas 
the  alumina  particles  are  also  homogeneously  distributed.  The  elastic  moduli  along  the  principal 
axes  of  the  composites  are  predicted  using  a  model  which  assumes  the  composites  to  be 
comprised  of  an  effective  matrix  consisting  of  the  pure  aluminum  and  the  homogeneously 
panicles.  The  effective  matrix  is  then  considered  to  be  reinforced  with  planar-randomly 
distributed  fibers.  The  predicted  as  well  as  the  measured  elastic  moduli  are  listed  in  Table  4.  The 
results  show  that  the  increase  in  the  elastic  moduli  due  to  the  reinforcement  is  reasonably  well 
predicted  by  the  model.  A  more  crucial  test  of  the  applicability  of  this  model  is  its  ability  to 
predict  the  elastic  anisotropy  between  the  directions  in  the  fiber-rich  plane  and  the  direction 
normal  to  the  plane.  Since  the  fiber  content  is  the  same  for  all  three  specimens,  the  anisotropy  is 
expected  to  be  more  pronounced  the  lower  the  stiffness  of  the  effective  matrix  i.e.  the  lower  the 
concentration  of  alumina  particles  is.  This  behavior  is  experimentally  observed  as  shown  by  the 
values  of  the  Young's  and  shear  moduli  in  Table  4.  The  difference  in  the  moduli  between  the  in¬ 
plane  directions  and  the  normal  direction  is  less  pronounced  the  higher  the  stiffness  of  the 
effective  matrix  becomes.  Quantitatively,  the  measured  anisotropies  are  well  predicted  by  the 
model.The  model,  however,  assumes  that  there  is  no  interaction  between  particles  and  fibers 
and  therefore  it  is  more  suitable  for  composites  with  low  particle  content. 
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Abstract.  Parameters  characterizing  the  elastic  nonlinearity  in  metal-matrix  composites 
tire  studied.  The  composites  consist  of  the  aluminum  alloys  K091  or  7064  containing 
silicon  carbide  particles  up  to  20%  volume  fraction.  Two  different  ultrasonic  measure¬ 
ments.  namely  the  acoustoelastic  effect  and  the  harmonic  generation,  are  used  for  the 
determination  of  acoustic  nonlinearity  parameters.  Their  dependence  on  the  content  of 
SiC  in  the  composite  is  investigated.  The  values  of  the  nonlinearity  parameter  arc  found 
to  decrease  with  increasing  volume  fraction  of  SiC-particles.  The  changes  are  explained 
in  terms  of  the  effects  of  SiC-particles  on  the  second  and  third  order  elastic  constants  of 
the  composites. 


Introduction 

Clastic  nonlinearity  is  responsible  for  the  deviation  ofji  -material's  stress-strain 
response  from  the  linear  relationship  represented  by  the  isothermal  form  of  the 
generalized  Hooke's  law.  This  law  can  be  written  as 

frii  =  C|jk|Kkl' 

where  <rt,  and  eu  are  the  stress  and  strain  tensors  respectively,  and  C\,u  is  the 
fourth  order  tensor  which  represents  the  second  order  elastic  constants.  The 
nonlinear  elastic  behavior  of  materials  can  be  determined  from  measurements 
of  the  stress  depcndancc  of  ultrasonic  velocities  as  well  as  the  distortion  of 
ultrasonic  waves  by  the  generation  of  higher  harmonics.  Consequently,  these 
effects  can  be  used  to  characterize  a  material’s  elastic  nonlinearity  nondestruc¬ 
tive^. 

In  recent  studies,  the  inllucnce  of  microstructurc  on  the  acoustic  nonlinear¬ 
ity  parameter  has  been  investigated  in  aluminum  alloys.  Razvi  cl  al.  HI  found  a 
linear  correlation  between  the  acoustic  nonlinearity  parameter  and  the  volume 
fraction  of  second  phase  precipitates  in  some  heat  treatable  aluminum  alloys. 
Another  correlation  was  found  by  Yost  et  al.  |2]  between  the  hardness  and  the 
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acoustic  nonlinearity  parameter  in  maraging  steel.  These  studies  directed  the 
attention  to  similar  investigations  of  the  nonlinearity  parameter  in  two  phase 
materials  such  as  metal-matrix  composites  (MMC).  These  composites  consist 
of  a  metallic  matrix  and  a  metallic  or  ceramic  reinforcement  which  can  he 
presented  in  the  form  of  particles,  short  or  long  libers.  The  properties  of  the 
composite  depend  on  those  of  the  constituent  phases  as  well  as  on  the  interac¬ 
tion  between  the  phases  in  the  composite.  The  interaction  is  determined  by  the 
geometrical  arrangement  of  the  reinforcement  as  well  as  the  bond  strength 
between  the  two  phases.  In  order  to  verify  the  desired  bulk  properties,  the 
nondestructive  characterization  of  linear  and  nonlinear  elastic  behavior  of  such 
composites  is  essential. 

In  this  investigation,  nonlinear  elastic  effects  were  determined  as  a  function 
of  the  SiC-particle  content  in  the  aluminum  alloys  8091  and  7064  by  measuring 
the  stress  dependancc  of  the  ultrasonic  velocity  and  the  harmonic  distortion  of 
ultrasonic  waves.  These  measurements  yield  the  effective  second  and  third 
order  clastic  constants  as  well  as  the  effective  acoustic  nonlinearity  parameter 
of  the  composite. 


Elastic  Nonlinearity  in  Solids 

All  elastic  properties  of  a  solid  can  be  derived  from  the  elastic  energy  represen¬ 
tation  of  that  solid.  In  general,  the  strain  energy  of  deformation  per  unit  vol¬ 
ume.  <t>.  can  be  expanded  as  a  power  series  in  the  elastic  strain,  e.  as 


*h(E)  +  C,iE,i  +  C uk|E||Eu  i"  “  (  iiklmnf'ifEk|Emn  ~  ■  •  •  •  t—l 

According  to  Brugger  [3],  <l>(,  is  the  strain  energy  per  unit  volume  in  the  unde- 
formed  stale  and  the  coefficients  C,,...  in  the  expansion  are  the  elastic  constants. 
In  a  linear  clastic  solid  there  arc  no  contributions  to  the  elastic  energy  from 
terms  higher  than  the  power  two  in  e.  All  materials,  however,  have  anharmonic 
elastic  potential  which  is  determined  by  the  ratio  of  the  third  and  the  second 
order  contributions  to  the  elastic  energy. 

In  an  isotropic  solid,  the  strain  energy  density  depends  only  on  the  invari¬ 
ants  / 1,  /:,  and  /x  of  the  strain  tensor  since  the  clastic  constants  are  invariant 
under  arbitrary  rotations.  Murnaghan  |4]  writes  for  an  isotropic  material  as 


<l>  =  <|)( 


A  +  2(jl  , 
+  <T  1 1  +  - Z -  /  T 


/  +  2»i 

3 


I]  -  :»»/,/;  +  M/v,  (3) 


where  cr  is  the  stress.  A  and  n  arc  the  second  order  Lamd  constants,  and  /,  m, 
and  n  arc  the  third  order  clastic  constants  in  Murnaghan  notation.  The  second 
and  third  order  elastic  constants  can  be  conveniently  obtained  from  ultrasonic 
measurements.  The  second  order  Lamd  constants  are  directly  related  to  ultra¬ 
sonic  velocities  as 


til 


li  =  i>V$  and  A  =  /»(Vf,  -  2V$). 


where  p  is  the  density  and  l's  and  Vy  are  the  shear  and  the  longitudinal  wave 
velocities,  respectively. 

A  formalism  which  allows  the  evaluation  of  the  third  order  elastic  constants 
/.  in.  and  it  from  the  strain  dependance  of  the  ultrasonic  velocity  was  developed 
by  Hughes  and  Kelly  |5).  In  an  uniaxially  strained  specimen,  the  velocities  of 
ultrasonic  waves  depend  linearly  on  the  elastic  strain'-.  Commonly,  three  wave 
modes  are  used:  the  two  linearly  polarized  shear  waves  and  the  longitudinal 
wave,  having  polarization  and  propagation  directions  as  indicated  in  Fig.  I.  The 
slopes  of  these  linear  relationships  normalized  with  respect  to  the  velocities  in 
the  strain-free  specimen  are  called  acoustoelastic  constants  and  related  to  the 
third  order  elastic  constants  /.  m.  and  n  as 


where  vy  are  the  ultrasonic  velocities  with  propagation  and  polarization  direc¬ 
tions  as  indicated  in  Fig.  I .  the  superscript  0  denotes  the  velocity  at  zero  strain, 
and  v  is  the  Poisson's  ratio. 

Another  method  for  the  determination  of  the  nonlinear  elastic  effect  is  the 
generation  of  higher  harmonics.  An  originally  sinusoidal  sound  wave  gets  dis¬ 
torted  while  propagating  through  a  nonlinear  material.  The  anharmonieilv  of 
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(he  interatomic  potential  function  leads  to  the  nonlinear  wave  equation  which 
can  he  written  for  a  longitudinal  wave  propagating  in  an  isotropic  solid  as  |6| 


P 


i)2u 

lit2 


-  (A 


2/  +  4  m 
A  +  2fx 


(6) 


where  n  is  the  displacement  and  /  and  v  are  the  time  and  space  variables, 
respectively.  This  nonlinear  wave  equation  differs  from  the  linear  wave  equa¬ 
tion  by  the  perturbation  term  on  the  light  hand  side  which  is  determined  by  (he 
ratio  of  the  third-  and  second-order  elastic  constants.  The  term  in  brackets  in 
I:q.  6  is  called  the  acoustic  nonlinearity  parameter  /j.  and  is  related  to  these 
elastic  constants  as 


-/i  =  3  +-  (2/  +  4 m)/(A  +•  2^)-  (7) 

It  is  interesting  to  see  from  I’q.  h  that  the  condition  for  linear  longitudinal  wave 
equation  is  not  the  third-order  elastic  constants  to  be  equal  to  zero,  but  the  ratio 
(2/  +  4w)/(  A  +  2/lx)  is  equal  to  -3. 

According  to  Thurston  and  Shapiro  |6],  the  solution  of  this  equation  for  u 
can  be  obtained,  using  a  perturbation  technique,  as 


n  =  A i  sintwr  -  kx)  +  A:  cos(2a ot  -  2 kx)  +  .  .  .  .  (8) 

Here,  m  is  the  frequency  of  the  ultrasonic  wave,  k  is  the  wave  number,  and  A  i  is 
the  amplitude  of  the  fundamental  wave.  The  generated  second  harmonic  will 
have  an  amplitude  A;  and  twice  the  fundamental  frequency.  The  amplitude  of 
the  second  harmonic  is  related  to  the  square  of  the  fundamental  amplitude  as 


A: 


=  jj  PAtA-.v. 


(9) 


where  .v  is  the  propagation  distance  of  the  wave. 

In  multi  phase  materials  such  as  metal-matrix  composites,  the  constituents 
may  have  very  different  elastic  and  mechanical  properties  which  result  in  the 
effective  macroscopic  properties  of  the  composite,  in  general,  these  effective 
properties  are  not  expected  to  follow  a  linear  law  of  mixture.  Cantrell  et  al.  |7| 
showed  that  the  effective  nonlinearity  parameter  of  multi  phase  materials  is  a 
nonlinear  function  of  the  volume  fraction  of  the  participating  phases.  Further¬ 
more,  different  phases  interact  with  each  other  at  the  interfaces  where  very 
high  residual  stresses  arc  likely  to  be  present.  These  are  either  compatibility 
stresses  due  to  differences  in  the  elastic-plastic  properties  of  the  constituent 
phases  or  thermal  stresses  arising  from  differences  in  the  coefficients  of  thermal 
expansion.  The  lattice  distortion  in  these  phases,  due  to  the  presence  of  resid¬ 
ual  stresses,  will  additionally  contribute  to  the  effective  acoustic  nonlinearity 
parameter  of  the  composite. 
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Table  I.  Chemical  composition  of  the  investigated  aluminum  alloxs  in  wt. 


Alloying  elements 


Alloy 

Si 

Fc 

Cu 

Mg 

Zr 

l.i  Zn  Ci  Co 

Al-7064 

0.05 

0.10 

2.00 

2.30 

0.20 

—  7.10  0.12  0.22 

Al-8091 

0.02 

0.01 

1 .90 

0.80 

0.20 

2.70  —  — 

Experimental 

Specimens 

In  the  present  investigation,  two  metal-matrix  composites  with  matrices  of  the 
heat  treatable  aluminum  alloys  Al-8091  and  AI-7064  are  examined.  Their  chemi¬ 
cal  compositions  are  shown  in  Table  1 .  The  aluminum  matrices  are  reinforced 
with  globular  silicon  carbide  (SiC)  particles  up  to  20%  volume  fraction  of  sizes 
ranging  between  I  and  5  pm.  The  specimens  were  received  as  extruded  rods  of 
I  inch  in  diameter.  Due  to  the  manufacturing  process,  the  particles  are  aligned 
along  the  extrusion  direction,  while  in  the  plane  perpendicular  to  the  extrusion 
direction  the  particles  are  randomly  distributed.  Micrographs  of  these  two 
planes  arc  shown  in  Fig.  2  for  the  composite  consisting  of  Al-8091  and  10%  SiC- 
particles.  For  harmonic  generation  measurements,  the  specimens  are  cut  to  a 
length  of  two  inches.  Opposite  faces  are  machined  parallel  and  then  polished  to 
optical  flatness.  After  measuring  the  acoustic  nonlinearity  parameter,  the  sides 
of  the  specimens  were  milled  to  a  square  shaped  cross  section  in  order  to 
measure  acoustoelastic  constants. 


Ultrasonic  Velocity  Measurements 

In  order  to  determine  the  third  order  elastic  constants,  the  strain  dependance  of 
the  ultrasonic  velocity  needs  to  be  measured.  This  effect  is  typically  in  the 
order  of  one  part  in  It)4  and  requires  very  precise  time  of  flight  measurements. 
The  system  used  in  this  investigation  utilizes  the  pulse-echo  overlap  method  for 
the  determination  of  the  time  of  flight  as  described  in  (8).  In  Fig.  3.  a  block 
diagram  of  the  system  components  is  shown.  Using  this  system,  an  absolute 
time  resolution  of  ±200  ps  can  be  measured.  During  the  time  of  flight  measure¬ 
ments.  compressive  strains  are  applied  to  the  specimen  along  the  extrusion 
direction  and  are  varied  systematically  in  the  elastic  regime.  The  lime  of  flight 
data  arc  acquired  by  a  computer  and  ultrasonic  velocities  are  calculated  as  a 
function  of  the  applied  strain.  Corrections  are  made  for  the  change  in  the  lateral 
dimension  of  the  strained  specimen  which  is  due  to  the  Poisson's  expansion. 
The  linear  relationship  between  the  ultrasonic  velocity  and  the  elastic  strain  is 
then  evaluated  by  a  least  squares  lit  algorithm.  The  slope  of  the  curve  is  nor¬ 
malized  with  respect  to  the  velocity  of  the  unstrained  specimen  and  is  called  the 
acoustoclastic  constant.  The  acoustoelastic  constants  of  the  three  wave  modes 
shown  in  Fig.  1  are  used  in  the  calculation  of  the  third  order  elastic  constants 
using  Fqs.  5. 
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Kir.  2.  Particle  distribution  in  extruded  MMC  specimens:  a  extrusion  direction,  b  transverse  direc¬ 
tion 


Absolute  Amplitude  Measurements 

For  the  direct  measurements  of  the  acoustic  nonlinearity  parameter  the  abso¬ 
lute  amplitudes  of  the  fundamental  and  the  second  harmonic  waves  are  deter¬ 
mined.  This  is  achieved  using  a  capacitive  detector  system  which  is  described 
in  detail  elsewhere  19].  The  block  diagram  in  Fig.  4  shows  the  experimental 
setup  schematically.  For  the  excitation  of  the  fundamental  longitudinal  wave. 
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Fig.  3.  Block  diagram  of  the  ultrasonic  time  of  Might  measurement  setup 


TEKTRONIX 
2465  B 
Oscilloscope 


HPX590  B 
Spectrum  Analyzer 


Fig.  4.  Block  diagram  of  the 
harmonic  generation  setup 
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Fig.  5.  Correlation  between 
fundamental  and  second  harmonic 
amplitude 


an  electric  10  MHz  signal  is  applied  to  a  narrow  band  undamped  lithium  niobate 
transducer  which  is  bonded  to  the  surface  of  the  specimen.  The  longitudinal 
wave  which  has  a  bandwidth  of  200  kHz  and  a  duration  of  5  /j.s  is  propagated 
along  the  extrusion  direction.  The  superposition  of  the  fundamental  wave  and 
its  second  harmonic  is  picked  up  by  the  capacitive  transducer  on  the  opposite 
surface.  During  the  measurements,  the  received  signal  is  amplified  and  the 
amplitudes  of  its  frequency  components  are  analyzed.  These  measurements  are 
repeated  for  different  driving  amplitudes.  The  amplitude  of  the  second  har¬ 
monic  is  plotted  vs.  the  square  of  the  fundamental  amplitude  in  Fig.  5  for  tw  o  of 
the  specimens  examined  in  this  investigation.  The  linear  relationship  between 
the  two  quantities  is  predicted  by  Eq.  9.  The  slope  of  a  linear  least  squares  tit 
through  the  measured  points  is  then  used  to  evaluate  the  acoustic  nonlineariH 
parameter  /j. 

Results  and  Discussion 

For  both  sets  of  composites,  the  second  and  third  order  elastic  constants  are 
calculated  using  measured  ultrasonic  velocities  and  their  strain  derivatives, 
respectively.  Their  values  are  listed  in  Table  2.  It  has  to  be  noted  that  the  third 


Table  2.  Second  and  third  order  elastic  constants  of  NlMCs 


SOEC  |GPa|  TO  EC  |GPa] 


A. 

li 

r 

1 

w 

n 

A 1-8091 

44.9 

4 1.0 

SO.  2 

0.30 

218 

5^8 

-4.45 

-  IOC?  SiC 

40.3 

35.8 

91.7 

0.28 

-  185 

-  365 

-  454 

-15%  SiC 

46.0 

38.1 

96.9 

0.27 

-  162 

-  460 

-423 

AI-7064 

59.3 

27.4 

73.4 

0.34 

-  324 

-397 

-403 

-15%  SiC 

57.0 

34.7 

91.0 

0.31 

-215 

-  397 

-4S4 

-20%  SiC 

57.6 

38.1 

99  1 

0.30 

■214 

-397 

-461 
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Table  3.  Nonlinearity  parameters  ol'  metal-matrix  composites  as 
calculated  from  elastic  constants  and  measured  from  harmonic  gen¬ 
eration 


fi  ftq.  12) 

Rel.  change 

0  Imeas.l 

Rel.  change 

A 1-8091 

15.2 

09? 

10.2 

09? 

-  109?  SiC 

12.6 

-  \T-, 

8.5 

179? 

-  IV?  SiC 

11.5 

249? 

"5 

-269? 

AI-7064 

16.6 

0r? 

8.6 

()9? 

-  IV?  SiC 

13.0 

._  ncj 

6  9 

-  20'? 

-209?  SiC 

12.1 

-  17rr 

>9 

-  3i  9; 

order  elastic  constants  of  the  composites  are  determined  using  formulas  w  hich 
were  derived  for  isotropic  materials,  although  the  micrographs  reveal  some 
anisotropic  particle  distribution .  This  may  affect  the  absolute  value  of  the  third 
order  elastic  constants  determined  in  these  composites,  ll  is  assumed,  how¬ 
ever.  that  changes  in  these  values  due  to  different  particle  content  is  not  af¬ 
fected. 

The  data  in  Table  2  indicate  a  considerable  increase  in  the  magnitudes  of 
the  Young's  modulus  £  and  the  shear  modulus  /u.  with  the  SiC-particle  content. 
At  the  same  time,  the  Poisson's  ratio  u  is  reduced.  The  experimental  error  in 
these  constants  is  estimated  to  be  within  27c .  One  can  also  see  from  Table  2  that 
the  changes  in  the  third  order  elastic  constants  are  less  pronounced.  The  inac¬ 
curacy  in  these  values  is  15 7c  for  /.  5 7c  for  m  and  3rf  for  n.  Only  in  the  case  of 
the  constant  /.  a  clear  influence  of  reinforcement  is  observed.  Its  magnitude 
decreases  when  SiC-particles  are  added  to  the  matrix.  The  value  of  the  con¬ 
stant  /,  however,  saturates  at  higher  particle  content  tn  the  composite.  For  the 
constants  m  and  n.  no  significant  trend  can  be  observed. 

The  values  of  the  second  and  third  order  elastic  constants  are  used  to 
calculate  the  nonlinearity  parameters  according  to  Eq.  7  and  their  values  are 
listed  in  Table  3.  The  inaccuracy  in  the  calculated  nonlinearity  parameter  is 
estimated  to  be  10 7c.  Also  included  in  Table  3  are  the  values  of  the  nonlinearity 
parameter  determined  directly  from  the  harmonic  generation  experiments. 
Here,  the  inaccuracy  is  estimated  to  be  of  the  order  of  15rf.  From  the  data  in 
Table  3.  one  can  see  that  both  calculated  and  experimentally  measured  nonlin¬ 
earity  parameters  decrease  considerably  with  increasing  particle  content. 

Table  3  also  indicates  that  the  experimentally  measured  values  of  the 
acoustic  nonlinearity  parameter  are  about  507c  smaller  than  those  calculated. 
This  difference  may  be  caused  by  the  anisotropy  in  the  elastic  properties  due  to 
the  particle  alignment  as  well  as  the  texture  in  the  aluminum  matrix  |10|.  The 
values  of  the  directly  measured  nonlinearity  parameter  are  obtained  using  lon¬ 
gitudinal  waves  which  propagate  along  the  extrusion  direction  while  those 
calculated  are  obtained  using  second  and  third  order  elastic  constants  mea¬ 
sured  in  the  transverse  direction.  Because  of  the  strong  dependence  of  the 
nonlinearity  parameter  on  the  crystallographic  orientation  of  single  crystals, 
texture  is  likely  to  introduce  significant  anisotropy  in  this  parameter.  Further¬ 
more.  inaccurate  absolute  values  of  the  third  order  elastic  constants,  due  to  the 
use  of  formulas  which  are  only  valid  for  isotropic  materials,  may  contribute  to 
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Fig.  6.  Relative  change  of 
measured  and  calculated 
nonlinearity  parameters  as  a 
function  of  SiC-particle  content 


the  observed  difference.  Good  agreement,  however,  is  found  in  the  relative 
changes  of  the  measured  and  calculated  nonlinearity  parameters.  As  can  be 
seen  from  Fig.  6.  the  relative  change  in  the  nonlinearity  parameter  with  the  SiC 
content  indicates  a  linear  relationship  between  the  two  quantities.  The  relative 
change  is  in  the  order  of  3 O'T  over  the  range  of  200  volume  fraction  of  SiC- 
particles  added  to  the  Al-alloys. 

An  explanation  for  the  decreasing  values  of  the  nonlinearity  parameter  in 
MMCs  may  be  found  by  separately  analyzing  the  contributions  of  the  second 
and  the  third  order  elastic  constants  shown  in  Fq.  ?.  The  second  order  elastic- 
constants  appear  in  the  denominator  of  Eq.  7  and.  accordingly,  their  increase 
with  particle  content  will  yield  a  decrease  in  the  nonlinearity  parameter.  The 
increase  in  these  constants  is  determined  by  the  elastic  moduli  of  the  constitu¬ 
ents  of  the  composite  as  well  as  the  interface  between  reinforcement  and  ma¬ 
trix.  Theoretically,  the  upper  and  lower  limits  for  the  elastic  moduli  of  two 
phase  materials  are  determined  by  the  isostrain  and  isostress  conditions  |llj. 
The  third  order  constants  in  the  numerator  of  Eq.  7  decrease  in  magnitude  as  a 
function  of  particle  content  and.  thus,  their  contributions  also  lower  the  values 
of  the  nonlinearity  parameter  in  the  composites.  However,  it  appears  that  the 
contributions  of  the  second  order  elastic  constants  dominate  the  magnitude  of 
the  effective  nonlinearity  parameter  in  the  composite,  especially  at  higher 
amounts  of  reinforcement. 

This  analysis  can  also  be  applied  to  the  results  obtained  by  Razvi  et  al.  on 
AI-7075  alloys  containing  different  amounts  of  second  phase  precipitates  w  hich 
are  shown  in  Fig.  7.  Here,  the  nonlinearity  parameter  was  measured  using 
harmonic  generation  and  found  to  increase  with  the  volume  fraction  of  second 
phase.  Also  included  in  Fig.  7  are  the  results  obtained  on  the  AI-7064  metal- 
matrix  composites.  Unlike  the  SiC-particle  reinforced  AI-7064  MMCs.  the 
change  in  the  second  order  elastic  constants  of  AI-7075.  due  to  the  presence  of 
precipitates,  is  found  to  be  very  small  1 1|.  On  the  other  hand,  the  presence  of 
the  second  phase  precipitates  is  responsible  for  a  high  degree  of  distortion  in 
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Fig.  7.  Comparison  of  measured 
nonlinearity  parameters  in 
precipitation  hardened  AI-7075  and 
AI-70M  based  MMC  specimens 


the  matrix,  which  is  the  effective  hardening  mechanism  in  such  alloys.  Conse¬ 
quently,  the  denominator  of  Eq.  9  remains  approximately  constant.  Unfortu¬ 
nately.  the  effect  of  second  phase  precipitates  on  the  third  order  clastic  con¬ 
stants  of  AI-7075  was  not  determined.  However,  measurements  of 
acoustoelastic  constants  on  steel  |I2]  and  heat  treatable  aluminum  alloys  113] 
indicate  that  these  constants  increase  as  the  amounts  of  second  phase  precipi¬ 
tates  in  these  alloys  are  increased.  These  results,  nevertheless,  show  that  the 
third  order  elastic  constants  in  these  alloys  are  likely  to'increase  as  a  function 
of  second  phase  precipitates,  and.  in  turn,  contribute  to  the  increase  of  the 
acoustic  nonlinearity  parameter  which  is  opposite  to  that  observed  in  metal- 
matrix  composites. 


Conclusions 

The  second  and  third  order  elastic  constants  as  well  as  the  acoustic  nonlinearity 
parameter  have  been  determined  in  metal-matrix  composites  consisting  of  SiC- 
particles  up  to  20rf  volume  fraction  and  AI-7064  or  -‘.i  8091.  The  results  show 
that  the  acoustic  nonlinearity  parameter  is  significantly  influenced  by  the  vol¬ 
ume  fraction  of  the  second  phase  reinforcement  present  in  the  composite.  The 
third  order  elastic  constants  m  and  n  remain  r.ichanged  while  the  constant  /  as 
well  as  the  nonlinearity  parameter  decrease  as  the  volume  fraction  of  SiC- 
particles  is  increased.  This  behavior  can  be  explained  in  terms  of  the  increase  ot 
the  effective  second  order  elastic  constants  and  the  decrease  of  the  effective 
third  order  elastic  constants  with  :  einforccmcnt.  This  argument  can  also  be 
used  to  explain  the  increase  in  the  acoustic  nonlinearity  parameter  as  a  function 
of  second  phase  precipitates  observed  by  Razvi  et  al.  in  the  aluminum  alloy 
7075. 
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